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1.0  BACKGROUND 


The  computer  code  described  in  this  report  -  The 
APL  RCS/Statistics  Code  -  represents  years  of  de¬ 
velopment  by  E.  Shotland,  J.  W.  Follin,  Jr.,  and  F. 
C.  Paddison  (APL)  and  J.  W.  Crispin,  Jr.,  A.  L. 
Maffett,  and  K.  M.  Siegel  (the  University  of 
Michigan).  The  results  are  in  the  form  of  first-  and 
higher-order  statistics  for  the  radar  cross  section 
(RCS)  of  a  complex  target  as  a  function  of  fre¬ 
quency,  polarization,  aspect  angle,  and  constituent 
parameters  of  the  target.  The  statistics  are  developed 
from  the  premise  that  a  complex  target  can  be  broken 
down  into  a  set  of  simple  scatterers  whose  RCS  can 
be  more  manageably  calculated  from  available  tech¬ 
niques  (thus,  the  program  can  be  conveniently  up¬ 
dated  to  include  the  most  accurate  techniques  cur¬ 
rently  at  hand)  and  whose  relative  phases  are  ran¬ 
domly  (i.e.,  uniformly)  and  independently  dis¬ 
tributed. 

Methods  for  determining  both  first-order  RCS  sta¬ 
tistics  for  a  complex  target  and  also  simple  com¬ 
ponent  scatterer  RCS  data  were  developed  by  mem¬ 
bers  of  Siegel’s  Radiation  Laboratory  at  the  Univer¬ 
sity  of  Michigan,  principally  Crispin  and  Maffett. 
The  extensions  of  these  methods  to  include  the 
(higher-order)  lobe  and  glint  statistics  were  developed 
by  Shotland  and  Follin  under  the  sponsorship  of  the 
Advanced  Research  Projects  Agency  in  an  extension 
of  the  Advanced  ALBIS  Program.  The  (higher- 
order)  beta  statistics  were  developed  by  Follin  during 
the  Cruise  Missiles  Observables  Program.  All  of 
these  statistical  methods  were  organized  into  a  coher¬ 
ent  program  by  Maffett,  and  the  result  was  coded  by 
H.  W.  Klimach.  Paddison  served  as  Program 
Manager  for  the  APL  efforts. 

The  techniques  upon  which  the  simple  scatterer 
RCS  calculations  are  based  are,  to  a  large  extent,  the 
methods  of  geometrical  and  physical  optics  with  their 
attendant  assumptions  and  restrictions  on  body 
smoothness  and  size  relative  to  illuminating  wave¬ 
length.  Included  in  the  code,  however,  are  other 
methods  that  are  appropriate,  for  example,  to  travel¬ 
ing-wave  phenomena  and  wedges  (with  straight 
edges). 

For  regions  where  the  ratio  of  body  length  to  wave¬ 
length  may  become  important  (particularly  if  radar 
absorbing  materials  have  been  used  to  reduce  large 
RCS  contributions),  diffraction  effects  have  not  been 


taken  into  account,  but  several  methods  are  being 
evaluated  to  determine  which  would  be  most  appro¬ 
priate  to  the  nature  of  the  APL  code.  The  methods 
include  the  geometrical  and  physical  theories  of 
diffraction,1'2  the  equivalent  current  technique,3  and 
the  numerical  electromagnetic  code.4  The 
geometrical  theory  is  based  on  the  tracing  of  rays  and 
can  be  described  by  currents  induced  on  the  illum¬ 
inated  portion  of  the  structure.  The  physical  theory 
adds  currents  induced  in  the  shadowed  area  by  dif¬ 
fraction  of  the  incident  field.  The  equivalent  current 
technique  adds  edge  currents  on  assumed  filamentary 
edges.  The  numerical  electromagnetic  code  is  an 
integral-equation  surface-current  determination  in 
which  the  structure  is  broken  up  into  small  cells,  and 
the  induced  charge  and  current  on  each  piece  are 
calculated  from  the  incident  field  and  from  the 
charge  and  current  of  each  other  cell.  It  is  planned  to 
include  one  or  more  of  these  methods  in  the  APL 
code  to  account  for  various  diffraction  effects  caused 
by  some  of  the  scattering  components  of  a  complex 
target. 

Bistatic  RCS  computational  capabilities  are  not  yet 
fully  incorporated  into  the  code.  The  complexity  of 
bistatic  RCS  estimation  over  monstatic  is  greatly  in¬ 
creased  by  its  additional  aspect  variables,  so  the  sta¬ 
tistical  approach  gains  still  more  importance  as  a  re¬ 
ducer  of  data  bulk.  The  bistatic  capability  will  be  in¬ 
cluded  in  the  APL  code  as  soon  as  possible. 

New  coordinate  system  arrangements  have  been  in¬ 
troduced  recently  to  simplify  the  execution  of 
conical-aspect  views  of  a  target  (see  the  Appendix). 


Vl.  15.  Keller,  "(.ieometrical  theory  of  Oil liaction,"  J.  Opl. 
W.  1»i..  52.  I  If)  I  50(1962). 

I’.  Va.  iJIrimises.  "Met  nod  ol  Edge  Waves  in  the  Physical 
I  henry  ol  Diffraction.'  /cd-to.  Sue.  Radio,  1-24.5  (1962); 
translation  by  U  S.  An  force  foreign  technology  Division 

U97I), 

I  I  .  Knoll  and  I.  15.  A.  Senior.  "Comparison  of  three  Higlt- 
I -requeue)  Diffraction  I ecliniqnes,”  Proc.  ULt  62,  1468- 
1474(1974). 

*Ci.  .1.  Burke  ami  A.  J.  Poggio,  \umenvul  Lleclruiiiaxnetic 
Colic  -  Method  oj  Moments.  Naval  Electronics  Systems  Com¬ 
mand  I  D  1 16  ( 1977). 
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2.0  GENERAL  DISCUSSION  OF  RCS  DATA 


Experimental  and  predicted  RCS  data  usually  are 
in  the  form  of  a  received  signal  referenced  to  the 
transmitted  signal,  versus  aspect  angle.  Typically, 
there  are  separate  tabulations  for  each  frequency  of 
interest.  The  radars  used  for  aircraft,  missile,  and 
space-object  detection,  tracking,  and  identification 
are  usually  at  wavelengths  that  are  small  in  compari¬ 
son  to  the  dimensions  of  the  targets.  Usually  the 
targets  are  complex  shapes.  At  some  aspects  such  as 
broadside,  the  RCS  of  a  target  may  be  many  orders 
of  magnitude  greater  than  at  other  aspects,  may  be 
dominated  by  a  single  component,  and  hence  can  be 
characterized  simply.  However,  at  most  aspects  there 
is  a  complicated  interference  pattern  with  maxima 
spaced  roughly  \/L  apart  in  ang.e  and  c/L  apart  in 
frequency  (where  X  is  the  wavelength,  c  is  the  velocity 
of  light,  and  L  is  a  typical  length  parameter).  If,  for 
example,  L  is  10  meters  (i.e.  a  small  missile),  roughly 
109  data  points  are  needed  to  cover  X  band  for  all 
polarizations.  A  statistical  description  is  needed  to 
reduce  the  data  storage  and  retrieval  problem. 


The  magnitude  of  RCS  data  is  not  the  only  reason 
for  an  interest  in  a  statistical  description.  The  main 
users  of  these  data  are  the  hardware  designer,  the 
radar  system  designer,  and  the  performance  analyst. 
The  latter  two  have  long  searched  for  a  statistical  des¬ 
cription  of  the  amplitude  scintillation  of  RCS.  The 
literature  contains  many  attempted  fits  of  experi¬ 
mental  data  with  single-parameter  statistical  des¬ 
criptors. 

The  parameters  of  interest  in  the  vicinity  of  any 
one  aspect  are  the  mean  RCS  and  its  probability  dis¬ 
tribution,  the  lobe  widths  in  angle  and  frequency, 
and  the  mean  centroid  and  its  variances. 

Over  the  years,  the  technique  of  breaking  down  a 
complex  target  into  a  finite  set  of  component  scatter¬ 
ed  has  been  evolved  together  with  statistical  descrip¬ 
tions  of  lobe  structure  and  the  appropriate  descriptor 
of  amplitude  scintillation,  the  two-parameter  beta 
distribution  function.  The  computer  code  described 
in  this  report  embodies  the  aforementioned  predic¬ 
tive  RCS  and  its  statistical  descriptors. 


3.0  THE  CAPABILITIES  OF  THE  CODE  FROM  THE  VIEWPOINT 

OF  THE  RCS  ANALYST 


The  APL  RCS/Statistics  Code,  designed  for  the 
PDP  11/60  computer,  has  been  checked  out  with 
several  missile  configurations.  Its  computational 
capabilities  fall  into  two  categories,  either  of  which, 
with  various  options,  can  be  executed  from  a  master 
code.  They  are: 

1 .  First-order  statistics  for  RCS,  and 

2.  Lobe  and  higher-order  statistics  for  RCS. 

At  this  stage,  the  first-  and  higher-order  statistics  are 
limited  to  monostatic  situations.  However,  the  lobe 
statistics  are  coded  in  such  a  way  as  to  be  applicable 
to  bistatic  situations  as  soon  as  that  capability  is  in¬ 
troduced. 


The  entire  development  is  based  on  the  premise 
that  a  complex  scattering  body  can  be  broken  down 
into  a  finite  set  of  component  scatterers  (usually 
chosen  to  be  elementary  scatterers,  some  of  which 
may  correspond  to  simple  geometric  shapes).'  If  the 
contributions  from  this  set  of  scatterers  are  combined 
with  proper  phases,  the  result  is  an  estimate  of  rela¬ 
tive-phase  RCS;  this  option  is  available  in  the  APL 
code.  Figure  1  is  an  example  of  relative-phase  RCS. 


.1 .  w  Crispin.  Jr.,  and  A.  I  Mat  ten .  ‘‘I- mi  mat  ion  o  l  RC  S  toi 
Simple  Shapes,'*  "I  Miinanoii  >»|  RC  S  tor  Complex 
Shapes.*'  /V'/t  .  Nil  (IWPL 
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Figure  1  —  Example  of  calculated  monostatic  relative- 
phase  RCS  1/24th  scale  Convair  990. 


If  it  is  assumed  that  phases  among  scattering  com¬ 
ponents  are  distributed  independently  and  randomly 
(i.e.,  uniformly,  over  the  interval  0  to  27t),  a 
hierarchy  of  statistics  can  be  developed.  Under  first- 
order  statistics  fall  the  usual  mean  RCS  and  root- 
mean-square  (rms)  spread  about  the  mean  (see  Fig. 
2).  Under  higher-order  statistics  fall  higher  moments 
from  which  skewness  and  kurtosis  may  be  derived  lo 
show  the  appropriate  beta  distribution  characteristics 
for  the  RCS  data  under  examination.'  "*  (These 
higher-order  statistics  can  also  be  used  to  analyze 


AJ.  W.  t  ollin,  Jr..  Statistical  /'ra/wwi  «/  Ki  S.  IHt  API  y\l 
81-115(1981). 

J.  W.  Follin,  Jr.,  and  A.  1  .  Malleti.  "RC  S  Scimillaiions  and 
Their  Statistical  Description.”  in  I’m c.  Second  Animal  lai- 
tical  Air  Surveillance  Control  Con/.,  Rome  Air  l)c\clopnicni 
Cenler  (1981). 

*J.  W.  Follin.  Jr.,  F.  C  Paddison,  and  A  I  Mallen.  "  I  lie  RC  S 
of  Two  Cruise  Missiles  and  their  Statistical  Description." 
DARPA  8th  Strategic  Space  Symposium  ( 1982). 
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Figure  2  —  Example  of  calculated  monostatic  mean  RCS 
with  rms  spread  (including  lobe  statistics);  1/25th  scale 
Convair  990. 


measured  RCS  data;  this  feature  is  not  part  of  the 
APL  code.)  Various  properties  of  the  autocorrelation 
of  RCS  can  be  used  to  describe  lobe  statistics  (in  both 
frequency  and  aspect).1'  Lobe  characteristics  are  illus¬ 
trated  in  Fig.  2  at  observation  angle  intervals  of  30°. 
The  terms  "major”  and  "minor”  refer  to  half  the 
long  and  short  axes  of  the  ellipse;  they  give  the  mean 
lobewidth  at  that  aspect  in  both  the  horizonal  plu,.c 
(minor)  and  the  normal  plane  (major),  “^ngle” 
refers  to  the  tilt  of  the  lobe  ellipse  with  the 
observation  axis.  N,  (GHz)  refers  to  the  number  of 
lobes  of  RCS  per  gigahertz  of  frequency. 

At  this  writing,  the  component  scatterer  analyses 
rest  largely  upon  physical  and  geometrical  optics 
methods.  However,  the  code  could  be  expanded 
readily  to  include  diffraction  treatments  and  even  the 
more  exact  integral-equation  analyses  exemplified  by 
the  numerical  electromagnetic  code.J 

The  program  as  written  and  described  here  was  in¬ 
tended  for  targets  of  modest  complexity  at  observing 
frequencies  in  the  centimeter  region.  We  plan  to  re¬ 
code  the  program  for  a  more  capable  computer  in  the 


i  .  Shoiland,  Stan.sinal  \nalvsis  n/  Nadar  lancet  S anniiaimn. 
IHt)  API.  PC.  955  (1967). 
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near  future  to  allow  the  handling  of  more  complex 
targets  up  to  the  mid-millimeter  range  of 
frequencies. 

The  effects  of  a  surface  treatment  with  radar  ab¬ 
sorbing  material  can  be  taken  into  account  through 
the  uniform  application  of  a  reflection  coefficient  to 
the  component  on  which  the  material  is  placed.  No 
account  is  taken  of  the  possibility  of  volumetric  scat¬ 
tering  for  either  radar  absorbing  material  or  dielec¬ 
tric  components. 

The  major  emphasis  of  the  code  is  on  a  statistical 
examination  of  the  RCS  of  a  target.  Considering  the 
number  of  variables  upon  which  RCS  depends,  a 
very  voluminous  collection  of  data  can  be  amassed 
for  a  single  target,  even  when  domains  of  variables 
are  rather  severely  restricted.  A  statistical  examina¬ 
tion  of  RCS  offers  a  condensation  of  facts  and,  thus, 
an  economy  of  output  not  available  with  a  determin¬ 
istic  treatment  of  similar  scope. 

All  of  the  RCS  statistics  are  essentially  unaffected 
by  small  changes  in  target  configuration  or  frequency 
that  may  drastically  affect  the  lobe  structure  of  deter¬ 
ministic  treatments  (especially  for  frequencies  of  1 
GHz  and  above). 

Last  but  not  least,  statistical  RCS  results  mesh  well 
with  the  requirements  of  detection  and  tracking  sys¬ 
tems  for  statistical  formulations. 


3.1  Inputs 

The  inputs  usually  required  for  the  code  are  wave¬ 
length;  polarization  (linear,  but  extendable  to  ellip¬ 
tical);.  plane  of  observation;  range  of  observation 
angle  in  the  observation  plane;  pitch,  roll,  and  yaw  of 
the  target;  and  tilt  of  the  target  rotation  plane  (to 
make  comparisons  with  RATSCAT  measurements). 

For  a  particular  complex  target,  a  suitable  mathe¬ 
matical  model  must  be  derived  in  the  form  of  a  set  of 
component  scatterers  matched  to  the  target  features 
with  the  help  of  drawings,  pictures,  and  artists'  con¬ 
ceptions.  The  set  of  component  scatterers  must,  of 
course,  be  limited  to  the  types  presently  available  to 
the  APL  code  (some  20  as  of  this  writing). 

Computations  of  RCS  are  performed  in  coordinate 
systems  generic  to  the  various  components  and  so 
must  be  transformed  to  the  coordinate  systems  in 
which  the  statistics  are  formulated. 


3.2  Outputs 

The  first-order  and  lobe  statistics  have  already- 
been  illustrated  in  Fig.  2.  Note  that  for  the  latter 
case,  for  any  line  through  the  lobe  ellipse  center,  the 


distance  from  the  center  to  the  ellipse  intercept  gives 
the  lobewidth  for  that  direction  of  observation.  For 
example,  at  the  30°  observation  angle,  the  ellipse  has 
a  major  axis  oriented  82°  from  the  observation  axis, 
indicating  thin  lobes  in  the  horizontal  plane  of  the 
target  and  fat  lobes  in  the  perpendicular  plane. 

The  output  for  the  higher-order  statistics  may  be 
either  printed  or  graphical.  Let  us  consider  the  latter. 
One  has  available  skewness  versus  width  or  kurtosis 
versus  width  (width,  skewness,  and  kurtosis  are  de¬ 
fined  in  Section  4  of  this  guioe).  Skewness  and 
kurtosis  are  plotted  in  Figs.  3  and  4,  respectively,  ver- 


Width 


Figure  3  —  Calculated  width-skewness  statistics  for  the 
RCS  scintillation  of  a  1/25th  scale  Convair  990 

sus  width  for  a  typical  target  at  10°  observation  angle 
intervals.  It  should  be  noted  that  the  data  presented 
in  Figs.  3  and  4  may  be  based  either  on  experimental 
data  or  on  theoretical  calculations.  The  code  can  pro¬ 
vide  the  theoretical  values.  The  points  are  plotted 
within  boundaries  that  denote  the  limits  for  beta-dis¬ 
tributed  component-RCS  values.  Locations  denoting 
normal,  lognormal,  chi-square,  and  Rayleigh  distri- 
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Figure  4  —  Calculated  width-kurtosis  statistics  for  the  RCS 
scintillation  of  a  1/25th  scale  Convair  990. 


for  the  beta  distribution  given  by 


0(x;a,b) 


T(0  +  b ) 

r(a)T(6) 


x“-' 


(1  -  x)b  1 


where  T(a)  denotes  the  gamma  function  of  an  argu¬ 
ment  a.  Thus,  the  component  scatterers  that  con¬ 
tribute  to  the  RCS  at  100°  from  nose-on  produce, 
under  the  radom-phase  assumption,  an  RCS  that  is 
distributed  according  to  the  above  beta  distribution. 
[Note:  The  skewness  can  be  calculated  to  be  0.88, 
which  agrees  very  well  with  the  skewness  for  100° 
found  at  width  0.7  in  Fig.  3.  The  beta  distribution  for 
the  above  values  of  a  and  b  is  skewed  to  the  left,  with 
a  long  tail  to  the  right.  ] 

The  beta-distributed  quantity  x  is  actually  a/  a max, 
where  a  is  the  RCS  for  some  fixed  aspect,  frequency, 
and  polarization.  So,  for  example,  since  omax  (100°) 
=  1.14  m2  and  since  the  mean  of  @(x;a,b)  is 
a/(a  +  b)  =  0.2,  the  mean  RCS  at  100°  is  a  = 
0.2(1.14)  =  0.23  m2  or  -6.4  dBsm  (dBsm  is  decibels 
relative  to  a  square  meter).  Note  that  this  value 
agrees  with  the  mean  RCS  at  100°,  which  can  be  read 
on  Fig.  2.  The  beta  distribution  of  RCS  at  a  fixed 
observation  angle  arises  from,  and  depends  on,  the 
assumption  that  the  phases  among  components  are 
randomly  distributed  (i.e.,  uniformly,  between  0  and 
2w)  for  every  observation  angle  value. 


butions  all  fall  outside  of  these  boundaries,  for  both 
width/skewness  space  and  width/kurtosis  space. 

For  example,  consider  Fig.  4  where  the  observation 
direction  100°  produces  in  width/kurtosis  space  the 
point  (0.7, 0.7).  These  statistics  yield  the  parameter 
values 

a  =  0.96,  b  =  3.82 


3.3  Glint 

Illustrations  of  the  output  of  the  program  predict¬ 
ing  glint  are  not  included  in  this  writing;  that  portion 
of  the  program  has  yet  to  be  validated.  However,  a 
discussion  of  the  theory  used  is  given  in  Section 
4.6. 
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4.0  CALCULATIONS 


4. 1  Relative-Phase  RCS 

Consider  a  class  of  complex  targets,  each  of  whose 
members  can  be  broken  down  into  a  finite  set  of 
component  scatterers.  Each  component  scatterer  is 
assumed  to  be  in  the  far  field  of  the  transmitter  so 
that  it  is  illuminated  by  a  plane  wave;  it  is  also  in  the 
far  field  of  the  receiver  so  that  the  observed  field  has 
a  plane  wavefront.  The  major  dimension  of  a  com¬ 
ponent  is  assumed  to  be  large  with  respect  to  the 
wavelength  \.  Polarizations  of  the  transmitter  and 
receiver  are  taken  to  be  horizontal  (H)  or  vertical  (V) 
in  this  discussion;  however,  they  could  be  circular  or 
elliptical.  The  possible  combinations  are  HH,  HV, 
VH,  and  VV:  only  the  first  and  last  are  considered 
here. 

Let  the  signal  scattered  by  the  ith  component  be 
denoted  by 

a,  exp  J  =  v  -T7  (1) 

where  the  amplitude  a,  (which  is  usually  taken  to  be 
v  ct„  where  a,  is  the  RCS  of  the  fth  component)  and 
the  phase  may  depend  on  aspect  (for  both  trans¬ 
mitter  and  receiver,  or  on  their  bistatic  separation  |3), 
frequency,  polarization,  and  the  constituent  paiam- 
eters  of  the  component.  The  scattered  amplitude  is 
then 

A  =  £  Vo,  exp(/t/',).  (2) 


The  relative-phase  RCS,  a,  for  the  complex  target  is 

S  , 

a  =  AA*  =  j  £  Vo,  exp(Ji, )  | 

1  =  1 

\ 

=  £  A  fok  expOV-,  )exp(  -j^k )  .  (3) 

i.k  =  t 

4.2  Random-Phase  RCS 


where  the  asterisk  means  complex  conjugate  and 
<  •  •  •  >  indicates  an  average  over  all  phases.  The  de¬ 
viation  from  the  mean  is 

F  =  AA'  E\a } 

=  2  £  a,ak  cos  (4,  -  ^  ) 


=  2  £  Vo, a k  cos(i/  -  \j/k)  (5) 


whence  the  variance  is 


M;  =  <F:>  =  2S;  =  2  £  o,ok 


~  ( £  °.)  -  S  (6) 

'=1  t=l 

4.3  Higher-Order  Statistics 

If  the  triple  and  quadruple  sums  are  written  as 
S,  =  £  0,0,0,  (7) 

,■  k<  , 

and 

■Sa  =  £  o,oka,a„,,  (8) 

I  •  ''  I  <  rti 

the  third  and  fourth  central  moments  are,  respective¬ 
ly/ 

Mi  =  <F’>  =  12S,  (9) 


For  the  random-phase  approximation,  the  mean 
RCS,  Ela],  over  all  random  phases,  is  anc* 

v,  (  *£[o]  1  =  {AA‘)  =  Eo„  (4)  Mj  =  <F*  >  =  6S,:  +  12S,S,  +  132S4  .  (10) 
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In  order  to  examine  the  nature  of  the  distributions 
governing  RCS  for  a  complex  target,  we  have 
modified  Pearson’s  method  somewhat.  Instead  of 
examining  target  RCS  behavior  in  skewness-kurtosis 
space,  we  have  set  up  two  spaces,  one  a  width-skew¬ 
ness  space,  the  other  a  width-kur;  -is  space,  where 
width  w  is  defined  as 


where 


y,'  +  2 

n  —  — - -  and  7 1  is  as  above, 

2w  —  7,' 


<  w<  1  ,  lim  7,  =  6  ; 


w  =  n2/v , :  .  (II) 

In  fact,  to  simplify  the  analysis  somewhat,  we  have 
defined  a  modified  skewness  as 

■>i  =  Pi/g;'’,  •  (12) 

Kurtosis  is  defined  in  the  standard  manner  as 


Lower  boundary: 


7,  =  —1.5,  0  <  w  <  -  : 

2 


Left  boundary: 

part  1 :  w  =  0 ,  -  1.5  <7,  <  0  ; 


From  the  nature  of  the  sums  Slt  i  =  1,  2,  3,  4, 
above,  it  can  be  shown  that  the  allowable  regions  in 
width-skewness  (actually,  read  modified  skewness) 
and  width-kurtosis  spaces  are  almost  triangular  and 
are  restricted  as  follows  (see  Figs.  3  and  4): 


part  2:  w  = 


t(i  +  2) 

U  +  l)5 
6 l(t  +  4) 
(I  +  2)T  ' 


(as  above) , 


4.3.1  Width-Skewness 

Right  boundary: 


7/  =  4 w  -  2,  -  <  w  <  I 

r,  2 


Lower  boundary: 


The  location  of  points  within  the  restricted  regions 
will  vary  with  the  aspect  at  which  a  complex  target  is 
illuminated  and  viewed.  The  corresponding  beta  dis¬ 
tribution  parameters  may  be  calculated  from  those 
coordinates.  The  beta  distribution  is  defined  as 

r(u  +  b) 

(J(x)  =  — — ~  X'  '(1  -  x)  ,  (20) 

r(a)l'(b) 


7 1  =  0,  0  <  w  < 

Left  boundary: 
at  +  2) 

w  =  - -  , 

( t  +  n- 

,  =  2  at  +  3) 

7|  ”  (/  +  i)(/  +  2)  ’ 


0  <  t  <  00 . 


4.3.2  Width-Kurtosis 

Right  boundary: 


6 n:  1  3  1 

- -  (  H’  +  ,  -  —  —  ■  ) 

(n  +  l)(n  +  2)  tv  w  n  n 


Os.vsl, 

where  r(u)  is  the  gamma  function  of  argument  w, 
and  .v  =  0/ a,,,,,, .  Then,  for  example,  the  parameters  a 
and  b  can  be  determined  from  width-kurtosis  space. 
The  result  can  be  used  to  find  skewness,  which  can  be 
checked  against  the  location  of  the  aspect  point  in 
width-skewness  space.  The  equations  are 


a(a  +  b  +  1) 

2 (b  -  a) 
a(u  +  b  +  2) 

6  (a  +  b  +  1) 

(a  +  b  +  2  )(a  +  b  +  3) 

/a  b  1  \ 

(  + - -  3)  . 

V)  a  a  +  b  +  1  ' 
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The  values  of  a  and  b  may  be  determined  from  any 
pair  of  these  equations,  but  the  first  two  should  be 
used  for  experimental  data.  If  one  takes,  for 
instance,  the  aspect  of  100°  (Fig.  2),  one  can  compute 
a  =  0.96  and  b  =  3.82. 

In  order  to  check  the  beta  distribution  theory 
against  measured  data,  a  slightly  different  viewpoint 
must  be  taken  because  the  output  of  a  measurement 
exercise  is  usually  in  the  form  of  RCS  versus  aspect. 
If  the  RCS  is  sufficiently  oscillatory  in  the  chosen 
aspect  interval  (usually  10°,  but  ranging  from  5  to 
30°,  depending  on  frequency),  the  random  phases 
have  had  the  opportunity  to  range  over  most  of  the 
intervals  of  definition,  and  so  the  RCS  in  the  aspect 
interval  should  be  distributed  in  essentially  the  same 
way  that  the  RCS  at  a  fixed  aspect  was  owing  to  the 
random  phase  fluctuations. 


4.4  Modified  Statistical  Distribution 

An  extension  of  the  statistical  theory  just  presented 
has  been  developed  where  the  requirement  of  con¬ 
stant  component  RCS  is  relaxed  and  an  individual 
variance  and  third  and  fourth  moments  about  the 
mean  are  assigned.  The  physical  regions  of  the  statis¬ 
tical  space  resulting  in  this  modification  will  be  pub¬ 
lished  shortly. 


where  k  is  the  wavenumber.  The  results  obtained  are 
the  number,  N. ,  of  lobes  of  RCS  per  hertz  of  fre¬ 
quency: 

N,  =  (i  ii,  !((',-  0(1  +  cos  .i) 

cS  v  “ 

+  (S  -  £, )  sin  ,3|;  )  ,  (23) 

where  c  is  the  speed  of  light  and  S  =  is  from  Eq. 
6.  Also  obtained  are  the  number,  N  ,  of  lobes  of  RCS 
per  radian  in  traverse  relative  to  a  rotation  axis 
oriented  at  an  angle  6  with  the  £  axis  in  the  plane: 

/V  =  ( £  A.  d.d.  )  ,  r,s  =1,2,  (24) 

Ao  '  v  ' 

where 

d,  =  sin  f> ,  d:  =  -cos  6,  (25) 

and 

A  =  E  o  nk  ((£.  -$.)(!  +  cos  d) 


4.5  RCS  Lobe  Structure 

If,  instead  of  taking  the  relative  phases  to  be  ran¬ 
dom  variables,  one  examines  their  variation  (and 
hence  the  variation  of  relative-phase  RCS)  as  a  func¬ 
tion  of  receiver  aspect  (in  a  bistatic  situation,  either 
aspect  or  bistatic  angle),  it  is  possible  to  obtain 
scintillation  statistics.  The  results  are  derived  from 
the  RCS  correlation  function,  using  the  zero  crossing 
theorem  of  Rice."’  They  represent  an  extension  of 
results  for  the  monostatic  case  obtained  by 
Shotland.1' 

Suppose  the  target  is  at  the  origin  of  a  £,  r },  f  coor¬ 
dinate  system  whose  £(' plane  contains  the  transmitter 
and  receiver.  Then  the  phase  term  ^  of  Eq.  1  is  given 
by 

=  k  |  f,  ( 1  +  cos  d)  +  £,  sin  dl  .  (22) 


-  (v  -  j'.  )  sin  dp' ' 


x  l(>)  ->?.)(  1  +  cos d)T  •  (26) 

The  numbers  ,V.  and  .V  ,  considered  as  a  function  of 
n  for  a  fixed  aspect  angle,  form  ellipses.  They  can 
also  be  converted  into  local  lobewidths  of  frequency 
and  aspect,  respectively,  by  considering  their  recipro¬ 
cals.  For  any  line  through  the  ellipse  center,  the  dis¬ 
tance  from  the  center  to  the  ellipse  intercept  gives  the 
lobe  rale  for  that  direction  of  rotation  axis. 


4.6  Clint 

The  apparent  electrical  location  of  the  target  re¬ 
flection  depends  on  the  type  of  radar  seeker  used  to 
measure  it.  For  the  relative-phase  return  (Eq.  2)  with 
the  receiver  location  at  X,  Y,7.;  X,  Y<  < Z;  we  have 


S.  ().  Rice,  Maihemaiical  AimKms  * »t  Random  Noinc.”  m  ('j\,  ft  k(X  —  c  )  X  —  c 

Selected  Papers  nn  \tnse  and  Snnfiusin  Pn nesses,  \cK*mi  -  -  -  (kr  )  =  -  -  =  k—  .  ,  (27) 

Wax.  cd.,  Dmer  Publications.  New  York  1 1954).  dX  dx  f  f 
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and 


and  the  variance  of  the  glint  is 


drh 

dz 


k-^-jd±k 

r, 


(28)  <*e  >  =  E  o,o,  (£,  +  £,  -  2ft2  .  (32) 

2c  jKj 


The  apparent  source  location  is 


£  =  rtan  1 


sdA/dXx 

\dA/dz) 


So  far  we  have  discussed  only  one  component  of 
glint,  but  it  is  clear  that  expressions  similar  to  Eqs.  30 
and  32  are  needed  for  -q.  In  addition,  the  variance  is  a 
tensor  operator  so  that  an  off-diagonal  term  is 
needed. 

Defining 


Zika,  £,  expQj, ) 
£  ika,  expOV, ) 


Bmn  =  £  o,o,  U,  +  £,  -  2£)4~ 
•*j 


£  a, a, l  exp(/'(^  -  4*,)) 


E  a, a,  e\p(J ( \p,  -  ift)) 


X  (v,  +  V,  -  2v)m+n-2 , 
and 

B  =  E  Bn,n  aman , 

m,n 


m,n  =1,2  (33) 


(34) 


a, a, l  cos  ty,  -  yj/j) , 


where  ax  =  cos  a,  a2  =  sin  a,  the  standard  deviation 
(29)  0f  gijnt  in  the  a  direction  is 


where  the  real  part  is  taken  in  the  last  step  since  the 
RF  signal  is  actually  proportional  to  the  real  part  of 
Eq.  2. 

The  division  is  normally  performed  by  an  auto¬ 
matic  gain  control  (AGC)  circuit  with  a  time  constant 
so  that  the  numerator  and  denominator  are  averaged 
independently.  The  result  in  the  random-phase  ap¬ 
proximation  is 


S„  = 


2 E(o) 


VB. 


(35) 


Further  in  analogy  to  the  derivation  of  Eqs.  23  and 
26,  we  define 

Cmn  =  E  o.o, +  £,  -  2£)4~'"“" 


I  =  ^  .  (30) 

La, 

With  a  fast  AGC,  the  fluctuations  of  Eq.  29  become 
larger  but  cannot  be  expressed  in  closed  form.  In 
addition,  the  result  depends  on  the  dynamic  range  of 
the  circuits.  Note  that  in  a  bistatic  geometry  the  glint 
is  measured  relative  to  the  receiver  and  is  indepen¬ 
dent  of  bistatic  angle  except  for  the  effects  of 
shadowing  of  components. 

Subtracting  Eq.  30  from  Eq.  29,  we  have  the  devi¬ 
ation  from  the  mean  for  the  slow  AGC  case. 


A£  =  ~  E  a,fl,(ft  -  ft  cos  M  -  i pj) 


=  '  E  aia,(ft  +  ft  -  2ft cos  -  tj),  (31) 

^  i<j  * 


X  (tj,  +  77,  -  2rj i)"""-2 


x  [(£  -  r,)(l  +  cosff)  +  (£,  -  ft)  sin#]2  , 

(36) 


Dmnn  =  E  a-a,(ft  +  ft  -  2ft4  '”-” 


X  (1,,  +Vj-  2rj)m+n~2 


x  [(ft  -  ft)(l  +  cos/3)  -  (r,  -  f,)  sin /3]4~f-J 
x  (U  -  Tjy)(l  +  cos^r*-2,  (37) 

C  =  E  C«»a»«. ,  (38) 
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and 


D=T,Y.DmnrsamandA-  (39) 

m,  n  r,s 

Then  the  number  of  glint  lobes  in  the  a  direction  per 
cycle  is 


N. 


"J  2cS 


(40) 


and  the  number  of  glint  lobes  in  the  a  direction  per 
radian  of  rotation  about  the  axis  oriented  in  the  6 
direction  is 


AT ,  = - VD. 

2  XS„ 


(41) 


Note  that  Eqs.  40  and  41  are  correct  even  though 
the  assumption  of  a  slow  AG C  has  been  used;  only 
the  amount  of  the  displacement  is  in  error. 


5.0  USER’S  GUIDE  TO  THE  CODE  FROM  THE  VIEWPOINT 
OF  THE  COMPUTER  ANALYST 


5.1  Introduction 

The  APL  RCS/Statistics  Code  is  a  computer  soft¬ 
ware  package  capable  of  computing  and  displaying 
the  collective  cross  section  from  a  group  of  simple 
components  using  physical  optics  as  the  main  model¬ 
ing  basis.  In  addition  to  average  cross  section,  it  is 
also  able  to  compute  and  display  relative-phase  cross 
section,  RCS  lobe  and  glint  lobe  characteristics,  and 
higher-order  beta  statistics  of  the  collective  object  as 
a  function  of  observation  angle,  in  either  printed  or 
graphical  form.  Input  for  the  program  is  taken  from 
a  data  base  in  batch  mode  process,  and  output  results 
are  displayed  on  a  graphics  terminal,  plotter,  or  line 
printer.  The  RCS  analytical  basis  is  given  in  Refs.  1 1 
and  12. 

The  program  is  currently  running  on  a  PDP  1 1  /60 
computer  using  the  RSX-1 1M  V3.2  operating  system. 
The  graphical  peripherals  available  on  this  system  are 
the  Tektronix  terminal  and  the  ZETA  plotter.  The 
program  is  executed  by  issuing  the  computer  com¬ 
mand  RUN  RCS.  It  responds  by  prompting  the  user 
for  the  name  of  the  input  data  base  file.  The  data 


11  A.  L.  Maffett,  Bistatic  Lobe  Statistics,  JHU/APL  M80-7  (Nov 
1980). 

i2A.  L.  Maffett,  Linear  Transformations  Relating  Theoretical  and 
Experimental  Coordinate  Systems  Used  in  RCS  Determina¬ 
tion,  JHU/APL  QM-80-049  (M80-4)  (Feb  1980). 


base  name  must  then  be  entered,  followed  by  any 
optional  switches.  The  available  options  are; 

/G  enable  graphical  debugging 

/GC  enable  graphical  debugging 

(combined  components) 

/P  enable  printed  debugging 

/Z  direct  graphics  output  to  plotter 

(default  is  terminal) 

The  program  then  performs  the  functions  defined  in 
the  data  base  and  returns  to  computer  command 
mode. 

This  section  of  the  report  provides  the  user  with  a 
basic  understanding  of  the  program  along  with  a 
comprehensive  reference  manual  for  actual  program 
operations.  The  internal  structure  and  program  capa¬ 
bilities  are  presented  in  Section  5.2.  Section  5.3  pre¬ 
sents  the  display  options,  including  the  available  de¬ 
bugging  output.  General  features  of  the  input  data 
base,  such  as  looping  and  command  priorities,  are 
described  in  Section  5.4.  The  actual  RCS  commands 
are  presented  in  detail  in  Section  5.5. 


5.2  RCS  Program  Operations 

The  program  is  designed  to  compute  the  combined 
RCS  (called  random-phase  RCS)  of  several  simple 
components  as  a  function  of  observation  angle, 
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wavelength,  and  polarization.  In  addition  to  average 
cross  section,  the  program  also  computes  relative- 
phase  cross  section,  RCS  lobes,  average  glint,  and 
glint  lobes.  Finally,  it  may  also  be  used  to  compute 
higher-order  beta  statistics. 

A  complex  object  to  be  analyzed  in  terms  of  RCS 
must  first  be  separated  into  its  constituent  com¬ 
ponents.  The  basic  components  available  to  the  pro¬ 
gram,  such  as  ogives,  cylinders,  corner  reflectors, 
etc.,  are  described  in  Section  5.5.  Each  component 
must,  at  a  minimum,  be  specified  by  its  name,  physi¬ 
cal  dimensions,  and  orientation  in  the  radar  coordi¬ 
nate  system.  The  location  of  each  component  (or 
scattering  center)  in  the  radar  coordinate  system  must 
also  be  specified  if  lobe  or  glint  output  is  required. 
The  geometry  describing  the  scattering  centers  of 
each  available  component  has  been  specified  in  Ref. 
13. 

The  program  may  be  divided  functionally  into  four 
areas:  input,  cross-section  computation,  statistics 
and  lobe/glint  computations,  and  output  display. 
The  input  function,  which  is  actually  in  continuous 
operation  throughout  the  computations,  is  described 
in  detail  in  Section  5.4  in  the  context  of  a  data  base. 
In  summary,  the  input  must  be  ordered  to  supply  the 
program  with  global  parameters,  followed  by  local 
component  parameters,  and  ending  with  display  pa¬ 
rameters.  Since  knowledge  of  the  input  function  or 
output  display  function  is  not  necessary  to  under¬ 
stand  the  program’s  computational  algorithms,  they 
will  not  be  discussed  further  in  this  section. 

The  program  first  zeros  an  array  of  accumulators 
and  then  adds  in  the  cross-section  contributions  of 
each  component  read  on  input.  This  process  con¬ 
tinues  until  the  advent  of  an  output  display  instruc¬ 
tion,  which  causes  the  accumulated  results  to  be  dis¬ 
played  as  specified.  The  RCS  program  operates  com¬ 
pletely  on  one  component  at  a  time  by  computing  the 
component  RCS  at  a  fixed  wavelength,  polarization, 
and  radar  phi  angle  and  then  varying  the  radar  theta 
angle  (phi  and  theta  are  the  radar  polar  coordinates). 
Each  cross-section  value  as  a  function  of  theta  is 
stored  in  a  separate  accumulator.  The  range  and 
number  of  theta  values  desired,  specified  on  input, 
must  not  require  more  than  181  accumulators  (a  pro¬ 
gram  limitation).  This  limitation  may  be  circum¬ 
vented  in  some  instances  as  described  in  Section  5.3 
and  defined  in  Section  5.5. 

In  addition  to  accumulating  average  cross-section 
responses,  the  program  also  saves  the  square  of  the 
values  of  the  responses  in  order  to  compute  and  dis- 


nA.  L.  Maffctt,  Component  Scanerme  (  enters,  I H L  API  M H I 
I  (7  Jan  1981). 


play  the  standard  deviation  of  the  overall  response  at 
each  theta  increment.  If  relative-phase  rather  than 
average  cross-section  output  is  requested,  it  will  use 
the  two  accumulator  vectors  to  save  the  complex 
components  of  cross  section  rather  than  the  detected 
sum  and  sum  squared  values.  All  results  generated 
thus  far  may  be  stored  in  computer  memory  since 
they  are  relatively  few  in  number.  The  lobe  and/or 
glint  options,  however,  require  the  temporary  storage 
of  large  volumes  of  intermediate  data.  Those  data, 
which  consist  of  the  cross  section  and  three-dimens¬ 
ional  position  of  each  subreflector  in  the  object  for 
each  theta  angle  of  interest,  are  stored  on  a  disk  file 
for  later  processing. 

The  following  paragraphs  summarize  the  RCS  pro¬ 
gram  procedure  for  computing  cross  section  for  a 
single  component.  Cross-section  calculations  do  not 
actually  begin  until  a  component  description  has 
been  read  from  the  input  stream. 

The  program  begins  to  process  a  component  by 
computing  a  transformation  matrix  from  component 
space  through  object  space  to  radar  space.  The  com- 
ponent-to-object  transformation  matrix  (TRANS) 
will  have  been  read  previously  as  input.  The  object- 
to-radar  transformation  matrix  will  have  been  com¬ 
puted  from  previously  read  values  of  radar  platform 
pitch,  roll,  yaw,  and  tilt  (PRYT). 

The  program  then  enters  the  THETA  loop  to  com¬ 
pute  the  cross  section  at  each  required  theta  angle.  In 
the  loop,  it  determines  if  the  current  radar  theta  and 
phi  angles  are  valid  for  that  component.  If  valid,  they 
are  transformed  into  the  component  coordinate  sys¬ 
tem.  At  that  point,  calculations  that  correspond  to 
the  specific  component  of  interest  are  performed. 

Although  they  differ  in  detail,  all  groups  of  com¬ 
ponent  calculations  (i.e.,  subroutines)  have  several 
things  in  common.  They  must  first  collect  their  ap¬ 
propriate  inputs  and  then  determine  if  a  valid 
response  is  possible  from  the  current  parameter  con¬ 
figuration.  Each  component  subroutine  must  then 
compute  its  RCS  as  a  function  of  geometry,  radar 
characteristics,  and  observation  angle.  If  lobe  output 
is  also  required,  the  component  subroutines  must 
also  compute  the  location  of  the  reflecting  source  or 
sources  for  multiple  reflectors.  Component  position 
information,  after  being  transformed  to  the  radar  co¬ 
ordinate  system,  is  stored  in  a  file  for  future  lobe 
and/or  glint  processing. 

When  the  component  cross  section  has  been  com¬ 
puted,  it  is  added  to  the  accumulator  bin  appropriate 
to  the  current  theta  angle.  The  powers  of  these  results 
are  computed  and  stored  in  another  file  if  optional 
beta  statistics  arc  desired.  Beta  statistics  are  com¬ 
puted  in  the  output  display  section  of  the  program 
using  the  relatively  simple  equations  from  Ref.  7.  In 
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any  case,  the  observation  angle  is  incremented  by 
repeating  the  THETA  loop.  After  the  component 
cross  section  has  been  computed  for  all  required 
values  of  theta,  program  execution  returns  to  the  in¬ 
put  section  to  obtain  new  component  information 
from  the  data  base,  and  the  process  is  repeated. 

Basic  RCS  computations  pause  when  a  display  in¬ 
struction  is  read  on  input.  This  causes  the  program  to 
display  its  cross-section  results  (described  in  Section 
5.3)  and  then  begin  lobe  and/or  glint  calculations  if 
required.  The  remainder  of  this  section  highlights  the 
general  procedure  for  computing  lobe  and  glint  char¬ 
acteristics.  A  complete  analytical  description  is  given 
in  Refs.  14and  15. 

The  RCS  lobe  calculations  begin  by  reading  tem¬ 
porary  results  stored  by  the  main  program.  Each 
component’s  contribution  to  RCS  at  the  observation 
angle  of  interest  (specified  by  DELTA,  a  subset  of 
THETA)  is  read  by  the  lobe  subprogram  along  with 
its  position  in  space.  The  first  parameters  computed 
from  these  data  are  a  simple  cross-section  sum  and 
standard  deviation.  Cross-section  products  of  com¬ 
ponent  cross  section  and  position  are  then  computed 
and  summed.  The  square  root  of  this  result  divided 
by  the  standard  deviation  and  the  speed  of  light 
yields  the  number  of  lobes  of  radar  cross  section  per 
hertz  of  frequency.  A  different  set  of  summed  cross 
products  is  used  to  compute  lobe  width. 

The  glint  calculations  begin  by  computing  average 
moments  for  the  components  read  in  the  lobe  sec¬ 
tion.  Sums  of  cross  products  using  component  cross 
section,  position,  average  moments,  and  an  intersect 
plane  defined  by  the  angle  alpha  are  used  to  compute 
the  glint  standard  deviation.  The  average  moment 
arm  length,  along  with  standard  deviation,  as  a  func¬ 
tion  of  observation  angle  (theta)  is  the  primary  out¬ 
put  from  the  glint  section. 

Glint  lobes  are  computed  from  the  same  set  of  data 
used  to  compute  glint  moments.  In  this  case,  an  ex¬ 
panded  set  of  summed  cross  products  is  computed 
and  used  in  conjunction  with  glint  standard  deviation 
to  yield  the  number  of  glint  lobes  per  hertz  of  fre¬ 
quency  with  respect  to  observation  angles  delta  and 
alpha.  An  even  greater  expanded  set  of  summed  cross 
products  is  then  computed  and  combined  with  glint 
standard  deviation  to  yield  glint  lobe  width. 

The  RCS  lobe  and  glint  lobe  computations  are  re¬ 
peated  for  each  required  observation  angle  defined 


.1.  W.  Crispin  ami  K  M  Siegel.  Methods  oj  Radar  (  Hiss  Set  mm 
Analysis,  Academic  Press,  New  York  (1968), 

1  A.  I  .  Maffctl,  lieometric  and  Physical  Optus  / ormolus  anil  Or 
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JHU/APL  MHO-3  (28  Jan  1980) 


by  the  angle  delta;  in  contrast,  RCS  and  glint 
moments  are  defined  by  the  observation  angle  theta. 
When  computations  for  all  values  of  delta  have  been 
completed,  the  results  are  displayed,  and  execution 
control  is  returned  to  the  main  program. 


§.3  Display  Options 

The  real  power  of  the  RCS  program  lies  in  its 
ability  to  display  cross-section  output  results  in  a 
complete  and  elegant  manner  in  both  printed  and 
graphical  form.  This  section  discusses  the  standard 
display  options  available  in  the  program  followed  by 
additional  debugging  output.  The  standard  display 
options  are  RCS,  RCS  lobes,  glint,  glint  lobes,  and 
higher-order  statistics. 

The  commands  that  initiate  the  display  function 
are  GRAPH  and  PRINT  (defined  in  Section  5.5).  As 
the  names  imply,  GRAPH  produces  graphical  output 
on  a  terminal  or  a  plotter,  and  PRINT  produces 
printed  output  on  a  terminal  or  a  line  printer.  Since 
the  results  displayed  by  these  commands  arc  general¬ 
ly  the  same,  only  differences  in  the  data  presented 
will  be  described  in  conjunction  with  the  overall  data 
available  for  display.  Program  results  are  not  altered 
because  of  their  being  displayed.  This  allows  multiple 
display  commands  to  be  included  in  the  data  base  in 
order  to  view  either  accumulated  results  or  new 
results  if  program  accumulators  are  cleared  (see  the 
CLEAR  directive). 

The  basic  program  shows  monostatic  full-scale 
RCS  as  a  function  of  observation  angle.  Both  aver¬ 
age  cross  section  (shown  with  the  optional  plus  or 
minus  standard  deviation)  and  relative-phase  cross 
section  may  be  displayed  in  square  meters  on  a 
logarithmic  scale  (i.e.,  in  decibels  relative  to  a  square 
meter).  Relative-phase  cross  section  may  also  be 
shown  in  square  meters  on  a  linear  scale.  Graphical 
RCS  output  also  lists  certain  program  parameters, 
such  as  frequency,  polarization,  and  percent 
shadowed,  on  the  edge  of  the  graph.  The  output 
from  this  section  is  always  displayed  as  part  of  the 
GRAPH  or  PRINT  operation. 

The  RCS  lobe  output  is  displayed  on  top  of  the 
basic  cross-section  output  in  graph  mode  and  follow¬ 
ing  it  in  print  mode  if  LOBE  is  active.  The  lobe  out¬ 
put  shows  the  number  of  RCS  lobes  per  hertz  of  fre¬ 
quency  and  the  number  of  lobes  per  spatial  increment 
around  a  circle.  In  graph  mode,  the  spatial  lobes  are 
displayed  inversely  to  yield  lobe  width.  Samples  of 
lobe  width  are  computed  around  a  circle  and  plotted 
as  an  ellipse  centered  at  the  observation  angle  of  in¬ 
terest.  The  major  and  minor  axes  of  the  ellipse  along 
with  the  orientation  angle  may  then  be  found  and  dis- 
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played.  Since  these  values  are  searched  out  rather 
than  computed,  lobe  width  must  be  sampled  often 
enough  (from  4  to  180  times  as  defined  by 
DELTA(4))  to  converge  on  the  required  accuracy. 

If  GLINT  is  active,  monostatic  glint  is  displayed  as 
a  function  of  all  the  independent  parameters  men¬ 
tioned  thus  far  and  also  'he  angle  alpha,  which  de¬ 
fines  an  observation  plane  that  cuts  through  the  lobe 
structu  -e.  Since  a  complete  graph  or  print  list  show¬ 
ing  glint  results  is  produced  for  each  value  of  alpha, 
only  certain  values  ,  namely  0  and  90°,  are  used.  The 
varying  independent  parameter  on  a  glint  graph  or 
print  list  is  the  observation  angle  theta.  The  depend¬ 
ent  parameter  represents  the  glint  moment  arm  dis¬ 
placement  in  degrees;  values  of  standard  deviation 
are  included  as  an  option. 

Monostatic  glint  lobe  output  is  displayed  in  a  man¬ 
ner  completely  analogous  to  cross-section  lobe  out¬ 
put.  It,  therefore,  similarly  shows  the  number  of 
lobes  of  glint  per  hertz  of  frequency  and  the  number 
of  lobes  per  spatial  angle  at  samples  evenly  spaced 
around  a  circle.  The  amount  of  glint  lobe  data  output 
is,  as  with  the  RCS  lobes,  completely  defined  by  the 
parameter  array  DELTA. 

The  final  standard  display  option  presents  higher- 
order  beta  statistics  computed  from  the  basic  RCS 
results.  The  type  and  amount  of  output  produced  by 
this  option  are  very  different  between  graph  and 
print  modes.  The  two  plots  generated  in  graph  mode 
show  skewness  and  kurtosis  as  a  function  of  beta 
width.  Distribution  acceptance  boundaries  are  drawn 
on  each  plot,  followed  by  the  actual  computed  data 
values,  which  are  annotated  with  their  corresponding 
observation  angle.  Many  additional  statistics  are 
computed  during  print  mode.  The  results  are  dis¬ 
played  within  the  normal  RCS  output  listing  rather 
than  at  the  end. 

Three  mutually  exclusive  debugging  display 
options  are  available  in  the  RCS  program.  All  three 
are  intended  to  provide  partial  results  corresponding 
to  each  separate  component  rather  than  the  com¬ 
bined  object.  When  the  cross-section  contribution  of 
a  component  is  computed,  it  is  immediately  displayed 
on  either  a  printer  (option  /P)  or  a  graphics  output 
device  (option  /G).  After  it  is  displayed,  the 
accumulator  vectors  are  reset  to  zero,  and  the  next 
component  is  computed  and  displayed  on  a  new 
listing  or  graph.  This  process  continues  until  the 
entire  data  base  has  been  read.  Prints  or  graphs  for 
which  there  is  no  component  contribution  are  not 
displayed.  Furthermore,  component  display  may  be 
aborted  in  graph  mode  by  entering  a  control  -C' 
from  the  terminal  keyboard  for  each  component  not 
desired.  The  final  debugging  option,  /GC,  instructs 


the  program  to  combine  all  component  curves  onto 
one  graph. 

Debugging  display  options  may  not  be  used  in  con¬ 
junction  with  lobe,  glint,  or  higher-order  statistics 
output.  Any  PRINT  or  GRAPH  commands  found  in 
the  data  base  are  simply  ignored  by  the  program  in 
this  mode. 


5.4  Data  Base  Description 

AH  inputs  to  the  RCS  program  must  be  suplied 
from  within  a  data  base.  A  data  base  is  a  file  that 
may  be  created  or  modified  by  an  editor.  It  contains 
all  the  instructions  necessary  to  control  the  program. 
Each  line  in  the  data  base,  except  comment  and  title 
lines,  constitutes  a  single,  complete  instruction  to  the 
program.  Each  instruction  line  begins  with  either 
four  or  five  characters,  defining  the  directive  type, 
followed  by  a  space  and  up  to  nine  numeric  parame¬ 
ters.  All  directives  (i.e.,  instructions)  recognized  by 
the  program  and  their  functions  are  defined  in  detail 
in  Section  5.5. 

In  this  section  we  will  discuss  directive  interactions 
and  usage  protocols,  the  concepts  of  global  and  local 
directives  followed  by  their  relative  ordering  within 
the  data  base,  and  then  looping  capabilities  within 
the  data  base. 

All  valid  directives  in  the  RCS  data  base  fall  into 
one  of  two  functional  categories:  component  defini¬ 
tion  or  program  control.  All  component  definition 
directives  are  considered  to  be  local  parameters.  The 
program  control  directives  may  be  either  global  or 
local.  A  global  parameter,  once  defined,  retains  its 
specified  value  until  explicitly  altered  later  in  the  data 
base.  A  local  parameter  only  retains  its  assigned 
value  until  the  occurrence  of  a  component  definition 
directive,  which  has  the  effect  of  resetting  all  local 
RCS  program  control  directives  and  then  canceling 
itself.  The  local  parameters  must  therefore  be  reas¬ 
signed  before  the  occurrence  of  the  next  component 
directive,  if  required.  The  local  RCS  program  control 
directives  and  their  defaults  or  reset  values  are: 
ORGN  (zero),  RAMDB  (zero  dB),  TRANS  (identity 
matrix),  and  VALID  (all  angles). 

As  previously  indicated,  local  parameters  must  be 
ordered  so  that  all  control  directives  related  to  a 
component  directive  appear  before  the  component  in 
the  data  base.  Global  RCS  directives  that  define 
radar  parameters  must  also  appear  before  the 
component  directives.  The  radar  directives  are: 
DELTA,  GLINT,  LOBE,  PH1R,  POLAR,  PYRT, 
RPRCS,  STATS,  THETA,  and  WLEN.  The  display 
directives  GRAPH  and  PRINT  (either  or  both  may 
be  specified)  must  appear  in  the  data  base  after  the 
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required  component  directives.  A  single  title  line 
must  always  be  placed  immediately  after  the  display 
directives.  The  two  directives  NEWTH  and 
NEWWL,  which  control  data  base  looping,  must 
always  be  placed  at  the  end  of  the  data  base,  if  used. 
If  both  directives  are  used,  NEWTH  must  precede 
NEWWL.  The  remaining  directives  may  appear 
anywhere  in  the  data  base  within  the  constraints 
given. 

The  RCS  program  may  be  instructed  to  reexecute 
the  entire  data  base  by  using  the  directives  NEWTH 
for  new  THETA  and/or  NEWWL  for  new 
wavelength  (WLEN).  The  maximum  number  of 
THETA  increments  allowed  on  one  pass  of  the  data 
base  is  181.  That  limit  may  be  extended  by  using  the 
NEWTH  directive  in  conjunction  with  the  parameter 
THETA(4).  The  actual  number  of  increments  used 
(NOI)  is  computed  as  1  +  (THETA(2) 

THETA(1))/THETA(3).  For  the  first  pass  through 
the  data  base,  the  observation  angle  ranges  from 
THETA(l)  to  THETA<2)  in  NOI  steps.  If  the 
directive  NEWTH  is  detected,  the  RCS  program 
stops  current  processing  and  starts  reading  the  data 
base  from  the  beginning.  If,  when  the  THETA 
directive  is  read,  THETA(2)  is  larger  than  or  equal  to 
THETA(4),  the  RCS  program  immediately  jumps  to 
the  first  line  following  the  NEWTH  directive.  The 
only  directives  that  may  follow  NEWTH  and 
NEWWL  are  comments.  If,  on  the  other  hand, 
THETA(4)  is  larger  than  THETA(2),  the  starling 
value,  THETA(l),  is  replaced  with  THETA(2),  and 
the  ending  value,  THETA(2),  is  replaced  with  the 
minimum  of  THETA(2)  +  NOI  and  THETA(4). 
This  processing  loop  continues  as  described  while 
THETA(4)  remains  larger  than  THETA(2).  All 
grapical  output  resulting  from  this  loop  is  placed  on 
the  same  graph.  Debugging  options  and  the  directives 
LOBE,  GLINT,  and  STATS  are  not  valid  if  THETA 
looping  has  been  specified. 

The  set  new  wavelength  directive,  NEWWL, 
instructs  the  RCS  program  to  reexecute  the  data  base 
in  a  manner  similar  to  NEWTH.  However,  in  (his 
case,  each  time  the  directive  WLEN  is  read,  the  next 
value  in  the  parameter  list  will  be  used  as  the  new 
wavelength.  For  example,  if  the  WLEN  directive  has 
three  associated  parameter  values  from  a  maximum 
of  nine,  the  wavelength  will  be  set  to  the  first  param¬ 
eter  value  on  pass  one,  the  second  value  on  pass  two, 
and  the  third  on  pass  three.  The  RCS  program  will 
not  loop  for  pass  four  because  a  fourth  parameter 
value  was  not  specified  on  the  WLEN  command  line. 
The  only  valid  lines  following  the  NEWWL  directive 
are  comment  lines. 


5.5  Data  Base  Directives 

This  section  presents  a  complete  list  of  valid  RCS 
directives  and  associated  parameters.  The  directives 
have  been  divided  into  two  groups:  RCS  control  and 
component  definition.  A  complete  description  of  the 
available  RCS  components  is  given  in  Refs.  1 1  and 
12.  All  angles  must  be  specified  in  degrees  and  dis¬ 
tances  in  meters. 

5.5.1  RCS  Control  Directives 


CLEAR 

Clear  RCS -clears  all  ac¬ 
cumulated  radar  cross-sec¬ 
tion  results. 

COMM  n 

Comment  -  ignores  the  next 
n  lines  in  the  data  base.  This 
instruction  may  be  used  to 
insert  a  specified  number  of 
comment  lines  in  the  data 
base. 

COME 

Comment  End -ends  com¬ 
ment  block  mode.  This  in¬ 
struction  will  cancel  the  ef¬ 
fect  of  COMS.  The  next  data 
base  line  will  be  processed 
normally. 

COMS 

Comment  Start  -  starts  com¬ 
ment  block  mode.  All  data 
lines  following  this  one  in  the 
data  base  will  be  ignored 
until  the  COME  directive  is 
found.  These  two  directives 
(COMS  and  COME)  provide 
the  user  with  an  efficient  tool 
for  ignoring  blocks  of  data. 

DELTA  «!,  .  .  . ,/?4  Delta  Angles  -  defines  the 


lobe  display  angles.  RCS 
lobe  and  glint  lobe  results 
will  be  displayed  at  angles 
ranging  from  «1  to  n2  in  in¬ 
crements  of  «3  degrees.  The 
number  of  lobe  samples,  be¬ 
tween  4  and  180,  is  defined 
by  m4.  Only  delta  angles  that 
lie  within  the  theta  range  will 
be  used. 

GLINT  n  Glint  Switch  -  activates  the 

glint  lobe  calculations  if  n  is 
set  to  1  or  2  and  LOBE  is  ac¬ 
tive.  Glint  data  will  not  be 
computed  otherwise.  If  n  is 
set  to  2,  the  standard  devia¬ 
tion  will  be  added  to  the  dis¬ 
placement  curve. 
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GRAPH  u',  .  .  .,n9  Graph  Display  -  collected  re¬ 
sults  will  be  output  in  graph¬ 
ical  form  (see  Section  5.3). 
For  average  cross-section 
output,  plus/minus  standard 
deviation  curves  will  be 
added  to  the  graph  if  the 
absolute  value  of  n  1  is  equal 
to  3  instead  of  1.  If  n  1  is 
negative,  terminal  graphs 
will  be  automatically  sent  to 
a  hard-copy  device.  The  x- 
axis  end  points  are  specified 
by  m2  and  m3.  The  number  of 
x-axis  labels  and  tick  marks 
per  label  are  specified  by  m4 
and  m5.  The  y-axis  end  points 
are  specified  by  m6  and  m7. 
The  number  of  y-axis  labels 
and  tick  marks  per  label  are 
specified  by  m8  and  m9.  The 
next  line  following  this  in¬ 
struction  is  always  used  for 
the  title.  This  instruction  will 
not  affect  any  RCS  results. 

GSIZE  nl,  .  .  .  ,n5  Graph  Size -locates  the 
graph  boundaries  on  the 
plotter  page.  The  lower  left- 
hand  corner  of  the  graph  will 
be  placed  at  «l,n2  inches 
from  the  paper  corner  and 
the  upper  right-hand  corner 
at  n3,n4  inches.  The  column 
of  data  displayed  at  the  right 
of  the  graph  may  be  further 
displaced  by  n5  inches.  The 
default  values  for  nl  through 
m5  are  2.25,  2.25,  8.75,  6.25, 
and  0  inches,  respectively.  If 
used,  GSIZE  must  appear 
before  the  GRAPH  di¬ 
rective. 

LOBE  n  Lobe  Switch  -  activates  the 

cross-section  lobe  calcula¬ 
tions  if  n  is  set  to  1.  Lobe 
data  will  not  be  computed 
otherwise. 

NEWTH  New  Theta  -  reruns  data 

base  with  new  theta  angle 
range.  This  instruction  will 
cause  the  RCS  program  to  re- 
execute  the  data  base  from 
the  beginning  using  the  next 
theta  range  computed  by  the 
THETA  directive.  The  in- 


NEWWL 


ORGN  x,y,z 


PAUSE 


PHIR  n 
POLAR  n 


PRINT  n 


PRYT  p,r,y,l 


struclion  will  be  ignored  if  no 
new  theta  range  is  required 
(see  THETA  and  Section 
5.4). 

New  Wavelength  -  reruns 
data  base  with  new  wave¬ 
length.  This  instruction  will 
cause  the  RCS  program  to  re- 
execute  the  data  base  from 
the  beginning  using  the  next 
wavelength  on  the  WLEN  di¬ 
rective.  The  instruction  will 
be  ignored  if  no  new  wave¬ 
length  is  available  (see 
WLEN  and  Section  5.4). 
Component  Origin -the 
component  following  this  in¬ 
struction  will  have  its  origin 
located  at  x,  y,  and  z  in  radar 
coordinates.  This  instruc¬ 
tion,  useful  for  lobe  and  glint 
computations  only,  will  be 
reset  to  zero  after  the  next 
component. 

Pause  Program- the  RCS 
program  will  stop  and  wait 
until  a  carriage  return  is  re¬ 
ceived  from  the  user  termi¬ 
nal.  This  command  should 
not  be  used  in  batch  mode. 
PHI  Radar -sets  radar  phi 
angle  to  n  degrees. 
Polarization  -  sets  system 
polarization  to  either  vertical 
(n  =  1)  or  horizontal  (n  = 
2). 

Print  Display  -  collected  re¬ 
sults  will  be  output  in  printed 
form  (see  Section  5.3). 
Printed  output  will  be  dis¬ 
played  on  the  user  terminal  if 
n  =  0.  It  will  be  sent  to  a 
nonspooled  line  printer  if  n 
=  1,  a  spooled  line  printer  if 
n  =  2,  and  a  file  if  n  =  3. 
The  next  line  following  this 
instruction  is  always  used  for 
the  title.  This  instruction  will 
not  affect  any  RCS  results. 
Pitch,  Roll,  Yaw,  Tilt  -  de¬ 
fines  radar  rotation  angles. 
These  angles  will  be  used  to 
construct  the  radar-lo-object 
coordinate  system  transfor¬ 
mation  matrix.  PRYT  will 
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default  to  the  identity  matrix 
if  not  defined.  If  all  angles 
arc  zero,  the  resultant  matrix 
will  have  its  x  and  v  com¬ 
ponents  switched.'" 

RAMDB  n  Radar  Absorption  Materi¬ 

al  -  modifies  the  effective 
cross  section  of  a  compo¬ 
nent.  The  RCS  of  the  next 
component  following  this  in¬ 
struction  will  have  n  dB 
added  to  it  before  it  is  used. 
Thereafter,  the  \alue  of  n 
will  be  reset  to  0  dB. 

RPRCS  nl,  nl  Relative  Phase  Switch -ac¬ 

tivates  the  relative  phase 
RCS  calculations  if  nl  is  set 
to  1.  Relative-phase  data  will 
not  be  computed  otherwise. 
Relative-phase  cross  section 
will  be  displayed  on  a  loga- 
rthmic  scale  unless  n2  is  spec¬ 
ified,  which  defines  the  maxi¬ 
mum  linear  value  to  be 
graphed  in  square  meters. 

STATS  n  Statistics  Switch  -  activates 

the  beta  higher-order  statis¬ 

tics  calculations  if  n  is  set  to 
1.  Higher-order  statistics  will 
not  be  computed  otherwise. 

THETA  nl,...,  n4  Theta  Angle  Range-sets 
radar  theta  range  and  incre¬ 
ment.  The  radar  observation 
angle,  theta,  will  vary  from 
nl  to  n2  degrees  in  increments 
of  n 3.  The  total  number  of 
discrete  values  may  not  ex¬ 
ceed  181.  If  the  total  maxi¬ 
mum  angle,  n4.  is  present, 
theta  will  vary  from  n\  to  n4 
degrees  in  groups  of  (n2  - 
nl)/n3  increments.  Since  this 
option  reruns  the  RCS  pro¬ 
gram  for  each  1HETA 
group,  the  number  of  dis¬ 
crete  observation  angles  al¬ 
lowed  becomes  unlimited. 
The  command  NEWTH 
must  be  included  at  the  end 
of  the  input  stream  (but  be¬ 
fore  the  directive  NEWWL, 
if  present)  to  force  the 
program  to  loop  to  the  next 
THETA  group  (see  Section 
5.4).  [f  the  parameter  n4  is 


specified,  the  following  com¬ 
mands  are  not  allowed: 
LOBE,  GLINT,  STATS,  and 
debugging  options. 

TRANS  nl,  .  .  .  ,n9  Transformation  Matrix  as¬ 
signs  component  orientation. 
The  nine  values  n  1  through 
n9,  specified  in  degrees,  de¬ 
fine,  by  row,  the  direction 
cosine  matrix  of  the  compo¬ 
nent  following  this  instruc¬ 
tion.  The  transformation 
matrix  will  be  checked  for 
validity  before  being  used  by 
a  component  and  then  reset 
to  the  identity  matrix  after 
being  used. 

VALID  n\ . /)  8  Valid  Range -sets  range  of 

valid  angles  for  component. 
The  values  of  //I  through  the 
maximum  r/8,  specified  in  de¬ 
grees,  will  define  observation 
angles  for  which  the  follow¬ 
ing  component  is  v  isible.  The 
parameter  values  must  be 
specified  in  pairs  and  imply 
the  following:  nl,  n2  =  first 
valid  range  of  THETA;  n3, 
n4  =■  associated  valid  range 
of  PH1R  (all  PH1R  if  not 
present);  n5,  n 6  =  optional 
second  valid  range  of 
THETA;  nl,  n 8  =  associ¬ 
ated  valid  range  of  PH1R. 
All  ranges  are  assumed  to  be 
counterclockwise  from  the 
first  element  in  a  pair  to  the 
second.  The  valid  range  will 
be  reset  to  all  angles  after  be¬ 
ing  used  for  a  component. 

WLEN  nl,  n2  .  .  .  Wavelength  -  assigns  pro¬ 
gram  wavelengths.  This  in¬ 
struction  sets  the  initial  pro¬ 
gram  wavelength  to  nl 
meters.  Each  following  pass 
through  the  data  base  in 
looping  mode  will  set  the 
program  wavelength  to  the 
next  n  on  the  command  line 
(pass  2  would  use  a  wave¬ 
length  of  n2,  etc.)  up  to  a 
maximum  of  n9.  The  instruc¬ 
tion  NEWWL  must  be  used 
at  the  end  of  the  data  base  if 
multiple  wavelength  passes 
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(i.e.,  looping)  are  desired  (see 
Section  5.4). 

5.5.2  Component  Definition  Directives 

CAVA  a  Cavity -defines  a  cavity  with 

equivalent  circular  area  a. 

CAVD  d  Cavity  -  defines  a  cavity  with 

circular  diameter  d. 

CFLAT  r  Circular  Flat  Plate -defines 

a  flat  plate  with  radius  r. 

CONE  A,n,a,  Truncated  Elliptical  Cone- 

L 1,  L2  defines  a  cone  with  lengths 
Z.1  and  L2  from  the  im¬ 
aginary  cone  tip.  Cone  tip  to 
LI  is  truncated.  The  elliptical 
cone  base’s  major  radius 
length  is  a  and  the  major-to- 
minor  ratio  is  n.  The  cone 
half-angle  is  defined  by  A. 

CYLIN  a,b,L  Elliptical  Cylinder  -  defines 

a  cylinder  with  major  radius, 

a ,  minor  radius  b,  and  length 
L. 

DIHED  a,b,c  Dihedral  -  defines  a  90°  di¬ 

hedral  with  sides  a  and  c  and 
height  b. 

ECORN  a,b,c  Elliptical  Corner -defines  an 

elliptical  corner  with  edges  a, 

b,  and  c. 

LWIRE  a  Loop  Wire -defines  a  cir¬ 

cular  wire  loop  with  radius 

a. 

OGIVE  R,a,b,A  Ogive -defines  an  ogive  with 

radius  of  curvature  R  and 


center  to  nearest  side  length 
a.  If  ogive  is  truncated,  pa¬ 
rameter  b  and  angle  A 
specify  the  amount  of  trun¬ 
cation  (b  is  zero  otherwise). 

PARAB  p  Parabola  -  defines  a  parabo¬ 

la  where  the  parabola  equa¬ 
tion  is  p  —  (,v**2  + 

y**2)/  -4z. 

RCORN  a,b,c  Rectangular  Corner- defines 

a  rectangular  corner  with 
edges  a,  b,  and  c. 

RFLAT  a,b  Rectangular  Flat  Plate -de¬ 

fines  a  flat  plate  with  sides  a 
and  b. 

SOLID  a,b,c  Solid  -  defines  an  ellipsoid 

with  radii  a,  b,  and  c. 

TCORN  a,b,c  Triangular  Corner  -  defines 

a  triangular  corner  with 
edges  a,  b,  and  c. 

TORUS  a,b  Torus  -  defines  a  torus  with 

inner  radius  a  -  b  and  outer 
radius  a  +  b. 

TWAVE  p,L,g  Traveling  Wave -defines  a 

traveling  wave  with  the  pa¬ 
rameters  p  (relative-phase  ve¬ 
locity),  L  (length),  and  g  (re¬ 
flection  coefficient). 

WEDGE  A,L,B  Wedge- defines  a  wedge 

with  angle  A  and  height  L. 
Parameter  B  is  used  to  limit 
the  cross  section  to 
S/WLEN**2. 

WIRE  L,a  Wire -defines  a  wire  with 

length  L  and  wire  radius  a. 
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6.0  PROGRAM  LISTING 


>>>  FILE  *  SYO : [ 1 1 1 ,  1 2 ]RCS . FTN 


27-SEP-83  14:32:52  <«< 


RRRR 
R  R 
RRRR 
R  R 
R  R 


RADAR 


CROSS 


SECTION 


NOTE!  THIS  PROGRAM  STORES  MATRICES  BY  ROWS  (NOT  COLUMNS) 

MAXRCS  CONTAINS  THE  CURRENT  MAXIMUM  RCS  BUFFER  SIZE  (1811 
REALMS  T1 ,T2,T3,T4lT5,T6lSQ,ZERO 

REALM  CMDSI26) .OBJSII8) , RCS (181) , RCS2 (181 ) .TRANS (9) , VAL 10(8) 
REALM  ARRI9) , PRYT 1 9 ) .MAT ( 9 ) .C0MPI9) , PRYTVI4) , PRYT2 (9 ) 

LOO  I  CALM  ATTN 

LOGICALM  F I LE (80) , BUF (80) , SLASH , AG ,AP . OBJECT , AZ , AC , AN , DOT , DRCS (4) 
REALM4  POLVH ( 2 ) , YYV (181 , 2  I , XV  1 1 81 ) . SLOBE ( 4 , 1 81 .2) 

EQUIVALENCE  (FILE. BUF I , (YYV, SLOBE) . ( XV , SLOBE ( 1 . 1 , 2 ) ) 

REALM4  I  DENT (9) . GRARR ( 9 ) , S IDE (3 . 4) . DEBOUT  1 6) , ORGN 1 3 ) .M80DB 
REALM4  DELTA (4) . ALPHA ( 4 1 , SCALE  1 6 ) , THETA ( 4 ) .GSIZEI5) 

INTEQERM2  LSIDE (4) .ROUTE 

INTEGERM2  DEBUG .COMB  I N ,GL I  NT 

LOG ICALM2  GOULD, LP , OPENS . LPF .OLDGRF .MQRAPH.GDAUTO 
COMMON /SW I TCH/  I OW, THETA , GOULD , LP , OPENE , LPF , OLDGRF .MGRAPH . 

1  GDAUTO, DEBUG, COMB  IN, GLINT, I  BUG, LBUG, I  STATS, I ZETA , STATS . NRCS , IOR 

COMMON /RCSCOM/CUTOFF , WLEN , P I , P I  2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 
9  C2PHP.S2PHP, RAD, POLAR, LOBE, 

9  ARR 

COMMON /RCS LOB/ I OL , ORGN , COMP , SLOBE , BETA , 

9  ANG , CANG , SANG , PHR , CPHR , SPHR , KLOBE 
COMMON  /RPRCS/  I  RPR , I RCS , RAM , RCS , RCS2 
COMMON  /GRCOM/  GRUB 

DATA  POLVH/’ VER’ HOR’ / .CLEAR/ ’ CLEA ’/ ,COME/ ’ COME ’/, SLASH/ ’/’ / 

DATA  CMOS/ ’ COMM' ,'COMS’ ,'WLEN’ ,’POLA’ ,’TRAN’ ,’VALI’ ,’PHIR’ ,’THET' , 
3  ’GRAP’ , 'PR IN’ , 'PRYT' . 'PAUS' , 'CLEA' ,  DEBU’ , ’NEWW' , ’ RAMD’ , 'LOBE' , 

9  'BETA'  ,  'ORGN'  ,  'DELT'  ,  ALPH'  ,  RPRC’  ,  'STAT'  ,  'GUN'  ,  'NEWT'  ,  GSIZ’/ 

DATA  OBJS/'OGIV ,’CYLI ' , 'CONE' , 'PARA' , ’ SOU ’ , 'TORU' , ’ WEDG’ . 

3  'WIRE' , ’ LWIR' , 'CAVA' , 'CAVD' , ’ CFLA ' , ’ RFLA ’ , 'TWAV , 'DIKE' , 'TCOR' , 

I  » RCOR 1  ' ECOR 1 / 

DATA  i DENT/1 . ,0 . , 0 . , 0 . , 1 . , 0 . , 0 , , 0 . , 1 . / , AG/ 'G' / , AP/ ’ P ’ / , AC/ ' C’ / 

DATA  DEBOUT/O. , ’  COM’ , ’ PONE ’ , ’ NT  O' , ' UTPU' , ’ T’ / , ALL/’ ALL  '/ 

DATA  GRARR/ -1 ,0,180,7,3,-60,60,7,1/, AZ/ ' Z’ / , AN/’ N' / , DOT/ ’ .’/ 

DATA  S I DE/’ STAR ’ ,’BOAR’ ,’D’ ,’PORT’ ,0,0, 'TOP' ,0,0,’BOTT’ ,’OM' ,0/ 
DATA  LSIDE/7 ,4,3 ,6/ ,DRCS/' , ’ , 'R' , 'C ,’ S’ /, ZERO/O . 00/ 

DATA  GSIZE/2.25,2.25,8.75,6.25,0./ 

DEBUG  OUTPUT 

I BUQal s>RCS  OUTPUT 
I BUGs2s> TRANS, PRYT, COMP 
I  BUQMa>  ANGLES  (RADAR!  PR  I  ME ) 
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600> 

C 

6 1  0  > 

C 

SET  RCS  CONSTANTS 

620> 

C 

630> 

NCMO-26 

640> 

NOBJ  =1 8 

oso> 

MAXRCS=1 81 

660> 

Pl=3. 1415926 

670> 

PI2SPI/2. 

680> 

RADsPI/180. 

690> 

CUTOFF=1 .E-10 

700> 

LINE=0 

710> 

G0ULD=. FALSE, 

720> 

C 

730> 

C 

740> 

C 

ASS  1 QF  LUNS  AND  ATTACH  INPUT  FILE 

750> 

C 

760> 

1  OR  =  1 

770> 

lOWs  2 

780> 

IOPs3 

790> 

100=4 

800> 

1  OF  =  5 

8 1  0  > 

1 0L  =  6 

820> 

10E  =  7 

830> 

C 

840» 

CALL  ASNLUN ( i OR , ' T 1  * ,0) 

850> 

CALL  ASNLUN I I0W, ' T 1 ' .0) 

860> 

CALL  ASNLUN! t OF , ’ SY  * ,0) 

870> 

CALL  ASNLUN! IOP, ’GO’ ,0) 

880> 

CALL  ASNLUN! IOL, ’SY* .0) 

890> 

CALL  ASNLUN! !OE, 'SY' ,0) 

900> 

OPINE=. FALSE. 

91 0> 

C 

920> 

WRITE! IOW.1) 

930> 

1 

FORMAT!’  ENTER  RCS  DATA  FILE  NAME’) 

940» 

READ  I  1  OR, 2, END  =  9999)  FILE 

950> 

2 

FORMAT I80A1  ) 

960> 

CALL  TR IM80 (NF .FILE) 

970> 

FILE! NF ♦ 1 ) =0 

980> 

C 

990> 

C 

INPUT  OPTIONS 

1000> 

c 

1010» 

c 

/z 

s  DRAW  GRAPHS  ON  ZETA  PLOTTER 

1 020> 

c 

/p 

*  PRINT  RESULTS  FOR  EACH  COMPONENT 

1 030> 

c 

/G 

s  GRAPH  RESULTS  FOR  EACH  COMPONENT 

1 040> 

c 

/QC 

s  DRAW  ALL  /G  OUTPUT  ON  ONE  GRAPH 

1 050> 

c 

1 080> 

DEBUQsO 

1070* 

COMB  1 N=0 

1 080» 

IZETAsO 

1 090> 

1  F  =  1 

1 1 00> 

IDOTsO 

1 1 1 0> 

DO  6  JsIF.NF 

1 120> 

JJ=J 

1 1 30> 

IF(FILEUJ)  .EQ.DOT)  1  DOT  =  J J 

1 1 40> 

IFIFILEIJJ) .EO. SLASH)  GO  TO  7 

1 1 50> 

6 

CONTINUE 

1 1 80> 

GO  TO  8 

1 170> 

7 

NFsJJ-1 

1 1 80> 

IFIFILEIJJ*)). EQ. AZ)  IZETA=200 

1 190> 

1 F (F  1 LE I J J*1 ) . EQ . AG )  DEBUG  =  1 
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1  200> 

1 21 0> 

1  2  2  0  > 
1230*  9 
1  2  4  0  >  C 
1  2  5  0  > 

1  2  6  0  > 

1 270> 

1 280>  1 
1 290>  9 
1 300> 

1  3 1  0  >  C 
1  320> 

1  3  3  0  > 

1 340> 

1  350> 

1  360> 

1  370>  C 
1  380>  C 
1  3  9  0  >  C 
1  400> 

1 41 0> 

1  4  2  0  > 

1  4  3  0  > 

1 440>  C 
1 450>  C 
1 460>  C 
1 470>  C 
1 480>  C 
1 490>  C 
1 500>  3 
1  5 1  0  > 

1  520> 

1 530>  C 
1  540> 

1  550> 

1 560>  4 
1 570>  C 
1 580>  C 
1 590>  C 
1 600> 
1610> 

1  020> 

1  6  3  0  > 

1 840>  5 
1 650>  C 
1 660>  C 
1 670>  C 
1 680> 

1 690>  C 
1 700>  C 
1  7 1  0  >  C 
1 720>  C 
1  7  3  0  >  C 
1 740>  C 
1 7  5  0  >  C 
1 760>  C 
1 770»  C 
1 780>  C 
1 790>  1 


I F I F I LE  <  J J  ♦  1 1 . EQ . AP )  DEBUG: 2 

IF (DEBUG. EQ.1 . AND . F I LE I J J *2 ) . EQ . AC )  COMBIN=-1 

I F  s  J  J  ♦  2 

I F 1 1  DOT. QT. 01  00  TO  9 
DEFAULT  FILE  EXTENT  I  ON  TO  .RCS 
DO  13  J  =  1 ,4 
NF=NF*1 

F  I  LE  I NF ) =DRCS ( J ) 

3  CONTINUE 

FILE! NF  *11=0 

OPEN ( UN  I T= I  OF, NAME =F I LE, TYPE  = 'OLD'  , READONLY ) 

T IME  =  SECNDS (0.0) 

CALL  OR  I N I T 1 0 , IZETA, I0G, I ER) 

CALL  QRHVI ' H' I 

CALL  GRPICI2. 25, 2.25,8.75, 6. 25) 

IF( IZETA. EQ.O)  CALL  GR0EV ( ’ TEK ’ , ’ HARD ’ , I ER I 

INITIALIZE  RCS  VECTOR  A  OTHER  STUFF 

CMD30 . 

1WL  =  1 

OLDGRFs. FALSE. 

GO  TO  70 


* NMNttftttttMNNNM * ***** ** N * **##### * 

READ  NEXT  LINE  IN  RCS  DATA  FILE 
NttNftMNNttKNMNttNMMNNtttflHHHHHHtMNMft 

CONTINUE 

IFIOBJECT.AND.DEBUG.EQ.il  GO  TO  45 
I F (OBJECT  AND. DEBUG. EQ. 2)  GO  TO  50 

'.  I  NE  =  L  I NE  +  1 

READ! IOF,4,ENO=9990,ERR=9993I  CMD.ARR 
FORMAT (A4,1X,9F10.0) 

DETERMINE  IF  CMD  IS  A  DIRECTIVE 

DO  5  J=1 , NCMO 
I F  I  CMD . EQ . CMOS  I J 1 1 

9  GO  TO  (10,12.15.20.25,30,35,40,45,50,55,50,70,80,85,65, 
8  500,510,520,540,550,560,570,580,590,6001,  J 
CONTINUE 

GO  TO  OBJECT  SEARCH  IF  NO  DIRECTIVE  MATCH  FOUND 
GO  TO  90 


EXECUTE  APPROPRIATE  DIRECTIVE 


"COMM”  -  SKIP  COMMENT  LINES 
0  NCOMM=ARR ( 1  I 
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1800» 

1 8 1 0  > 

1 820> 

1 930> 

11 

1 840> 

1 850> 

C 

1 860> 

12 

1 870> 

1 880> 

1 890> 

1  900> 

C 

1 9 1 0  > 

C 

1 920> 

C 

1 9  3  0  > 

IS 

1 940> 

1  9  S  0  > 

1  960> 

1 970> 

C 

1  980> 

C 

1  990> 

C 

2000» 

20 

2010> 

2020> 

C 

2030> 

C 

2040> 

C 

20S0> 

25 

2060> 

2070> 

26 

2080> 
2090> 
2 1  00  > 

C 

2110> 
21  20> 
2 1 3  0  > 

24 

2 1  40> 
21  50> 

27 

2 1  6  0  > 

28 

2 1  7  0  > 
21  80> 

29 

2 1  9  0  > 
2200> 
2  2 1  0  > 
2220> 

C 

2230> 

C 

2240> 

C 

22S0> 

30 

2260> 

2270> 

31 

2280> 

2290> 

C 

2300> 

C 

2  3 1 0  > 

C 

2320> 

33 

2330> 

2340> 

36 

23S0> 

2360> 

C 

237D> 

C 

2380» 

C 

2390> 

40 

IF (NCOMM. LE .0)  GO  TO  3 
00  11  ,1  =  1, NCOMM 
L INE=L INE  +  1 

READ  I  I  OF, 4, END  =  99901  CMD 
GO  TO  3 

"COMS"  -  SEARCH  FOR  "COME* 

RE AD (I OF, 4, END  =  9993)  CMD 
LINEUINE  +  1 

IF (CMD. EQ. COME)  GO  TO  3 
GO  TO  12 

“WLEN*  -  SET  WAVELENGTH 
WIEN=ARR I IWL I 

IF IWLEN.LE.O.AND. IWL.EQ.1 )  GO  TO  99P1 
I F I WLEN . LE . 0 )  GO  TO  9990 
GO  TO  3 

"POLAR*  -  SFT  POLARIZATION  VECTOR 

POLAR=AMAX1  11 .  ,  AM  INI  (2.  .ARRI1M) 

GO  TO  3 

"TRANS"  -  COMPUTE  COORDINATE  TRANSFORMATION  MATRIX  1 

CONTINUE 
DO  26  J  =  1  ,9 

TRANS  I J ) =COS (ARRIJ)MRAD) 

CALL  TRANSP ( TRANS , MAT) 

CALL  MATMUL I  TRANS , MAT .MAT) 

MARE  SURE  TRANS  IS  VALID 
DO  24  J=1 ,9 

IF ( ABS (MAT (J)-IDENT(J) ) .QT.1 .E-3)  GO  TO  27 
CONTINUE 
GO  TO  3 
WR  I TE ( 2 , 26 1 

FORMAT  I ’  RCS>  --  BAD  TRANS  DIRECTION  COSINE  MATRIX’/’  RCS>'I 
WR I T E (2 ,29)  ARR 
FORMAT  I’  RCS>  '.3F10.3I 
WR I TE  (2 ,29 ) 

WR I TE ( 2 , 29 1  MAT 
GO  TO  9993 

"VALID"  -  SET  VALID  ANGLE  REGIONS 

CONTINUE 
DO  31  J=1 ,3 
VAL I D ( J ) =  ARR  t  J ) 

GO  TO  3 

"PHIR"  -  SET  RADAR  PHI  ANGLE 
PHR=ARR ( 1 ) 

I F I ABS  IP HR) . QT . 1 80 . )  WRITE  I IOW.36)  PHR 
FORMAT!’  RCS>  --  PH  I R= ’ , F6 , 1 , ’  I S  OUT  OF  BOUNDS’) 

GO  TO  3 

"THETA"  -  SET  THETA  RANGE  (MIN, MAX, INC)  A  TOTAL  MAX 
IF (OLDGRF)  GO  TO  43 
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2400 

2410 

2420 

2430 

2440 

24S0 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 


THETA (1 )sARR(1 ) 

IF (ARRI3) . LE . 0 . )  GO  TO  42 
THETA  12 )=AMAX1 I  THETA (1) , ARR  12) ) 

THETA  13 ) =ARR (3 ) 

THETA  14) =  ARR (41 
MGRAPHs. FALSE. 

I F ( THETA (4) . GT . THETA  1 2 ) )  MGRAPHs . TRUE . 

C 

NRCS=(ARRI2)-ARR11 ) )/ARR<3)*1 .001 
I F (NRCS . LE .MAXRCS)  GO  TO  3 
WRITE! IOW, 41 ) 

41  FORMAT)’  RCS>  --  TOO  MANY  OUTPUT  POINTS’) 

GO  TO  9993 
C 

42  THETA (2 1 sARR !  1 1 

THETA ( 3 ) si  . 

NRCSsi 
GO  TO  3 
C 

43  T  =  THETA  1 1 ) 

THETA  1 1 ) =THETA (2 ) 

THETA (2) sAMI N1 ( 2 . HTHETA ( 1 ) -T . THETA (4 ) ) 

GO  TO  3 
C 

C  ’’GRAPH"  -  GRAPH  RESULTANT  RCS  VECTOR  s  CMD 1 9 ) 

C 

45  CONTINUE 

C  I F I  I  STATS . EQ. 1 )  GO  TO  4542 

IFINRCS.LE.1 .OR. DEBUG. EG. 2)  GO  TO  3 
I F 1 1 ZETA . EQ. 0 .AND .COMB IN . LE . 0 .AND . .NOT.OLDGRF) 

1  CALL  GRMODEC  ERASE’  ,3) 

ATTNs. FALSE. 

IF ICMD.EQ.CM0SI9) )  GO  TO  4546 
C 

C  GRAPH  DEBUG  MODE 

M60DB=1 ./10**6 
DO  4541  Jsl , NRCS 
IFIRCSIJ) .GT.M60DBI  GO  TO  4540 
4541  CONTINUE 

GO  TO  4558 
C 

4540  CALL  TRAPCCI IOR.ATTN) 

DO  4545  J=1 .9 

4545  ARR I  J) =GRARRI J ) 

4546  CALL  GRSCL I ARR 12 ) , ARR 1 6 ) , ARR (3 ) , ARR  1 7 ) ) 

IF  l  IRPR.GT.O. AND. RPRCSG.GT.O. J CALL  GRSCL (ARR 12) .0. , ARR (3 ) , RPRCSG) 
SCALE (1 )sARR(2) 

SCALE  1 2 ) =ARR 1 6 ) 

SCALE  1 3 ) =ARR 1 3 ) 

SCALE (4) =ARR ( 7 1 
SCALE ( 5 ) sARR 1 4 ) 

SCALE  16) sARR 1 5 ) 

I F (COMB IN.EQ.1 . OR . OLDGRF )  GO  TO  4547 
C 

YSQSIZEC2) 

IF IKLOBE . LE .0)  GO  TO  44 
YsY-,95 

CALL  GRP LOT (GS I ZE (3 1 , GS I ZE (2  I ,01 
CALL  GRP LOT  I GS I ZE ( 1 ) , GS I ZE 1 2) , 1 ) 

C 
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3000 

3010 

3020 

3030 

3040 

3050 

3000 

3070 

3080 

3000 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3100 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3200 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3300 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3475 

3480 

3400 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 


44  CONTINUE 

CALL  GRAX I S ( 0 , GS I ZE ( 1 )  ,Y,GSIZE(3)-GS!ZE11)  ,0.  ,ARR<2)  ,ARR(3)  , 

1  IFIX(ARRU) ),  IFIX(ARRISI)  ,  ’OBSERVATION  ANGLE', 00,1) 
DY=GSIZE(4)-GSIZE<2) 

I F ( IRPR . EQ . 0 .OR . RPRCSG . EQ. 0 . )  GO  TO  4560 
C 

CALL  GRAX I S (3 ,G$ I ZE ( 1 ) .GSIZE12) ,DY,80. ,0. .RPRCSG, LRPRCS, 2, 

1  'RC3  ISM) ' ,8,1) 

GO  TO  4547 
C 

4560  CALL  GRAX I S (3 ,GS I ZE (1) , GS I ZE (2) , DY , 00 . , ARR (6 ) , ARR 1 7 ) . IF  IX (ARRIS)) , 
*  IF  IX (ARR 10 ) ) , ' RCS  IDBSM) ' ,10,1) 

C 

4547  ANQsTHETA 1 1 ) 

DO  46  J  =  1 , NRCS 
XV  U)  SANG 
ANG=ANG+THETA(3) 

YYYsRCSIJ) 

IF! IRPR.EQ.O)  GO  TO  4561 
YYYsRCS (J ) **2*RCS2 ( J ) **2 
I F I RPRCSG . EQ . 0 . )  GO  TO  4561 
YYVU,1)sYYY 
GO  TO  46 

4561  YYVIJ ,1 ) = 1 0 . #ALOG1 0 ( AMAX1 (CUTOFF, YYY) ) 

46  CONTINUE 

CALL  GRWI ND ( -NRCS , XV, YYV, 16,0, .1,0. ,1) 

I F (COMB I N . EQ . 1 )  GO  TO  4558 

IF(AB$(ARR(1J) . LE.1. OR. DEBUG. EQ.1. OR. IRPR.EQ.il  GO  TO  49 
S I GN  =  - 1  . 

DO  47  K=1 ,2 
DO  48  4*1 , NRCS 

S^SQRT ( AMAX1  ( 0 .  , RCS ( J ) M2-RCS2 ( J ))) 

48  YYV ( J ,K ) =1 0 . *ALOQ1 0 ( AMAX1 (CUTOFF , RCS (J I *S IGNNS) ) 

CALL  GRWI ND( -NRCS, XV, YYV (1 ,K) ,IB,0, .1 ,0. ,1) 

47  SIQNal. 

49  I F ( IZETA.EQ.O)  CALL  GAMODE ( ’ ALPHA ’ , 2 ) 

I F (CMD . EQ . CMOS (9 ) I  GO  TO  4548 
DEBOUT  1 1 ) =CMD 

I F ( COMB  I N . EQ. -1 )  DEBOUT ( 1 J =ALL 
CALL  GRTXT (DEBOUT , 21 ) 

GO  TO  4549 

4548  READ ( IOF,2,END=9993)  BUF 

L I NE=L I NE ♦ 1 
I F (OLDGRF )  GO  TO  4559 
CALL  TR IM80 (NB , BUF ) 

CALL  GRTXT (BUF ,NB) 

QDAUTOs. FALSE. 

IF (ARR ( 1 1 . LT.O. )  GDAUTO- . TRUE . 

4649  XMI D= (GS I ZE (3 ) *GS t  Z  E ( 1 ) ) / 2 . 

CALL  GRSXYIXMID,GSIZEI4)*.75,XX,YY) 

CALL  GRPRNT (0,XX,*1,YY,.15,0.) 

CALL  GRTXT ( 'MONOSTAT 1C  FULL-SCALE  ’,0) 

I F ( IRPR.EQ.O)  CALL  GRTXT (' RCS ’, 0) 

I F ( IRPR. EQ.1 )  CALL  GRTXT (’ RPRCS’ , 0) 

CALL  GRSXY(XMID,GSIZEI4)*.5,XX,YY) 

CALL  GRPRNT (0,XX,-1,YY,.15,0.) 

C 

C  SET  TERMINAL  TO  ALPHA  MODE  AND  SCALE  TO  INCHES 
IF( IZETA.EQ.O)  CALL  GRMODE ( ' ALPHA ' , 3 ) 

CALL  GRSCL (GS I ZE ( 1 ) ,GS I ZE 1 2 ) , GS I ZE ( 3 ) ,QS I ZE (4) ) 
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3590> 

X=GS 1 ZE I  3 ) +1 . 05  +  GS 1 ZE ( 5 ) 

3600> 

C 

361  0> 

YsQSIZE (41 ♦  .75 

3620> 

FREQs. 3/WLEH 

3630> 

CALL  GRTX7( ’FREQUENCY  ’ ,0) 

3640> 

CALL  GRPRNTI1 ,X,0,Y, .12,0.) 

3650> 

CALL  GRNUMI ’ F ’ ,-1 ,3, FREQ) 

3660> 

CALL  GRTXT I ’  GHZ’ ,0) 

3670> 

CALL  GRPRNT ( -1 , X , 0 , Y , . 1 2 , 0. ) 

3680> 

YSY-.2 

3690> 

CALL  GRTXTI ’POLARIZATION  ’,0) 

3700> 

CALL  GRPRNT 1 1 ,X,0,Y, .12,0.) 

371 0> 

CALL  GRTXT (POLVH ( 1 F 1 X (POLAR) ) , 3 ) 

3720> 

CALL  GRPRNT 1 -1 , X , 0 , Y , . 1 2 , 0 . ) 

3730> 

IF (MGRAPH)  GO  TO  4542 

3740> 

YsY-,2 

3750> 

PERSHsi 00.-1 00. MFLOAT ( NVAL 1 D) /AMAXO (1 .NTOTAL) 

3760> 

CALL  GRTXT ( ’%  SHADOWED  ’,0) 

3770> 

CALL  GRPRNT ( 1 , X , 0 , Y , . 1 2 , 0 . ) 

3780> 

CALL  GRNUMI ’ p ’ ,-1,1 .PERSH) 

3790> 

CALL  QRPP.Nn-1  ,  X ,  0 ,  Y ,  .12,0.) 

3800> 

4542 

Y=Y- .2 

3  8 1  0  > 

CALL  GRTXTI 'PITCH  ' ,0) 

3820> 

CALL  GRPRNT ( 1 ,X,0,Y, .12,0.) 

3830> 

CALL  GRNUMI ’ F ’ ,-1,3, PRYTVl 1 ) ) 

3840> 

CALL  GRPRNT I -1 ,X,0,Y, .12,0. ) 

3850> 

Y  =  Y-  .2 

3860> 

CALL  GRTXTI’ ROLL  ' ,0) 

3870> 

CALL  GRPRNT 1 1 ,X,0,Y, .12,0.) 

3880> 

CALL  GRNUMI ’ F ’ , *1 ,3 , PRYTV 12)) 

3890> 

CALL  GRPRNT 1-1 ,X,0,Y, .12,0.) 

3900> 

Y=Y- .2 

391  0> 

CALL  GRTXTI’ YAW  *  ,0) 

3920> 

CALL  GRPRNT 1 1 ,X , 0, Y , .12,0.) 

3930> 

CALL  GRNUMI ’ F ’  , -1 , 3 , PRYTV 1 3 ! ) 

3940> 

CALL  GRPRNT I -1 , X , 0 , Y , .12,0.) 

3950> 

YsY-,2 

3980> 

CALL  GRTXTI  'TILT  ’  ,0) 

3970> 

CALL  GRPRNT 1 1 , X , 0 , Y , . 12,0.  ) 

3980> 

CALL  GRNUMI ' F ’ , -1 , 3 , PRYTV ( 4) ) 

3990> 

CALL  GRPRNT 1-1 ,X,0,Y, ,12,0. I 

4000> 

Y=Y-.2 

401 0> 

ISIDEsO 

4020> 

IFIPHR.EQ.90. )  1 S 1 D E  =  1 

4030> 

IFIPHR.EQ.-90. )  1 S 1 DE=2 

4040> 

IFIPHR.EQ.O. )  1 S 1 DE-3 

4050> 

IFIPHR.EQ.180. )  ISIDE&4 

4060> 

IFI ISIDE.EQ.O)  GO  TO  4556 

4070> 

CALL  GRTXT IS  IDE  11 , 1 S 1 DE I , LS 1 DE 1 1 S 1 DE ) ) 

4G80> 

CALL  GRPRNT (0 , X , 0 , Y , .12,0.) 

409  0> 

GO  TO  4557 

41 00> 

4556 

CALL  GRTXTI 'PHIR  ’ ,0) 

41 1 0> 

CALL  GRPRNT 1 1 , X , 0 , Y , .12,0.) 

41  20> 

CALL  GRNUMI ’ F ’ ,-1 ,3,PHR) 

41 30> 

CALL  GRPRNT I -1 ,X,0,Y, .12,0.) 

41 40> 

4557 

CONTINUE 

4 1  5  0  > 

1 F IMGRAPH )  GO  TO  4559 

41 60> 

IF (DEBUG. NE .0)  GO  TO  4550 

4 1  70  > 

C 

41  80> 

C 
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4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4280 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 


C 

IFCKLOBE.LE.O)  GO  TO  4551 
CALL  GRPL0T (11. ,0. ,0) 

IGRPR:  1 

WR I TE ( I OL ' 1 )  IGRPR , KLOBE .GLINT, NRCS , NLOBE , THETA .ALPHA .DELTA , 

1  BETA , SCALE , IZETA , I PCODE ,WLEN , LBUQ , NB , BUF , GDAUTO , 

2  FREQ, POL VH( I  FIX (POLAR) ) , PERSH, PITCH, ROLL. YAW, T I LT.PHR.G8I2E 
CLOSE (UN  I Ts IOL ) 

CALL  SPAWNS 

OPEN  I  UN  IT: I OL, NAMES’ LOBE. TMP’ ,TYPE=’0LD’ , RECORDS  I ZEsI  81  , 

2  ACCESS: ' DIRECT’ .SHARED) 

C 

4551  IF  I ISTATS.EQ.O)  GO  TO  4550 
C 

C  GRAPH  HIGHER  ORDER  BETA  STATISTICS  IF  REQUIRED 

C 

CALL  HOSBDI1 .DUMMY, IOE, DUMMY, SLOBE, BUF, NB, 

1  FREQ, POL VH( IF  IX (POLAR) ) .PERSH, PITCH, ROLL, YAW, TILT, PHR) 

C 

C  CLEANUP  GRAPHS 

C 

4550  IF1QDAUTO. AND. IZETA. EQ.O. AND. DEBUG. EQ.O)  CALL  GRTEKG 

IFIDEBUG.EQ.1 .AND. COMBIN.EQ.O. AND. .NOT. ATTN)  CALL  GRTEKG 
4559  IFI IZETA. EQ.O)  CALL  GRPLOTIO. ,8 .4,0) 

IF  I IZETA. EQ.O)  CALL  QRMODEI ’ALPHA’ ,3) 

IFI IZETA. GT.O. AND. .NOT. MGRAPH)  CALL  GRPLOT (1 1 . ,8 . 5 , -1 ) 

IFI IZETA. QT.O. AND. MQRAPH)  CALL  QRUPPN 
IF ICMD.EQ.CMDSI9) )  GO  TO  3 
4558  COMBI N= I ABS ICOMB IN) 

I F (ATTN)  CALL  DETACH! IOR, IER) 

CMD=CLEAR 

CALL  GRPLOTIO. ,8.4,0) 

CALL  QRMODEI ’ALPHA’ ,3) 

GO  TO  70 

4544  READ  I IOF , 2 )  BUF 

L I NE=L I NE+1 
GO  TO  3 
C 

C  "PR I  NT"  -  PRINT  RESULTANT  RCS  VECTOR 

C 

SO  CONTINUE 

IF (DEBUG .GE . 1 . AND .CMD . EQ.CMDS I 1 0) )  GO  TO  4544 
10= IOW 

C  FORCE  OUTPUT  TO  LINE  PRINTER  FOR  DEBUG 
I PCODE=ARR 1 1 ) 

IFIARRI1 ) .GE.1 . 0.OR.CMD. NE. CMOS  1 10))  10* IOP 
ANQsTHETAM  ) 

ANGINC=THETA(3) 

LP=. FALSE. 

LPFs. FALSE. 

IFIARRI1) .GE. 2.0. AND. DEBUG. EQ.O)  LPs . TRUE . 

IFIARRI1) .GE. 3.0. AND. DEBUG. EQ.O)  LPF* . TRUE . 

IFILP)  CALL  ASNLUN ( I  OP , ’ SY  * , 0 ) 

IFILPI  OPEN (UN IT= IOP, NAME: ’LP.LST’ ) 

IFI IO.EQ. I  OP .AND. .NOT. IGOULD.OR.LP) )  CALL  ATTACH! 10, IER) 

IFI 10. EQ. I  OP. AND. .NOT.LP)  GOULD:. TRUE. 

IF (CMD . EQ.CMDS 1 101 1  GO  TO  5050 
C 

C  DEBUG  OUTPUT  PRINT  MODE 
M60DBS1 ./10*«6 
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4790> 

4800> 

481 0> 
4820> 

5054 

4830> 

4840> 

5055 

4850> 

C 

4860> 

4870> 

5053 

4880> 

4890> 

56 

4900> 
491  0> 
4820> 
4930> 

5050 

4940> 

4950> 

4960> 

54 

4970> 

4980> 

4990> 

5057 

5000> 

C 

501 0> 

5051 

5020> 

5030> 

5040> 

51 

5050> 

5060> 

57 

5070> 

C 

5080> 
5090> 
5)  00> 
51 1 0» 

5 1 2  0  > 
51 30> 

C 

51 40> 
51 50> 

5059 

5 1 60> 

C 

51 70> 
51 80> 
51 90> 
5200> 
521 0> 
5220> 

5058 

5  2  3  0  > 

53 

5240> 

5250> 

58 

5260> 

C 

5270> 

C 

5280> 

C 

5290> 

5300> 

C 

531 0> 

C 

5320> 

S330> 

52 

5340> 

C 

5350> 

C 

5360> 

5370> 

5380» 

C 

DO  S0S4  J=1 , NRCS 

IFIRCSIJ)  .GT.M60DB)  GO  TO  5053 

CONTINUE 

WRITE!  10,505$)  CMD.LINE 

FORMAT! '  0*  ,A4,’  HAS  NO  RCS  CONTRIBUTION  (DATABASE  LIHE’,I3,T) 
GO  TO  S0S6 

CONTINUE 

WRITE!  10,56)  CMD.UNE 

FORMAT  (/' O' ,A4, '  COMPONENT  OUTPUT  FROM  DATABASE  LINE  MS) 

GO  TO  50S1 

READ! I  OF ,2 ,END=9993)  BUF 
L I NE  =  L I NE  +  1 
CALL  TRIM80(NB,BUF) 

WRITE! 10,54)  (BUF ( I ) . I *1 ,NB) 

FORMAT (/’O’ , 1 20A1 ) 

FREQs. 3/WLEN 

WRITE! 10,5057)  FREQ 

FORMAT!'  FREQUENCY (GHZ) *• ,F5.1) 

I BGS I  BUG . AND . 1 
I F ! CMD . NE . CMOS (10))  IBQ=1 
I F ( I  STATS .EQ. 1 )  IBQsO 
IF ( IRPR.EQ.O)  WRITE! 10,51 ) 

FORMAT! '0  RCS  OUTPUT ' /2X ANGLE ', 5X ,’ RCS' . 5X , 

S  ’  RCS-S’ , 5X , ' RCS*S ’ /IX) 

IF! IRPR.EQ.1 I  WRITE! 10,57) 

FORMAT!  0  RPRCS  OUTPUT ’ /2X ,* ANGLE ', 2X ,’ RCS IDBSM) ', 2X , 

1  ’ RCS ( SM) ' ) 


DO  52  Jsl.NRCS 
IF  I IRPR.EQ.O)  GO  TO  505B 
YsRCSU)»«2*RCS2IJ)fttt2 
YYYsl0.*A LOG 10 ( AMAX1 (CUTOFF, Y )) 

WRITE! 10, 5059 )  ANG.YYY.Y 
FORMAT IF7 .1 , 2F1 0 . 4 1 
GO  TO  58 

RCSC=AMAX1 (CUTOFF, RCS(JJ) 

S=SQRT ( AMAX1 ( 0 . , RCS (J ) HH2-RCS2 ( J ) ) ) 
DBSM=10.*ALOG10(AMAX1 (CUTOFF , RCSC) ) 
DBSMMSslO.HALOGIOIAMAXI (CUTOFF, RCSC-S) ) 
DBSMPSsIO.hALOGIOIAMAXI (CUTOFF , RCSC+S) ) 

WRITE! 10,53)  ANG , DBSM , DBSMMS , DBSMPS 
FORMAT (F7.1 , 3F1 0 . 4 , 3X , 2F1 0 . 4 ) 

IF! ISTATS.EQ.O)  GO  TO  52 

IFtSTATS. QE.ANGINC. AND. AMODiANG, STATS) .NE.O)  GO  TO  52 

COMPUTE  HIGHER  ORDER  STATISTICS  FOR  BETA  DISTRIBUTION 
HOSBD  DOES  NOT  USE  STATS  INTERNALLY  WITH  PRINT 

CALL  HOSBD (0 , J , IOE , 10, S LOBE ) 

ANGSANG+ANGINC 

IF (DEBUG. NE. 01  GO  TO  5056 

COMPUTE  AND  PRINT  LOBE  STATISTICS  IF  REQUIRED 

IF (KLOBE . LE . 0)  GO  TO  5550 
IGRPRsO 
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5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5430 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5560 

5590 

5600 

5610 

5620 

5630 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5760 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

3900 

5910 

5920 

5930 

5940 

5950 

5960 

3970 

5980 


WRITE ( I0L’ 1 )  I GRPR , KL0BE .GLINT, NRCS , NLOBE , THETA . ALPHA . DELTA , 

1  BETA , SCALE , I ZETA , I PCODE . WLEN , LBUG 
CLOSE (UN  I T: I OL ) 

CALL  SPAWNS 

OPEN  I  UN  IT: IOL .NAME: ’ LOBE . TMP ’ , TYPE* ’ OLD ’ ,REC0RDSIZE=181 , 

2  ACCESS:' DIRECT’ .SHARED) 

C 

5550  I F (GOULD)  WRI TE ( IOP , 5052) 

5052  FORMAT  I ’1') 

5056  IFIGOULD. AND. DEBUG. EQ.O)  CALL  DETACH (10, IER) 

IF (GOULD. AND. DEBUG. EQ.O)  GOULD: . FALSE . 

IF  I LP. AND. .NOT.LPF)  CLOSE (UN  I Ts IOP , D ISPOSEs ’ PR  I  NT’ ) 

IF (LP.AND. LPF)  CLOSE (UN I T= IOP , D I SPOSE* ’ SAVE’ I 
IF (LP)  CALL  ASNLUNI IOP, ’GD’ ,0) 

I F ( CMD . EQ . CMOS (10))  GO  TO  3 
CMD=CLEAR 
GO  TO  70 
C 

C  "PRYT”  -  MASTER  ROTATION  ANGLES 

C 

55  CONTINUE 

222  FORMAT ( ' O' / ( 1 X ,3F8 .31) 

C  COMPUTE  YAW=>RESULT 
PRYTV(3)=ARR(3) 

C=CO$(ARR (3 )#RAD) 

S  =  S I  N  ( A R R  ( 3  HtRAD ) 

CALL  MATPUTIPRYT.1 . ,0. ,0. . 0 . , C , -S ,0 . . S ,C) 

C  MULTIPLY  BY  ROLL=>ROLLMRESULT:RESULT 
PRYTV ( 2 ) =ARR 1 2 ) 

C=COS ( ARR ( 2 ) *RAD 1 
S=S I N ( ARR (2 1 MRAD) 

CALL  MATPUT(MAT,C,-S,0. , S,C,0. ,0. ,0. ,1 . ) 

CALL  MATMUL (MAT, PRYT, PRYT) 

C  MULTIPLY  BY  PITCH=>PITCH*AESULTsRESULT 
PRYTV ( 1 1 =ARR ( 1 1 
C=C0S(ARR(1 IMRAD) 

S-SINIARRI1 )NRAD) 

CALL  MATPUTIMAT.C.O. , S , 0 . , 1 . ,0 . , -S , 0 . ,C) 

CALL  MATMULIMAT, PRYT, PRYT) 

C  MULTIPLY  BY  T I LT=>RESULT#TILT=RE5ULT 
PRYTV ( 4 ) :ARR (4) 

C=COS(ARRf 4)*RAD) 

S=SIN(ARR(4)mRAD) 

CALL  MATPUTIMAT.C.O. ,-S,0. ,1 . ,0. ,S,0. ,C) 

CALL  MATMUL (PRYT, MAT, PRYT J 
C  MULTIPLY  BY  B:>B*RESULT=RESULTs>PRYT 

CALL  MATPUT (MAT , 0 . ,-1.,0.,1.,0.,0.,0.,0.,1.l 
CALL  MATMULIMAT, PRYT, PRYT) 

GO  TO  3 
C 

C  "PAUSE”  -  WAIT  FOR  A  CARRIAGE  RETURN 

C 

60  CONTINUE 

READI I  OR , 6  I ,END  =  9990)  J 

61  FORMAT (A1) 

GO  TO  3 
C 

C  "RAMDB"  -  MODIFY  COMPONENT  RCS 

C 

65  RAM310.««IARRM)/10. ) 
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5990> 

GO  TO  3 

6000> 

C 

601 0> 

C 

"CLEAR”  -  CLEAR  RCS  VECTOR 

6020> 

C 

6030> 

70 

CONTINUE 

6040> 

1 PASS®0 

6050> 

OBJECT:. FALSE. 

6060> 

NVALIDsO 

6070> 

NT0TAL=0 

6080> 

00  71  Jsl.MAXRCS 

6090> 

RCS2IJ1S0. 

61 00> 

RCS(J) =0. 

6 1 1  0  > 

YYVIJ.1 1=0. 

6 1  20> 

YYVIJ ,21=0. 

6 1  3  0  > 

XV! J 1=0. 

6 1  4  0  > 

C 

EXTRA (J  J so . 

61 50> 

1 F ( OPENE )  WR 1 TE 1 1 OE ’ J ) DZE no , DZERO , DZERO , OZERO , DZERO , DZERO 

6 1  60> 

71 

CONTINUE 

61 70> 

IF ICMO.EQ. CLEAR)  QO  TO  3 

6 1  8  0  > 

C 

6 1  9  0  > 

C 

RESET  SYSTEM  VECTORS 

6200> 

C 

62 1  0> 

POLAR:) . 

6220> 

C 

6230> 

C 

SET  PYRT  TO  IDENTITY  BY  DEFAULT 

6240> 

C 

6250> 

DO  74  J=1 ,9 

6260> 

74 

PRYT ( J  J  =  1  DENT ( J I 

6270> 

C 

8280> 

C 

RESET  LOBE  PARAMETERS 

6290> 

C 

6300> 

KLOBE=0 

631 0> 

ILOBEsO 

6320> 

ML0BE=0 

6330> 

NLOBEs-3 

6340> 

BETA=0. 

6350> 

C 

6360> 

C 

RESET  RPRCS  PARAMETERS 

6370> 

C 

6380> 

IRPRsO 

6390> 

C 

6400> 

C 

RESET  VALID  VECTOR  AFTER  EACH  OBJECT 

641 0> 

C 

6420> 

72 

VALID! 1 1 =-360 . 

6430> 

VAL 1 D 1 2 1 =360 . 

6440> 

DO  73  J=3 , 6 

6450> 

73 

VAL ID! J) =0. 

6460> 

RAM=1 . 

6470> 

DO  75  4=1 ,9 

6480» 

75 

TRANS ( J )  =  1  DENT ( J ) 

6490> 

DO  76  J=1 ,3 

6500> 

76 

ORQN ( J ) =0 . 

6510> 

IFdPASS.EQ.1l  QO  TO  99 

8520> 

QO  TO  3 

6530> 

C 

6540> 

C 

DEBUG  OUTPUT 

6550> 

C 

6560> 

80 

1  BUG* IF  1 X 1 ARR 1 1 1 1 

6570> 

LBUQ= IF  IX (ARR 12 1 1 

6580> 

QO  TO  3 
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6990 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

8800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6860 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6970 

6980 

6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 


C 

C  NEW  WAVE  LENGTH 

C 

85  I  WLs I WL ♦ 1 

I F ( IWl.EQ.10)  GO  TO  9990 
0LDGRF=. FALSE. 

86  REWIND  I  OF 

I F (KLOBE . GT . 0)  CLOSE (UNITS I0L) 

I F  (GOULD I  CALL  DETACH! I0P  ,  IERI 
GOULD:. FALSE. 

LINEsO 
GO  TO  70 
C 

C  THE  FOLLOWING  COMMANDS  ARE  NEW 

C  ADDITIONS  TO  RCS  FOR  THE  LOBE 

C  STATISTICS 

c 

C  "L03E”  -  SET  LOBE  SWITCH 

C 

500  CONTINUE 

C  IF ( I  RPR . EQ . 1 . AND . ARR ( 1 i . EQ. 0 . )  GO  TO  3 

C  I F ( I  STATS .EQ. 1 .AND .ARR 1 1 1 . GT . 0 . )  GO  TO  573 

LLOBEsARRM) 

IFILLOBE.GT.O. AND. KLOBE. LE. 01  CALL  NOLOCKIIOL) 

IF (LLOBE.GT.O. AND. KLOBE. LE.OI 

1  OPEN(UNIT=IOL, NAMEs’ LOBE. TMP' , TYPES ’ UNKNOWN’ , RECORDS  I ZEsI 81 , 

2  ACCESS:’ DIRECT’ .SHARED) 

IFILLOBE.GT.O. AND. KLOBE. LE.O)  WRITE ( IOL *  1 1  DUMMY 
IFILLOBE.GT.O)  MLOBEsLLOBE 
iriLLOBE.GT.O)  KLOBE=LLOBE 
GO  TO  3 
C 

C501  WRITEI2.502) 

C502  FORMAT!’  RCS>  —  RPRCS  AND  LOBE  CANNOT  BOTH  BE  ACTIVE’) 

C  GO  TO  9993 

C 

C  -BETA’’  -  DEFINE  BISTATIC  ANGLE 

C 

510  BETAsARRM) 

IFIBETA.LT. 5.)  BETAsO. 

GO  TO  3 
C 

C  "ORGN”  -  COMPONENT  ORIGIN  VECTOR 

C 

520  ORGN(l)sARRII) 

ORGN  1 2 ) sARR ( 2 ) 

ORGN  1 3 ) sARR  1 3 ) 

GO  TO  3 
C 

C  "DELT”  -  FIRST  LOBE  ANGULAR  RANGE 

C 

540  DO  541  Jsl,4 

541  DELTA  I J ) sARR I J ) 

DELTA ( 2 ) sAMAXI (DELTA  1 1 ) , DELTA ( 2 ) > 

DELTA  13 JsAMAXI (1 . , DELTA  13 ) ) 

IF (DELTA (4) .LT.4.0.0R.DELTAI4) .GT.180. )  DELTA(4)s48. 

GO  TO  3 
C 
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7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 

7600 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 


C  "ALPH"  -  SECOND  LOBE  ANGULAR  RANGE 

C 

550  DO  551  Jsl ,  4 

551  ALPHA ( J ) =ARR ( J ) 

GO  TO  3 

C 

C 

C  END  OF  LOBE  COMMANDS 

C 

CtHtKMKNNNttXKXXXttttttKXXXKXXKXIHIKttMKtttt 

c 

C  "RPRCS"  -  RELATIVE  PHASE  RCS 

C 

560  CONTINUE 

C  I F IMLOBE . GT . 0 .AND . ARR 1 1 ) . GT . 0 . ) GO  TO  501 

C  IF ( I  STATS. EQ. 1 . AND. ARR 1 1 ) .GT. 0. )  GO  TO  571 

I RPR=MI NO  1 1 ,MAX1 10. . ARR ( 1 ) ) ) 

I F 1 1  RPR . EQ. 0)  GO  TO  3 
RPRCSG=ARR(2) 

LRPRCS=ARR (3 1 
IFILRPRCS.LE.il  LRPRCS=5 
LOBE-1 
LLOBEsi 
MLOBEst 
GO  TO  3 
C 

C  "STATS"  -  HIGHERER  ORDER  STATISTICS 

C 

570  CONTINUE 

C  IF  1 1  RPR . EQ . 1 .AND .ARR ( 1 1 .GT . 0 . 1  GO  TO  571 

C  IFIMLOBE.GT.O. AND .ARRI1 I . GT . 0 . I  GO  TO  573 

I STATS=M I  NO  1 1 ,MAX1 10. .ARR II  III 
STATS-ARRI2) 

I F (STATS . LE . 0 . I  STATSslO. 

I F ( I  STATS . EQ. 0 .OR .OPENE I  GO  TO  3 

OPEN (UN  I T= IOE , NAME= ’ RCSTMP . DES’ , TYPE* ’ SCRATCH’ . 

1  ACCESS: ’ D I RECT ’ , RECORDS  I ZE - 1 6 , MAXREC  =  1 81  I 

OPENE:. TRUE. 

C 

DO  575  J=1 ,181 

575  WR I TE I IOE ’ J I DZERO, DZERO , DZERO , DZERO , DZERO, DZERO 

GO  TO  3 
C 

C571  WR I TE ( 2 , 572  I 

C572  FORMAT! *  RCS>  -•  RPRCS  AND  STATS  CANNOT  BOTH  BE  ACTIVE’ I 

C  GO  TO  9993 

CC 

C573  WR I TE (2 , 574) 

C574  FORMAT ( '  RCS>  ->  LOBE  AND  STATS  CANNOT  BOTH  BE  ACTIVE’) 

C  GO  TO  9993 

C 

C  "GLINT"  -  SET  GLINT  SWITCH 

C 

580  GL  I  NT:M  I  NO  1 2  ,MAX1  (0.  .ARR  Ml  I) 

GO  TO  3 
C 

C  "NEWTH"  -  REPEAT  LOOP  FOR  NEW  THETA  ANGLE  RANGE 

C 

590  CALL  CNEWTHIKLOBE, ROUTE) 
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7790> 

7800» 

C 

7810> 

C 

7820» 

c 

7830’ 

800 

7840> 

601 

7850> 

7860> 

7B70> 

c 

7880> 

C 

7890> 

C 

7900> 

90 

791 0> 

7920> 

7930> 

7940> 

91 

7950> 

C 

7960> 

7970> 

92 

7980> 

7990> 

C 

8000> 

C 

801 0> 

C 

8020’ 

C 

8030’ 

C 

8040’ 

C 

8050> 

C 

8060’ 

C 

8070’ 

100 

8080> 

C 

8090> 

8100’ 

99 

8110’ 

8120’ 

8 1 3  0  > 

81 40> 

8 1 5  0  > 

8160’ 

81 70> 

8180’ 

8190’ 

8200> 

C 

8  2 1  0  > 

8220’ 

8230’ 

8240’ 

8250’ 

C 

8260’ 

C 

6270’ 

c 

8280’ 

c 

6290’ 

c 

8300’ 

c 

8310’ 

c 

8320’ 

c 

8330’ 

c 

6340’ 

c 

8390’ 

c 

8360’ 

c 

8370’ 

101 

8380’ 

GO  TO  (3,86,99931 .ROUTE 

”QS I ZE"  -  SET  GRAPH  SIZE 

DO  601  J  =  1 ,5 
QSIZE (Jl-ARRIJ) 

CALL  GRPICIGSIZEI1 ) ,QS I ZE ( 2 ) .QSIZE (3) .QSIZE (41 ) 
GO  TO  3 

DETERMINE  IF  COMMAND  IS  AN  OBJECT 

CONTINUE 

DO  91  J:i , NOB J 
IOBJSJ 

I F (CMD . EQ.OBJS1 J 1 1  GO  TO  100 
CONTINUE 

WRITE! IOW, 92)  CMD 

FORMAT C  RCS>  —  INVALID  FUNCTION  ,A4, ’ "* ) 

GO  TO  9993 


ADD  OBJECT  CONTRIBUTION  TO  RCS  VECTOR 


IFIWLEN.LE.O. )  GO  TO  9991 

NOTE:  IPASS  IS  ONLY  USED  WITH  RFLAT  IN  CONJUCTION  WITH  LOBE 

I PASS=-1 
IPASS=IPASS*1 
OBJECTS. TRUE. 

ANQsTHETA ( 1 ) 

ANGMAXsTHETA ( 2 ) *  .  001 
AMGINCSTHETAI3) 

CPHR=COS ( PHRNRAO ) 

SPHRsS IN(PHRMRAD) 

IRCS=1 

CALL  TRANSP (PRYT .PRYT2 I 
CALL  MATMULIPRYT2, TRANS, COMP) 

IF  I ( I  BUG. AND. 2) .EQ.O)  GO  TO  101 
WRITEI2.222)  TRANS 
WRITEI2.222)  PRYT 
WRITE(2,222)  COMP 

ANG  IS  RADAR  THETA  ANGLE 
PHR  IS  RADAR  PHI  ANGLE 

IF  END  OF  ANGLE  RANGE  FOR  OBJECT,  OUTPUT  PARTIAL 
LOBE  DATA  IF  REQUIRED,  THEN  GET  NEXT  COMMAND. 

KLOBEs>$ET  FROM  LOBE  ENTRY  POINT 
MLOBE  =  > LOBE  IS  SET  FOR  RUN 
LLOBEs>LOBE  IS  SET  FOR  OBJECT 
LOBE s> LOBE  IS  SET  FOR  CURRENT  THETA 

IF  1 1 RCS . LE . NRCS)  QO  TO  131 

IFIIRPR.EQ.1 .AND. KLOBE. EQ.O)  GO  TO  13$ 
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8300 

8400 

8410 

8420 

8430 

8440 

8450 

8460 

8470 

8480 

8490 

8500 

8510 

8520 

8530 

8540 

8550 

8560 

8570 

8580 

8590 

8600 

8610 

8620 

8630 

8640 

8650 

8660 

8670 

8680 

8690 

8700 

8710 

8720 

8730 

8740 

8750 

8760 

8770 

8780 

8790 

8800 

8810 

8820 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8920 

8930 

8940 

8950 

8960 

8970 

8980 


I F I L LOBE . LE . 0 )  GO  TO  72 
C  OUTPUT  LOBE  INFO  TO  TEMPORARY  FILE 
NLOBEsNLOBEM 
00  139  J s 1  ,4 

139  WRITE) I0L ’ NLOBE* J )  (SLOBI (J , I , 1 ) , I =1 ,NRCS) 

DO  133  1=1  ,  NRCS 

I F ( SLOBE  1 1 , 1 ,2) .NE.O. )  GO  TO  134 

133  CONTINUE 
GO  TO  135 

134  NLOBE=NLOBE+4 
DO  137  J:1  ,4 

137  WRITE) IOL’NLOBE*JI  (SLOBE I J , I , 2 ) , I *1 , NRCS) 

C 

C  MAKE  SECOND  PASS  THROUGH  RFLAT  FOR  LOBE 

135  t  F (CMD . EQ.OBJS I 1 3 ) .AND . I  PASS . EQ . 0)  IPASS  =  1 
GO  TO  72 

C 

C  DETERMINE  IF  VALID  ANGLE  FOR  LOBE  CONTRIBUTION 

C 

131  LOBE=LLOBE 

IF (LOBE . LE .0)  GO  TO  132 
SLOBE II , IRCS,1)sO. 

SLOBE (1 , IRCS.2IS0. 

C 

C  SKIP  ANGLE  IF  NOT  VALID  FOR  CURRENT  OBJECT 
C 

132  DO  102  J:1  ,5,4 

IF (VALIDIJ) -VALIDIJ+1 ) I  103,102,104 

103  IFIANG. LT. VALIDIJ) .OR. ANG .GT . VAL I D I J ♦ 1 1 1  GO  TO  102 

IFfVALID(J+2)-VALID(J*3) )  105,109,106 

105  IFIPHR.LT. VALIDIJ*2) . OR . PHR . GT . VAL I D I J *3 ) )  GO  TO  102 

GO  TO  109 

106  IF(PHR.GT.VALID(J+3) . AND .PHR.LT .  VAL I D I J*2 ) )  GO  TO  102 

GO  TO  109 

104  IFIANG. GT.VALID(J*1) .ANO.ANG.LT. VALIDIJ))  GO  TO  102 

IFIVALIDIJ+2) -VAL I D I J*3  )  I  107,109,108 

107  IFIPHR.LT. VALIDIJ+2) . OR . PHR . GT . VAL I D I J +3 ) )  GO  TO  102 

GO  TO  109 

108  IF(PHR.GT.VALID(J*3) .AND. PHR.LT. VALIDIJ*2) )  GO  TO  102 

GO  TO  109 
102  CONTINUE 
C 

C  CURRENT  ANGLE  IS  NOT  VALID.  INCREMENT  ANGLE. 

C 

GO  TO  130 
C 

C  CURRENT  ANGLE  IS  VALID.  FIND  PRIMED  ANGLES  FROM  RADAR  COORDINATES. 
C  NOTE!  THIS  PROGRAM  STORES  MATRICES  BY  ROWS  (NOT  COLUMNS) 

C 

100  CONTINUE 

CANGsCOS (ANGffRAD) 

SANGSSINIANGKRAD) 

CTHP  =  COMP I  3 ) *SANQMCPHR  +  COMP 1 6  I *SANG*SPHR+COMP 1 9 ) *CANG 
THP=ACOSK(CTHP) 

STHPsSINITHP) 

CPHPS= (COMP  1 1 )«SANG#CPHR+COMP  <4)*SANG#SPHR*COMP 17) *CANG) 

SPHPSs (COMP (2)*SANG*CPHR+COMP (5)*SANG«SPHR*COMP (8)«CANQ) 

PHPsO. 

IFIABSISPHPS) .LT. IE-5. AND. ABSICPHPS) .LT.1 .E-5)  GO  TO  1111 
PHP*ATAN2 I SPHPS ,CPHPS) 
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89«0> 

1111 

CPHPsCOSIPHP) 

9000> 

C2PHP*CPHP#CPHP 

901  0> 

SPHPsSINIPHP) 

9020> 

S2PHP=SPHP»SPHP 

9030> 

223 

FORMAT (4F8. 31 

9040> 

C 

9050> 

C 

9060> 

C 

9070> 

C 

9080> 

C 

9090> 

C 

ADO  OBJECT  CONTRIBUTION  FOR  CURRENT  ANGLE 

91  00> 

C 

9 1 1 0  > 

C 

91 20> 

C 

91 30> 

RCSC=0 . 

91  40> 

IFII0BJ.LE.6.0R. IOBJ.GE.10)  GO  TO  110 

91  50> 

C 

91  60> 

C  FIND  UNPRI  1ED  ANGLES  FOR  OBJECTS  WHICH  USE  POLARIZATION 

91 70> 

C 

91 80> 

CTH  =  Prt  r2(3)KSANG»CPHR«PRYT2(6)ftSANGMSPHR»PRYT2(9J«CANG 

91 90> 

THsACOSK(CTH) 

9200> 

STHaSINITH) 

9  2 1 0  > 

CPHS=(PRYT2I1 )ft$ANGNCPHR+PRYT2 (4)*SANG*SPHR*PRYT2(7)ftCANQ) 

9220> 

SPHSa (PRYT2 ( 2 J#SANG#CPHR+PR YT2 ( 5 ) «SANG«SPHR*PR YT2 ( 8 1 «CANG ) 

9230> 

PH*0 . 

9240> 

IF(ABSISPHS) . LT . 1 . E-5 . AND . ABS (CPHS) .LT.l.E— 51  GO  TO  1112 

9250> 

PH:ATAN2(SPHS,CPHSI 

9260> 

1112 

CPHsCOSIPH) 

9270> 

SPHsSINIPH) 

9280> 

IF  1 1 (BUG. AND. 4) . EQ.Oi  GO  TO  110 

9290> 

DTHaTH/RAD 

9300> 

DPHrPH/RAD 

931 0> 

WRITEI2.223J  DTH.DPH 

9320> 

C 

9330> 

C 

GO  TO  SPECIFIED  OBJECT 

9340> 

C 

9350> 

110 

GO  TO  (200,210.220,230,240,250.260,270,280,290.300.310,320,330. 

9360> 

9  340,350,360,3701 ,  IOBJ 

9370> 

C 

9380> 

C 

9390> 

C  " 

OGIVE"  -  OGIVE  A  TRUNCATED  OGIVE 

9400> 

200 

CALL  OGIVE 

941 0> 

GO  TO  120 

9420> 

C 

9430> 

C  * 

CYLIN’  -  ELLIPTICAL  CYLINDER 

9440> 

C 

INPUTSaA.B.L 

94S0> 

210 

CALL  CYLIN 

9460> 

GO  TO  120 

9470> 

C 

9480> 

c  - 

CONE"  -  TRUNCATED  ELLIPTICAL  CONE  A  CIRCULAR  CONE  A  STANDARD  CONE 

9490> 

220 

CALL  CONE 

9500> 

GO  TO  120 

951  0> 

C 

9520> 

C  " 

PARAB"  -  PARABOLA 

9530> 

230 

CALL  PARAB 

9540> 

GO  TO  120 

9550> 

C 

9560> 

C  * 

SOLID-  >  ELL IPSIOD  A  SPHEROID 

9570» 

C 

INPUTSaA.B.C 

9580> 

240 

CALL  SOLID 
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9590> 

60 

TO  120 

9600> 

C 

961 0> 

C  "TORUS" 

-  TORUS 

9620> 

250  CALL 

TORUS 

9630> 

GO 

TO  120 

9640> 

C 

96S0> 

C  "WE0GE” 

-  WEDGE 

9660> 

260  CALL 

WEDGE (TRANS, CTH , STH ,CPH , SPH) 

9670> 

60 

TO  120 

9680> 

C 

9690> 

C  "WIRE" 

WIRE 

9700> 

270  CALL 

Wl  RE (TRANS, CTH, STH, CPH.SFH) 

9710> 

60 

TO  120 

9720> 

C 

9730> 

C  “LWIRE" 

-  LOOP  WIRE 

9740> 

280  CALL 

LWI RE (TRANS, CTH, STH, CPH. SPH) 

9750> 

GO 

TO  120 

9760> 

C 

9770> 

C  "CAVA” 

-  CAVITY  AREA 

9780> 

290  CALL 

CAVA 

9790> 

GO 

TO  120 

9800> 

C 

9  8 1  0  > 

C  "CAVD" 

-  CAVITY  DIAMETER 

9820> 

300  CALL 

CAVD 

9830> 

60 

TO  120 

9840> 

C 

9C50> 

C  "CFLAT" 

-  CIRCULAR  FLAT  PLATE 

9860> 

310  CALL 

CFLAT 

9870> 

60 

TO  120 

9880> 

C 

9890> 

C  "RFLAT" 

-  RECTANGULAR  FLAT  PLATE 

9900> 

320  CALL 

RFLATI IPASS) 

991  0> 

IF  I  IPASS.EQ.2)  RCSC  =  0. 

9920> 

60 

TO  120 

9930> 

C 

9940> 

C  "TWAVE" 

-  TRAVELING  WAVE 

9950> 

330  CALL 

TWAVE 

9960> 

60 

TO  120 

9970> 

C 

9980> 

C  “0 1 HED” 

-  DIHEDRAL 

9990> 

340  CALL 

DIHED 

10000> 

60 

TO  120 

10010> 

C 

1 0020> 

C  "TCORN" 

-  GENERA'.  TRIANGULAR  CORNER 

1 0030> 

350  CALL 

TCORN 

1 0040> 

60 

TO  120 

1 0050> 

C 

1 0060> 

C  "RCORN” 

-  GENERAL  RECTANGULAR  CORNER 

1 0070> 

360  CALL 

RCORN 

1 0080> 

60 

TO  120 

1 0090> 

C 

10100> 

C  "ECORN” 

-  GENERAL  ELLIPTICAL  CORNER 

1 01 1 0> 

370  CALL 

ECORN 

1 01 20> 

60 

TO  120 

1 01 30> 

C 

1 01 40> 

C 

1 0 1 5  0  > 

C  ADD 

DCS  CONTRIBUTIONS  TO  VECTORS 

1 0 1 6  0  > 

C 

1 0 1 7  0  > 

120  CONTINUE 

1 01 80> 

IFIRCSC.LT. 0. 1  WRITE  1 I0W, 121)  CMD.ANG.RCSC 
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1 01 90* 

121 

FORMAT  1 ’  RCS*  --  ’ , A4 , F7 . 1 , El 5 . 6 ) 

1 0200> 

1 F  (RCSC . GT . CUTOFF  I  NVAL 1 D=NVAL 1 0*1 

1 02 1 0  > 

RCSC=RCSCttRAM 

1 0220  > 

IFIIRPR.EQ.il  GO  TO  136 

1 0230> 

RCS (  IRCSIsRCSI  ircsurcsc 

1 0240> 

RCS2 ( 1 RCS) -ACS2 ( 1 RCS) * RCSC* RCSC 

1 0250> 

136 

IF  1 ISTATS.EQ.O)  GO  TO  130 

1 0260> 

C 

1 0270> 

READ ( 1 OE ’ 1 RCS )  T1 , T2 , T3 , T4, T5 . T6 

t  0280> 

T1  =T1 +DBLE l RCSC) 

1 0290> 

T2=T2+DBLE(RCSCJ**2 

1 0300> 

T3:T3+DBLE(RCSC)N*3 

1 03 1 0> 

T4=T4+DBIE (RCSC)**4 

1 0320> 

SQs SORT (RCSC) 

1 0330> 

T5sT5+SQ 

1 0340> 

T6=DMAX1 ( T  6 , SQ ) 

1 03S0> 

VR 1 TE ( 1 OE ’ 1 RCS)  T1 , T2 , T3 , T4 , T5 , T6 

1 0360> 

C 

I0370> 

C 

INCREMENT  ANGLE 

1 0380> 

C 

1 0390* 

130 

CONTINUE 

1 0400> 

ANG=ANG*ANGINC 

1 041 0> 

1 RCS  = IRCS+1 

1 0420> 

IF  1 1  PASS. HE .2)  NTOTALsNTOTAL*1 

1 0430* 

GO  TO  101 

10440> 

C 

1 0450> 

C 

1 0460> 

9991 

WRITE! IOW, 9992) 

1 0470> 

9992 

FORMAT  I’  RCS*  --  WAVELENGTH  UNDEFINED’) 

1 0480> 

C 

10490* 

9993 

WRITE! IOW, 9994)  LINE 

1 0500> 

9994 

FORMAT! ’  RCS*  ->  LAST  DATA  BASE  LINE  READ  IS’,14) 

10510* 

C 

10520* 

9990 

CONTINUE 

10530* 

C 

10540* 

C 

TERMINATE  RCS  PROGRAM 

10550* 

C 

10560* 

IFIKLOBE.GT.O)  CLOSE (UN  1 T= 10L , D 1 SPOSE* ’ DELETE ’ ) 

10570* 

IF(OPENE)  CLOSE (UNI T  = 1 OE ) 

10580* 

IF (GOULD)  CALL  DETACH! IOP , IER) 

10590* 

T IME=SECNDS ( T 1 ME ) / 6  0 . 

10600* 

1 F ( IZETA.EQ.O)  CALL  QRMODE ( ’ ALPHA ’ , 3 ) 

10610* 

WR 1 TE (2 , 9998 )  TIME 

10620* 

9998 

FORMAT!’  RCS*  --  ELAPSED  TIME  =’,F5.1,’  MINUTES’) 

10630* 

C 

10640* 

9999 

CALL  EXIT 

10850* 

END 
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>>>  FILE  s 
10>  C 
20> 

30>  C 
40> 

50>  C 
60> 

T0> 

80> 

90> 

100>  C 
1 1 0> 

1 20> 

1 30> 

1 40> 

1 50> 

1 6  0  > 

1 70> 

1  80>  C 
1  9  0  >  1 
200> 

21 0> 

220> 

230>  C 
240>  591 
250>  592 
260> 

270> 

280> 

290> 


SY  0 : C 1 1 1 1 1 23RCSCMD . FTN  27-3EP-83  14:34:38  <<< 

SUBROUTINE  CNEWTHIKIOBE, ROUTE) 

I NTEGERM2  ROUTE 

INTE0ERM2  DEBUG, COMB  I N.GLI NT 

LOG ICAL*2  GOULD , IP , OPENE , LPF , OLDGRF .MGRAPH , GDAUTO 
COMMON/SWITCH/  IOW,THETA(4) , GOULD, LP, OPENE, LPF, OLDGRF , MGRAPH, 

1  GDAUTO, DEBUG, COMB  IN, GLINT, IBUG.LBUG, ISTATS, IZETA 

ROUTED 

IF (THETA (4) .GT. THETA (2) )  GO  TO  1 

IF ( .NOT. MGRAPH)  RETURN 

IF ( IZETA. EQ.O. AND. .NOT. GDAUTO)  RETURN 

IF  I IZETA. EQ.O)  CALL  GRTEKG 

IF « IZETA. NE.O)  CALL  GRPLOT 1 1 1 . , 8 . 5 , -1 ) 

RETURN 

IFIKLOBE. EQ. 1. OR. GLINT. EQ.1. OR. ISTATS. EO.1. OR. DEBUG. NE.O)  GO  TO  591 
OLDGRF:. TRUE. 

ROUTED 

RETURN 

WRITE! IOW, 592 ) 

FORMAT! ’  RCS>  --  MULTIPLE  GRAPH  SEGMENTS  ARE  NOT  ALLOWED  WITH:’/ 

1  LOBE,  STATS,  OR  DEBUG!’) 

ROUTED 

RETURN 

END 
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>>>  FILE  -  SYO : [ 1 1 1 , 1 2 ] RCSOBJ . FTN  27-SEP-83  14:34:44  <<< 

0>  C 

5>  C  <tt><tt><«><«><«><M><«><M><*><*><*><it><N><ll><N><»><lt><N><N> 

10>  SUBROUTINE  OGIVE 

20>  C 

30>  C  WITH  LOBE 
40>  C 

50>  IMPLICIT  REALM  (A-Z) 

60>  INTEGERS  LOBE 

70>  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP . SPHP , 

80>  S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

90>  S  R1 , A , B , BETA 

100>  C 

1 1 0>  THPTsTHP 

1 20>  IFITHP.GT.PI2)  THPTsPI-THP 

1 3 0>  IF (LOBE .GT.O)  L2  =  SQRT (R1ttR1-(B+R1 -A) #*2) 

1 40>  IF IB.GT.O. )  GO  TO  S 

1 5 0 >  3  ALP  =  ACOSK(1 .-A/R1 1 

1 60>  RCSC=R1»R1#SIN(ALP)##2/  14.MP I ) 

1 7  0  >  NMA  =  P 12-ALP 

1 80>  IFITHPT.LT. NMA)  GO  TO  2 

190>  C 

200>  C  90-ALP  <  THETA’  <  90 

21 0>  S IGMA=P l«R1ttR1tf( 1 . - (R1 -A) / IRImSTHP) I 

220>  I F ( LOBE . LE . 0)  GO  TO  4 

230>  IFIRCSC.GT. SIGMA)  GO  TO  6 

240>  CALL  TOLOBE ( 1 , S I GMA . AtfCPHP , AffSPHP , 0 . I 

245>  4  RCSCsSIGMA 

250>  RETURN 

260>  C 

270>  C  THETA’  s  90-ALP 

280>  6  C A L P 2  =  ( 1  .  -  A/R1  lint 2 

290>  SALP  =  S I N I  ALP ) 

300>  ZALP=R1#SALP*( CALP2-SQRT (CALP2MCALP2-SALPMSALP ) ) 

31 0>  WALP=ZALP*TAN( NMA) 

320>  CALL  TOLOBE M , RCSC, WALPNCPHP , WALPNSPHP . ZALP ) 

330>  RETURN 

360>  C 

370*  C  THETA’  «  90-ALP 

380>  2  Zsl6.»PmCOS(THPTHn»6«(1  .-TANIALPHH»2»TANITHPTHh»2HHI3 

390>  IFILOBE. GT.O. AND. Z.LT. CUTOFF)  GO  TO  6 

400>  IFIZ.LT. CUTOFF)  RETURN 

41 0>  S  IGMA=WLEN#WLENl»TAN  ( ALP)#»4/Z 

420>  IFILOBE. LE.O)  GO  TO  1 

430>  IFIRCSC.LT. SIGMA)  GO  TO  6 

440>  CALL  TOLOBEI1 .SIGMA, 0. ,0. ,L2) 

450>  1  RCSC=AMIN1 I RCSC .SIGMA) 

480>  RETURN 

470>  C 

480>  C  TRUNCATED  OGIVE  (BETA  IS  ALWAYS  <  ALPHA) 

490>  C 

S00>  5  BET=BETA«RAD 

51 0>  IF (THPT .GT . P 12-BET)  GO  TO  3 

520>  C 

530>  C  THETA’  s  o 

540>  RCSC=P lNBNBttTANIBET)tf*2 

5S0>  IFISTHP.LT. CUTOFF. AND. LOBE. GT.O)  CALL  TOLOBEI1 .RCSC.O. ,0. ,L2) 

560>  IFISTHP.LT. CUTOFF)  RETURN 

570>  C 

380>  C  BETA  <  THETA’  <  90-BETA 
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590> 

600> 

6 1 0> 
620>  C 
630>  C 
640> 
650> 
660> 
670> 
680> 
690>  C 
700>  C 
705>  C 
71 0> 
720>  C 
730>  C 
7 40>  C 
750> 
760> 
770> 
780> 
790> 
800>  C 
8 1  0  > 
820> 
830> 
840> 
850> 
860>  C 
870> 
880> 
890>  C 
900>  C 
9 1 0> 
920» 
930> 
940> 

9  50> 
960> 
970>  C 
980>  C 
990>  1 
1000> 

1  01 0> 

1 020>  C 
1 030>  C 
1 03 5  >  C 
1 040> 

1 050>  C 
1 060>  C 
1 070>  C 
1 080> 

1 090> 

1 1  00» 

1 1 1 0> 

1 1  2  0  > 

1 1 30>  C 
1 1  40> 

1 1  50> 

1 1 60» 


S I GMA=WLEN»B«TAN ( TH»T*BET 1  XX2/ ( 8  . XP IXSTHP 1 

IF (LOBE .GT.OI  CALL  T0L0BE ( 1 , SIGMA , BxCPHP , BxSPHP , L 2 ) 
I F (BET . LE . THPT . AND . THPT . LE . P 1 2-BET)  GO  TO  1 

THETA’  <  BETA  OR  THETA’  >  90-BETA 

SIQsWLENXBXTAN(THPT-BET)XX2/(8.XPlXSTHP) 

I F I  LOBE . GT . 0)  CALL  TOLOBE 12 , SIG .BxCPHP , BxSPHP , -L2 ) 

S IGMA-S IGMA+S IG 
GO  TO  1 
END 


<XXKXKXXXKXXXKXXXXXXXXXXXXXKXXXXXK> <XXX> 

SUBROUTINE  CYLIN 
WITH  LOBE 

IMPLICIT  REALX4  (A-Z) 

INTEGERM2  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN . P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

8  C2 P HP , S2PHP , RAD , POLAR , LOBE , 

9  A.B.L 

IFISTHP.LE.O. )  RETURN 

A2=A*A 

B2?BxB 

Z*  IA2xCPHPxCPHP+B2xSPHPXSPHP)XXl  . 5 
RCSC=2.xPlxLXLXA2xB2/(WLENXZ) 

0*8  .XP  IXCTHPXCTHPXZ 
IFIABSIOI .LT. CUTOFF)  GO  TO  1 

THETA’  •  90 

SIGMA=2.xWLENxA2xB2xSTHP/0 
RCSC=AMIN1 (RCSC, SIGMA) 

I F I  LOBE . LE . 0 )  RETURN 

CALL  TOLOBE 1 1 , RCSC/2 . , AxCPHP , BxSPHP , 0 . ) 

CALL  TOLOBE (2 , RCSC/2, , AnCPHP ,B#SPHP , L ) 

RETURN 

THETA’  =  90 

I F ( LOBE . GT . 0)  CALL  TOLOBE ( 1 , RCSC , AttCPHP , BttSPHP , L/2 . ) 

RETURN 

END 


OtXftXKXMXttXMXXXttXNXftXttXlIXftXttXttXftXlIXftXIO 
SUBROUTINE  CONE 

WITH  LOBE 

IMPLICIT  REALM  (A-Z) 

INTEGERS  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 
*  C2PHP , S2PHP , RAD , POLAR , LOBE , 

9  ALP.NU.A.L1 ,L2 

I  F (CTHP . GE . 0  .  . AND . S I N ( ALPxRAD ) .GE.STHP)  RETURN 

NU2:NU«NU 

TALP *TAN (ALPxRAD) 
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1 1 7  0>  D  =  SQRT ( S2PHP +NU2#C2PHP ) 

1 1 80>  POLAHQ=ATAN  ( -  MU/  I TALPMD)  I 

I190>  RCSC =8 . NP ltt(L2ttM1 . 5 - L 1 «« 1  , 5 ) «#2#TALP**4 

1200>  RCSC=RCSC/(9 .«WLEN*NU2ftABS(COS(POLANG)«*3) ) 

1 21 0>  Z=ISTHP*TALP»D*NU*CTHP)##2 

1 220>  Z*Z*8 .#P l#STHP*D 

1 230>  IF (Z .GT .CUTOFF )  GO  TO  1 

1 240>  I F ( LOBE . LE . 0)  RETURN 

1 250>  L= ( LI +L2 ) /2 . 

1260>  CALL  TOLOBE ( 1 , RCS , L#TALP#CPHP , L#NU#TALP#SPHP , L ) 

1 270>  RETURN 

1 280>  C 

1 290>  C  OFF  -  NORMAL  INCIDENCE 
1 300>  1  SIGMA: I (STHP/DI - <CTHP»TALP/NU))»*2 

1 3 1 0  >  S IQMA?S IGMA*WLEN*NU2ttNU*TALP/Z 

1 320>  RCSCsAMINI (RCSC , $ IGMAal LI +L2)) 

1 330>  I F I  LOBE . LE . 0)  RETURN 

1 340>  CALL  TOLOBE 1 1 , AM INI IRCSC/2 . , L1*S IGMA) , LlMTALPttCPHP , 

1  350>  1  LlttNUttTALPttSPHP , LI  I 

1 360>  CALL  TOLOBE 1 2, AM  INI IRCSC/2. ,L2»SIGMA) , L2*TALP»CPHP . 

1 370>  1  L2MNUNTALPNSPHP  ,  L2 1 

1380>  RETURN 

1 390>  END 

1 400>  C 
1 4 1 0  >  C 

141 5 >  C  <tt><»><N><»x«><»><M> <*><«> <«><*><«> <*><»><«> <*><«> <»><*> 

1 420>  SUBROUTINE  PARAB 

1 430>  C 

1 440>  C  WITH  LOBE 
1 450>  C 

1 460>  IMPLICIT  REALM  (A-ZI 

1470>  INTEGERS  LOBE 

1 480>  COMMON /RCSCOM/CUTOFF , WLEN . P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

1 490>  8  C2PHP.S2PHP, RAD, POLAR, LOBE, 

1 S00>  8  P 

1 51 0>  C 

1 5 2 0 >  C  RCS  IS  UNDEFINED  AS  THP  APPROACHES  PI 

1530>  IF IABSITHP) .GT.70.MRAD)  RETURN 

1 540>  RCSC=4.*PI*P*P/CTHP«M 

1 550>  IF  I  LOBE . LE . 0)  RETURN 

1 560>  TTHP=TAN ITHP I 

1 570>  CALL  TOLOBE 1 1 . RCSC , 2 . *P*TTHP*CPHP , 2 . »P»TTHP»SPHP , -P«TTHP»TTHP ) 

1 580>  RETURN 

1 590>  END 

1 600>  C 
1 6 1  0  >  C 

1 6 1 5 >  C  <»><!»» <»><«><«><«> <«><«> <*><«> <«><«><*><«>  <«><<!> <«><«><«> 

1 820»  SUBROUTINE  SOLID 

1 630>  C 

1 840»  C  WITH  LOBE 
1 850>  C 

1 660>  IMPLICIT  REALP4  (A-Z) 

1 6 7 0 >  INTEGERM2  LOBE 

1 680>  COMMON /RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

1 690>  8  C2PHP.S2PHP, RAD, POLAR, LOBE, 

1700>  8  A , B ,C 

1710>  C 

1 72  0  >  RsAMBNC 

1  730>  D: I AmSTHPNCPHP  MH*2 

1740>  DsD*IB«STHP»»SPHPHH»2 
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1 7$0>  D=D+(CNCTHP)»»2 

1  7  6  0  >  RCSC  =  PHtR»R/(D*D) 

1 770>  IF (LOBE .QT.O)  CALL  TOLOBE ( 1 , RCSC, A»STHP*CPHP , B»STHP*SPHP , 

1 780>  1  C»CTHP) 

1790>  RETURN 

1 800>  END 

1 81 C>  C 
1  8  2  0  >  C 

1  825 >  C  <»><»><*><*><*><»> <*><*> <*><*> <»><#> <*><#> <»><*> <*><*><»> 

1 830>  SUBROUTINE  TORUS 

1 840>  C 

1 8 5 0 >  C  WITH  LOBE 
1 860>  C 

1 870>  IMPLICIT  REALM  (A-Z) 

1 8 8 0 >  I NTEQER*2  LOBE 

1 890>  COMMON/RCSCOM/CUTOFF ,WLEN,PI,PI2, RCSC , THP , CTHP , STHP , PHP , CPHP . SPHP . 

1 900>  8  C2PHP.S2PHP, RAD, POLAR, LOBE, 

1 9 1 0 >  *  A.B 

1920>  C 

1 930>  RCSC =8 .  NP l»*3*BtfAMA/WLEN 

1 940>  IFISTHP.LT. CUTOFF. AND. LOBE. QT.O)  CALL  TOLOBEI1 .RCSC.O. ,0. ,0. ) 

1 950>  IFISTHP.LT. CUTOFF)  RETURN 

1 960>  S I Q1 sP I*  (BttA/STHP  +  BwB ) 

1  9  7  0  >  S  I G  2  =  0  . 

1 9  8  0  >  I F I  LOBE .QT.O)  CALL  TOLOBE 1 1 , AM) HI (RCSC . S 101 ) . A»CPHP , BWSPHP , 0 . ) 

1 990>  IF (ABS(CTHP) .LE.B/I2.NAI)  GO  TO  1 

2000>  SIG2=Plw (BwA/STHP-BftB) 

201 0>  IF  I  LOBE .QT.O)  CALL  TOLOBE (2 .AMIN1 (RCSC . SIG2) . -AwCPHP , -BwSPHP .0 . ) 

2020>  1  RCSC=AMIN1 ( RCSC , S I G1 *S IG2 ) 

2030>  RETURN 

2040>  END 

20S0>  C 
2060>  C 

2065>  C  <«><«> <«><«> <«><)(> <«><«><«><*><*><»> <«><»><*><»><»> <»><»> 

2070>  SUBROUTINE  WEDGE ( COOR , CTH , STH , CPH , SPH ) 

2072>  C 

2074>  C  WITH  LOBE 
2080>  C 

2090>  IMPLICIT  REALM  (A-Z) 

21 00>  INTEGERS  LOBE 

21 1 0>  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

21 20>  9  C2PHP.S2PHP.RAD, POLAR, LOBE, 

2 1 30>  S  ANGLE. L.B 

2 1 40>  REAL  COOR 1 9 ) 

21 50>  C 

2 1 80>  WEDMAXsB/ IWLENttWLEN) 

2 1 7  0  >  ANG=ANGLE«RAD 

2 1 8  0  >  M=(PlMLWL/(4 . ftANQWANQ) )MS I N IP  I MP I / (2 . NANG) l««2 

21 90>  TsCOS IP l«P I / (2 . NANG) ) 

2200>  CsCOS IP lnPHP/ANG) *T 

221 0>  Dsl.-T 

2220>  IFIABSICJ .LE. CUTOFF. OR. ABS(D) .LE. CUTOFF)  RETURN 

2230>  C 

2240>  I F I  POLAR . EQ . 2 . )  GO  TO  2 

2250>  RCSC  =  2 .» I -COOR 1 3 ) »SPH  +  COOR 1 8 ) *CPH ) *R2 -1 . 

2260>  RCSC=Ml*(flCSC/C  +  1  ./DINW2 

2270»  GO  TO  4 

2280>  C 

2290>  2  RCSC  =  2 . M (COOR ( 3 ) wCTH*CPH+COOR 1 6 ) WCTHWSPH-COOR 1 9 Iff 

2300>  t  STHI«Ht2-1  . 
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251 0>  RCSC=MM (RCSC/C+1 ,/D)mh2 

2320>  4  DU=ABS (P 1 2-ACOSK I -COOR 13  JMSTHMCPH -COOR 1 6 ) hSTHh 

2330>  I  SPH-COOR 1 9 ) MCTH) I 

2340>  RCSCsAMINI IRCSC.WEDMAX) 

23S0>  IF IDUIIDU.LT. CUTOFF)  GO  TO  6 

2360>  $ IQMAsRCSCH (WLEN/ IP IMLMDU) )MM2/8 . 

2370>  RCSCsAMINI (RCSC, SIGMA) 

2372>  I F ( LOBE . LE . 0)  RETURN 

2374>  CALL  TOLOBE 1 1 , RCSC/2 . , 0 . , 0 . , 0 . ) 

2376>  CALL  TOLOBE (2 .RCSC/2 ., 0 ., 0 ., L) 

2380>  RETURN 

c 

23B4>  6  IF  I  LOBE . LE . 0)  RETURN 

238B>  CALL  TOLOBE II , RCSC , 0 ., 0 ., L/2 . ) 

2388>  RETURN 

2390>  END 

2400>  C 
241 0>  C 

241 5>  C  <*><»><»><»><«><!•><»><«><«><*><*><*>  <«><*> <*><)(> <*><*><*> 

2420>  SUBROUTINE  WIRE ICOOR ,CTH . STH ,CPH , SPH) 

2430>  C 

2440>  C  WITH  LOBE 
24S0>  C 

2460>  IMPLICIT  REALM4  (A-Z) 

2470>  INTEGER*^  LOBE 

2480>  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

2490>  S  C2PHP.S2PHP, RAD .POLAR, LOBE, 

2500»  8  L  ,  A 

2 5 1 0 >  REAL  COORI9) . P (3) , ZK ( 3 ) 

2S20>  C 

2530>  IFISTHP.LE.O. )  RETURN 

2540>  GAMMA =1.78 

25S0>  ZK (1) =COOR 1 8  )M  I -CTH) -COOR 19 ) M I -STHMSPH) 

2580>  ZK 12) =- ICOOR I7)M|  -CTH) -COOR I9)M I -STHhCPH) ) 

2570>  ZK 1 3 ) =COOR 1 7 ) Ml -STHmSPH ) -COOR 1 8 ) M I -STHHCPH) 

2580>  IF (POLAR. EQ. 2.)  GO  TO  2 

2590>  P I 1 ) s-SPH 

2«00>  P(2) *CPH 

261 0>  PI3)=0. 

2620>  GO  TO  3 

2630>  2  PM)=CTHHCPH 

2640>  P (2) sCTHMSPH 

2650>  P 13 1 s-STH 

2@0O>  C 

2670>  3  SPSsZK  1 1 ) HP  1 1 )  +ZK  ( 2 ) HP  ( 2 )  «-ZK  ( 3  )HP  1 3 ) 

2680>  SPSbSPS/SQRT I ZK (1) HZK ( 1 ) +ZK 1 2 ) MZK 1 2 ) +ZK ( 3 ) HZK 1 3 ) ) 

2690>  PS=ASINIAMAX1 1-1 . , AM  I N 1 ( 1 . . SPS) ) ) 

2700>  RCSC=P IMLHLHCOS IPS)HN4/ IP I2MPI2+ALOQ IWLEN/ (GAMMAHP IMA) )MM2) 

2710>  IFISTHP.LT. CUTOFF. AND. LOBE. GT.O)  CALL  TOLOBE I 1 ,RCSC ,0 ., 0 . ,0 . ) 

2720>  IFISTHP.LT. CUTOFF)  RETURN 

2730>  SIGMA: ( WLENHTAN ( THP ) MCOS (PS) HM2 ) M»2 

27  40>  S IGMA  =  S I GMA/ 1 1 6 .HP IM (P 1 2MP 1 2*ALOG IWLEN/ (GAMMAHP IMAHSTHP) )HR2)) 

2730>  RCSCsAMINI (RCSC, 2. HSIGMA) 

2780>  I F I  LOBE , LE . 0)  RETURN 

2770>  CALL  TOLOBE 1 1 .RCSC/2 ., 0 . ,0 ., L/2  . ) 

2780>  CALL  TOLOBE 12 , RCSC/2 ., 0 ., 0 ., -L/2  . ) 

2790>  RETURN 

2800>  END 

281 0>  C 
2820>  C 
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282 S>  C  <«><#> <#><*><»><*><»><»><»><*><»><*> <«><»><»><*> <#><!»><»> 

2830>  SUBROUTINE  LWI RE (COOR ,CTH , STH ,CPH , SPHI 

2840>  C 

2850>  C  WITH  LOBE 
2860>  C 

2870>  IMPLICIT  REAL84  IA-Z) 

2880>  INTEGERa2  LOBE 

2880>  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP . CTHP , STHP , PHP , CPHP . SPHP , 

2900>  8  C2PHP , S2PHP , RAD , POLAR , LOBE , 

2  9 1  0  >  8  A 

2920>  REAL  COOR (9) 

2930>  C 

2940>  I F (POLAR . EQ . 2  . )  GO  TO  2 

2950>  GAMMA =ACOSK < -COOR 1 2 1 aSPH+COOR ( S ) NCPHI 

2960>  GO  TO  3 

2970>  C 

2980>  2  GAIMA : ACOSK ( COOR ( 2 1 aCTHRCPH+COOR ( 5 1 RCTHaSPH-COOR  1 8 ) ftSTH ) 

2990>  C 

3000>  3  KA2= 1 4 . aP I aA/WLEN J  aSTHP 

301 0>  CG2=COS( GAMMA )aa2 

3020>  SG2C=S IN (GAMMA Jaa2aCTHPaCTHP 

3030>  CALL  BESJ (KA2,0, JO, . 1 , IER) 

3040>  IF(IER.NE.O)  GO  TO  9 

3050>  CALL  BESJ(KA2,2,J2, . 1 , IER) 

3060>  IF  I IER.NE . 0)  GO  TO  9 

3070>  RCSC = (SG2C+CG2)aJ0+(SG2C-CG2)aJ2 

3080>  RCSC=P laAaAaRCSCaRCSC 

3090>  I F I  LOBE . LE . 0)  RETURN 

3100>  IF (STHP. LE. CUTOFF)  CALL  TOLOBE ( 1 , RCSC , 0 . , 0 . , 0 . ) 

31 1 0>  IF (STHP. QT. CUTOFF)  CALL  TOLOBE ( 1 .RCSC .AaCPHP .AaSPHP , 0 . ) 

3 1 2  0  >  RETURN 

3130>  9  RCSsO. 

31 40»  RETURN 

31 50>  END 

31 60>  C 
31 70>  C 

3 1 75 >  C  <*><a><a><*><a><a><a><a><a><a><a><a><a><a><a><a><axa><a> 

31 80>  SUBROUTINE  CAVA 

31 90>  C 

3200>  C  WITH  LOBE 
321 0>  C 

3220>  IMPLICIT  REALa4  (A-ZI 

3230>  I  NTEGERa2  LOBE 

3240>  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 12 , RCSC , THP , CTHP . STHP , PHP , CPHP . SPHP , 

3230>  8  C2PHP.S2PHP, RAD, POLAR, LOBE, 

3260>  8  A 

3270>  C 

3280>  IF (ABS(THP) . GT.70.aRAD)  RETURN 

3290>  T= 14 .aP laA/ IWLENaWLEN) laal .5 

3300>  RCSC= . 4aTaWLENaWLEN 

331 0>  TTHP=TAN(THP) 

3320>  IFITTHP.LT. CUTOFF. AND. LOBE. GT. 01  CALL  TOLOBEI1 .RCSC.O. .0. ,0. ) 

3330>  IFITTHP.LT. CUTOFF)  RETURN 

3340>  SIGMA: .8aP laA/TTHP 

3350>  RC5C3AMIN1 (RCSC , SIGMA) 

3360>  I F I  LOBE . LE . 0 )  RETURN 

3370»  DsSQRT(4.aA/PI ) 

3380>  DS23D8STHP/2. 

3390»  CALL  TOLOBE ( 1 , RCSC , DS2aSTHPaCPHP , DS2aSTHPaSPHP , DS2aCTHP ) 

3400>  RETURN 
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3410>  END 

3420>  C 
3430>  C 

3435>  C  <9><9><9><9><9><9><9><9><9><9><9><9><N><9> <9><9><9> <9><9> 

3440>  SUBROUTINE  CAVD 

3450>  C 

3460>  C  WITH  LOBE 
3470>  C 

34S0>  IMPLICIT  REAL94  (A-Z) 

3490>  INTEGER92  LOBE 

3S00>  COMMON/RCSCOM/CUTOFF ,WLEN, P I , P 12 ,RCSC , THP , CTHP ,STHP ,PHP ,CPHP , SPHP , 

3510>  9  C2PHP.S2PHP, RAD, POLAR, LOBE, 

3520>  9  D 

3530>  C 

3S40>  IFIABS(THP) .GT.70.9RAD)  RETURN 

3$S0>  RCSC= .  059 (2 .9P  I 9D/WLEN19939WLEN9WLEN 

3560>  TTHPsTANITHP) 

3570>  IFITTHP.LT. CUTOFF. AND. LOBE. QT.O)  CALL  TOLOBEI1 ,RCSC,0. ,0. ,0. ) 

35B0>  IFITTHP.LT. CUTOFF)  RETURN 

3590>  SIGMA=  . 059 ( 2 .9P  I9D) 992/TTHP 

3B00>  RCSCsAMINI (RCSC.SIQMA) 

361 0>  IF  I  LOBE . LE . 0)  RETURN 

3630>  DS2=D9$THP/2. 

3640>  CALL  TOLOBE ( 1 , RCSC , DS29STHP9CPHP , DS29STHP9SPHP , DS29CTHP ) 

3650>  RETURN 

3660>  END 

3670>  C 
3690>  C 

3635 >  C  <9><9><9><9><9><9><9><9><9><9><9><9><9><9><9> <»><«> <9><9> 

3690>  SUBROUTINE  CFLAT 

3700>  C 

3710>  C  WITH  LOBE 
3720>  C 

3730>  IMPLICIT  REAL94  (A-Z) 

3740>  INTEQERK2  LOBE 

3750>  COMMON/RCSCOM/CUTOFF , WLEN . P I . P 1 2 , RCSC , THP . CTHP . STHP , PHP , CPHP , SPHP , 

3760>  9  C2PHP.S2PHP, RAD .POLAR, LOBE, 

3770>  9  A 

3790>  C 

3790>  IFIABSITHP) .GE.PI2)  RETURN 

3800>  RCSCM.«Pltt(PI«A«A)»tt2/IWLEN«WLEN) 

3310>  D*8.*P iMSTHPttTANITHP )ttK2 

3620>  IFID.LT. CUTOFF. AND. LOBE. QT.O)  CALL  TOLOBE ( 1 , RCSC ,0. ,0. ,0. I 

3830>  IFID.LT. CUTOFF)  RETURN 

3840>  SIQMA:A«WLEN/D 

3850>  RCSCSAMINI IRCSC, 2. NSIQMA) 

3860>  IF (LOBE . LE . 0)  RETURN 

3870>  CALL  TOLOBE 1 1 , RCSC/2 . , AttCPHP , AtfSPHP , 0. ) 

3880>  CALL  TOLOBE (2 , RCSC/2 ., -ANCPHP , -AUSPHP , 0 . ) 

3890>  RETURN 

3900>  END 

391 0>  C 
3920>  C 

3925>  C  <»><»> <»><9><»><»><ll> <»><»> <9 >««>«»><»><»><»> <» ><9> <«><«> 

3930>  SUBROUTINE  RFLATI (PASS) 

3940>  C 

3950>  C  WITH  LOBE 
3960>  C 

3970>  IMPLICIT  REAL94  (A-Z) 

3980>  INTEQER92  LOBE.IPASS 
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3990> 

4000> 

401 0> 

4 020>  C 

4030> 

4040> 

4050> 

4060> 

4070> 

4080> 

4090> 

41 00> 

4 1 1 0  > 

41 20>  1 
41 30> 

41 40>  C 
41  50> 

41 60> 

41  70» 

41 80> 

4 1  9  0  > 
4200> 

4  2 1  0  >  C 
4220>  C 
4230>  2 
4240> 
4250> 
<*260> 
4270> 
4280> 
4290>  5 
4300> 
431 0>  C 
4320>  C 
4330>  3 
4340» 
4350> 
4360> 
4370> 
4380» 
4390> 
4400>  C 
441 0>  C 
4420>  4 
4430> 
4440> 
4450> 
4460> 
4470>  C 
4480»  C 
4485 >  C 
4490> 
4492*  C 
4494>  C 
4500>  C 
4510> 
4520> 
4530> 
4540> 
4550> 


COMMON /RCSCOM/CUTOFF .WLEN.PI , P i 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 
S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

9  A ,  B 

IF (ABS(THP) . GE ,P 12)  RETURN 
RCSC=PI»I8.*A*B/WLEN)*#2 
S2THP*STHP*STHP 
S4THP=S2THP#S2THP 
IFIS4THP.LT. CUTOFF!  60  TO  4 
IF(S2PHP.LT. CUTOFF)  GO  TO  2 
IF(C2PHP.LT. CUTOFF)  60  TO  3 

S I QMAsWLENNWLEN/ (64 .«P l**3«S4THP*S2PHP*C2PHP ) 

SIGMAs2.ttSI6MA 
RCSC=AMIN1 I RCSC, 2. NS  1 6MA) 

I F ( LOBE . LE . 0 )  RETURN 

/  a 

I F ( (PASS . NE .2)  CALL  TOLOBE ( 1 , RCSC4 , A , B, 0 . ) 

I F 1 1  PASS . NE . 2)  CALL  TOLOBE 12 , RCSC4 , -A , B. 0 . ) 

IF) JPASS.EQ.2)  CALL  TOLOBE ( 1  , RCSC4 , A , -B , 0 . ) 

I F ( I  PASS . EQ . 2 )  CALL  TOLOBE 12 , RCSC4 , -A , -B ,0 . ) 

RETURN 

PHI’  :0  OR  180 

IF  I IPASS.EQ.2)  RETURN 
SIGMA:BMB/(4.«PI»S2THP) 

I F I  LOBE . LE . 0)  GO  TO  1 
S I G  =  AM  I N1 (RCSC/2. , SIGMA) 

CALL  TOLOBE (1 , S IQ, A , 0. ,0. ) 

CALL  TOLOBE 1 2  ,  S I Q ,  -  A  ,  0 . ,0.) 

RCSCsAMINI I RCSC, 2 .ttSIQMA) 

RETURN 

PHI’  s  90  OR  -90 

IF ( IPASS.EQ.2)  RETURN 
SIGMAsA»A/(4.*PI*S2THP) 

I F I  LOBE . LE . 0 )  60  TO  1 
S I G  =  AM  1 N 1 (RCSC/2. .SIGMA) 

CALL  TOLOBE (1,SIG,0.,B,0.) 

CALL  TOLOBE (2,SIG,0. , -8,0.1 
GO  TO  3 

THETA’  s  o 

I F ( IPASS.EQ.2)  RETURN 

I F ( LOBE . LE . 0 )  RETURN 

CALL  TOLOBE ( 1 ,RCSC,0. ,0. ,0. ) 

RETURN 

END 


<«><«>  <«><«><<(><*>  <*><*>  <»><#>«»><*><»><*><*>  «*><*>  <»><*> 

SUBROUTINE  TVAVE 

WITH  LOBE 

IMPLICIT  REALM  (A-Z) 

INTEGERM2  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP ,CPHP , SPHP , 
S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

8  P , L , GAMMA , A 
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4S60 

4570 

4550 

4590 

4600 

4610 

4620 

4630 

4640 

4645 

4650 

4660 

4670 

4680 

4685 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4775 

4780 

4790 

4795 

4800 

4810 

4820 

4830 

4840 

4841 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4935 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 


C 


C 

C 

C 

C 


C 


C 


1 

2 

4 


3 


C 

C 

c 

c 

c 

c 


IF (POLAR. EQ. 1 . )  RETURN 
KLS2.MPINL/WLEN 
I F  (KL  .  LE  .  4 . 1  RETURN 
Qs2.35*AL0G10(2.*KL)-.55 
T*1 .-P*CTHP 
IFIT.LE. CUTOFF)  RETURN 
FtSTHPMSlN(KLftT/l2.«Pn/T 
RCSC: (GAMMANFttF*WLEN/Q) **2/P I 

IF  I  LOBE . QT . 0)  CALL  TOLOBE 1 1 ,RCSC , A*CPHP , AKSPHP , -L/2 . ) 
RETURN 
END 


<*><«><«>«*><«><*><*> <*><«> <*><»><*> <«><»><«> <»><»> <«><*> 

SUBROUTINE  DINED 

IMPLICIT  REALM  (A-ZI 
INTEGERS  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP . STHP , PHP , CP  HP . SPHP , 
9  C2PHP.S2PHP, RAD, POLAR, LOBE, 

8  A.B.C 

IFIPHP.LE.O. .OR.PHP.GE.PI2I  RETURN 
PI4=PI/4. 

GAMMA=ATAN2 1 B , A I 
DU*PI2-THP 

COMPUTE  NORMAL  COMPONENT 
IFIPHP.GT. GAMMA)  GO  TO  1 
RCSCM .  HP  I  *  I ANB/VLEN  J  M2M .  NS2PHP 
GO  TO  2 

RCSCM .  HP  I N 1  C«B/WLEN  I  ««2M .  HC2PHP 
I F (DU**2 . GT .CUTOFF!  GO  TO  4 
I F ( LOBE . LE . 0)  RETURN 
IFISTHP.LT. 0.51  GO  TO  3 
SIQMA=WLEN»WLEN«RCSC/(8.M(Pi«BNDU)N»2) 

RCSC=AMIN1 I RCSC. SIGMA  I 

RETURN 

RCSC=0. 

RETURN 

END 


<«><«>  <«><»><»> 

SUBROUTINE  TCORN 
WITH  LOBE 

IMPLICIT  REALM  IA-ZI 
INTEGERS  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 
9  C2PHP.S2PHP, RAD, POLAR, LOBE, 

8  A.B.C 

I F I THP . LE . 0 . 0 . OR . THP . GE . P 1 2 1  RETURN 
IFIPHP.LE.O.O.OR.PHP.GE.PI2)  RETURN 
CALL  ORDER  I STHPNCPHP , STHPKSPHP , CTHP , L ,M, N) 

CALL  ORDER  I A , B ,C , T1 , T2 , T3  J 

LASL/T3 

MB&M/T2 
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5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5180 

5200 

5210 

5220 

5230 

5235 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 


5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 

5520 

5525 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

5610 

5620 

5630 

5640 

5650 

5660 

5670 


1 

2 


C 

C 

C 

C 

C 

C 


C 


1 

2 


C 

C 

C 

C 

C 

C 


C 


NC=N/T1 

I F I LA+MB . GT . NC )  Q0  TO  1 
A  RE  AM  .  nL*M*T1  /  I  LA+MB+NC ) 

GO  TO  2 

AREA: ( LAttLA+MBMMB+NCttNC) / ( LA+MB+NC) 

AREAMttBttCttl  LA+MB  +  NC-2  ,#AREA) 

RCSCM  .  #P  I  MARE  AH  ARE  A/  (WLENttWLEN) 

IF(IOBE.LE.O)  RETURN 

D*1 . /SORT ( 1 . /A«M»2*1 . /BNN2+1 ,/C*»2) 

CALL  T0L0BEI1 . RCSC , D#STHP*CPHP , D#STHP*SPHP , DaCTHP J 
RETURN 
END 


<MXKXttXftX«XNXtfXttXaX«XMXtfXKXMX*XNXIIXttX»> 

SUBROUTINE  RCORN 
WITH  LOBE 

IMPLICIT  REALM  (A*Z> 

I NTEQER*2  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 
8  C2PHP.S2PHP, RAD, POLAR, LOBE, 

*  A.B.C 

I F I THP . LE . 0 . 0 . OR . THP . QE . P 1 2 )  RETURN 
I F IPHP . LE . 0 . 0 . OR . PHP . GE . P 1 2 )  RETURN 
CALL  ORDER ( STHPMCPHP . STHPttSPHP , CTHP , L ,M , Nl 
CALL  ORDER(A,B,C,T1 ,T2,T3) 

LA=L/T3 

MBSM/T2 

NC=N/T1 

I F (MB . LT . NC/2 . )  GO  TO  1 
ARE A  s  L*T2# ( 4 . »T 1 -N/MB I 
GO  TO  2 

AREAM.*L»M»T1»T1/N 
RCSC =4 . *P  lit  ( AREA/WLEN)  ttM2 
IF (LOBE . LE .0)  RETURN 
D=1 . /SORT  11 ./A«»2* I . /C»»2) 

CALL  TOLOBEI1  , RCSC , D»STHP»CPHP , D*STHP»SPHP , D»CTHP I 
RETURN 
END 


<MX«XMX«Xil>  <*X«>  <MXKXKX<0<  PX  »X»XKXNX»XttX»> 

SUBROUTINE  ECORN 
WITH  LOBE 

IMPLICIT  REALM  (A-Z) 

INTEGERM2  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 
t  C2PHP.S2PHP, RAD, POLAR, LOBE, 

*  A.B.C 


IFITHP .LE.O.O.OR.THP.GE.PI2)  RETURN 
I F (PHP . LE . 0 . 0 . OR . PHP . GE . P 1 2 1  RETURN 
CALL  ORDER  I STHPhCPHP , STHP9SPHP ,CTHP , L ,M,N) 
CALL  ORDER (A , B , C , T 1 ,T2,T3I 
LA2=IL/T3HH»2 
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5680> 

5690> 

5700> 

57t0> 

5720> 

5730> 

5740> 

5750> 

S760> 

5770>  1 

5780> 

57#0>  2 

5800> 

581 0> 

5820> 

5830> 

S840> 

5850> 


MB2* IM/T2 ) #M2 
KC2=(M/T1 |N«2 
I F (NC2 . GE . LA2+MB2 I  QO  TO  1 

ABEA= I N/T1 IMATAN I ( I LA2+MB2 ) MM2-NC2MNC2 ) / (4.MNC2M I L/T3 )M (M/T2) ) ) 
ARE As ARE A+ (M/T2 INATAN I ( I LA2+NC2 ) MM2-MB2NMB2 ) / 

3  (4.MMB2MIL/T3IMIN/T1 ) ) ) 

ARE As ARE At ( L/T3 ) MATAN I  1 (MB2  +  NC2 IMM2-LA2MLA2 I / 

I  (4.MLA2MIM/T2IMIN/T1 ) ) ) 

QO  TO  2 

AREAS (M/T2)ftATAN(2.K(L/T3)tt(N/Tl) / (MB2+NC2-LA2 ) ) 

AREA s AREA + IL/T3) MATAN 12 .*IM/T2IM(N/T1 )  / 1 LA2-MB2+NC2 1 1 
AREAsAMBMCMAREA 
RCSC=4 .mP |MIAREA/WLEN)mm2 
IF  I  LOBE . LE . 0)  RETURN 
D=1 ./SQRTI1 ./AMM2+1 ./BMM2+1 ./CM*2) 

CALL  TOLOBE 1 1 , RCSC , OMSTHPMCPHP , DMSTHPMSPHP , DmCTHP ) 

RETURN 

END 
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FILE  a  SYO : [ 1 1 1 , 1 2 ] RCSSUB . FTN 


27-SEP-83  14:35:4®  ««« 


10 

20 

30 

40 

50 

51 

52 

53 

54 

55 
60 
70 
75 
60 
90 

100 

110 

120 

130 

170 

190 

190 

220 

222 

224 

226 

230 

240 

250 

260 

270 

260 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

450 

460 

470 

480 

482 

484 

490 

500 

510 

520 

530 

540 

550 


THIS  ROUTINE  ASSUMES  MATRICES  ARE  ROW  STRUCTURED 
SUBROUTINE  TOLOBE 1 1C , RCSL , X , Y , Z) 

INTEGERM2  DEBUG .COMB  I N ,GL I  NT 

LOG ICALN2  GOULD , LP .OPENE , LPF .OLDGRF .MGRAPH .GDAUTO 

COMMON /SWITCH/  IOW, THETA  14) .GOULD , LP , OPENE , LPF , OLDGRF .MGRAPH , 

1  GDAUTO, DEBUG, COMB  IN, GLINT, I  BUG , LBUQ , ISTATS, IZETA 

REALM  SLOBE  14,181  ,2) 

REALM  0RGN13) , COMP ( 9 ) .MTRAN (3 ,31, TMP (3 ) 

COMMON /RCSL OB/ I OL , ORGN , COMP , SLOBE , BETA , 

9  THR  , CTHR , STHR , PHR , CPHR , SPHR , KLOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP . CTHP , STHP , PHP , CPHP , SPHP , 
8  C2PHP.S2PHP, RAD, POLAR .LOBE, 

S  ARR 


C 

C 

C7 

C 


REALM  RCSI181 1  .RCS21181 ) 

COMMON  /RPRCS/  I  RPR , I RCS , RAM , RCS . RCS2 

IF (LBUG.NE .01  WRITEI2.7)  I C  ,  RCSL , X , Y , Z 
FORMAT C  TOLOBE= ' , 11 , 4E1 1 . 3 ) 

TMP  ( 1.1  -X 
TMP 12) =Y 
TMP ( 3 )  *Z 

CALL  MATVEC I  COMP , TMP , TMP ) 

TMP  Mis  TMP  m+ORGNM) 

TMP ( 2 ) = TMP ( 2 ) +ORGN 1 2 ) 

TMP ( 3 ) s TMP ( 3 ) ♦ORGN 1 3 ) 

CALL  MATPUT (MTRAN . CTHRNCPHR ,CTHR»SPHR , -STHR , 

1  SPHR, -CPHR. 0. . 

2  -STHRaCPHR , -STHRUSPHR , -CTHR) 

I F (BETA . LE . 5 . )  GO  TO  1 

CALL  MATVEC (MTRAN. TMP, TMP) 

IFI  (LBUG.AND.2) .EQ.OI  GO  TO  5 

WR I TE ( 2 , 4 )  I C , THR , PHR , X , Y , Z , ORGh .MTRAN 

FORMAT (IX, I2.8E10.3/3E11 , 3/3E1 1  . 3/3E1 1  . 3 ) 

I F ( IRPR.EO.O)  GO  TO  3 
ARQ:4 . »P I  * TMP ( 3 ) /WLEN 
SQsSQRT (RCSLftRAM) 

RCS ( I RCS) 8RCS ( I RCS ) ♦ SQmCOS ( ARG ) 

RCS2 ( I RCS) =RCS2 ( I RCS I *SQ«S IN (ARGI 

I F 1 1 LBUG . AND . 1 ) . HE . 0 ) WR I TE ( 2 , 6 ) I RCS , RCSL . RAM . TMP . RCS 1 1 RCS ) , RCS2 1 1 RCS ) 
FORMAT!’  RCSL* ’ , I3.7E11 .3) 

IF (KLOBE . LE . 0)  RETURN 

SLOBEd  ,  I  RCS ,  I C )  8RCSL0RAM 
SLOBE ( 2 , I RCS , 1C i = TMP 1 1 ) 

SLOBE ( 3 , I RCS , IC)sTMP(2) 

SLOBE (4 , I RCS , I C ) =  TMP ( 3 ) 
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560> 
570>  2 
580>  C 
590> 
600> 

81 0>  C 
620>  C 
630>  C 
640> 
650> 
660> 
670> 
680> 
690> 
700> 

71 0> 
720>  C 
730>  C 
7 40>  C 
750>  C 
760> 
770> 
780> 
790> 
800> 
810> 
820> 
830> 
840>  3 
8S0>  2 
860>  1 
870> 
880>  4 
890> 
900> 

91 0>  C 
920>  C 
930>  C 
940>  C 
950> 
960> 
970> 
980> 
990> 
1000> 

1 01 0>  3 
1 020>  1 
1 030> 

1 040>  4 
1 050> 

1 060» 

1 070>  C 
1 080>  C 
1 090> 

1 1 00> 

1 1 1 0> 

11 20> 

1 1 3  0  > 

1 1  40> 

1 1  50> 


IF ( ( LBUG. AND .  1 ) . NE . 0)  WRITE  1 2, 2)  1C,  IRCS,  ISLOBEIJ ,  IRCS,  10,1=1,4) 
FORMAT ( ’  TOLOBE>  ' , 2  1 2 , El  1 . 3 , 3X , 3F8 . 3 I 

RETURN 

END 


SUBROUTINE  ORDER (A1 , A2 , A3 , B1 , B2 , B3 I 

TIsAMINI ( A 1 , A2 , A3  I 

T3=AMAX1 C  A 1 ,A2 . A3 1 

B2=A1+A2+A3-T1-T3 

BIsTI 

B3  =  T3 

RETURN 

END 


THIS  ROUTINE  ASSUMES  MATRICES  ARE  ROW  STRUCTURED 

SUBROUTINE  MATMUL I A , B , C) 

REALM  A 1 9 ) ,B(91 ,T(9) , C ( 9 ) 

DO  1  J  =  1 ,3 
J 1 3= ( J-1 ) K3 
DO  2  K  =  1  ,3 
JK5J13+K 
T (K+ J 1 3 ) =0 . 

DO  3  L  =  1  ,3 

TIK*J13)sTIK*J13) »A  I J 1 3  *L  )  *B  IK  ♦  I  L  - 1 1  »3 1 
CONI .NUE 
CONTINUE 
DO  4  J=1,9 
C ( J I =T I J ) 

RETURN 

END 


THIS  ROUTINE  ASSUMES  MATRICES  ARE  ROW  STRUCTURED 

SUBROUTINE  MATVEC ( A , B , C ) 

REALM  A 1 9 )  ,  B 1 3 1  ,  C 1 3 )  ,T(3) 

DO  1  J  =  1  ,3 
J13= ( J-1 1*3 
TIJisO. 

DO  3  L=1 ,3 

T(J)=T(J)+A(J13+L)#B(L) 

CONTINUE 
DO  4  J=1 ,3 
CU)sT(J) 

RETURN 

END 


SUBROUTINE  MATPUT ( A , A 1 1 , A1 2 , A 1 3 , A2 1 , A22 , A23 , A3 1 , A32 , A33J 

REALM  A  f  9 1 

ACIIsAII 

A  12 IsA  12 

AI3ISA13 

A  1 4 ) =A21 

A(5)sA22 
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1 1 60> 

A(6I=A23 

1 1 70  > 

A(7 1 =A31 

1 1 80> 

A  ( 8 )  =  A  3  2 

1 1 90> 

A (9 ) =A33 

1 200> 

RETURN 

1 2 1 0  > 

END 

t  220> 

C 

1 2  3  0  > 

C 

1  2  40  > 

SUBROUTINE  TRANSPIC.E) 

1 2  5  0  > 

REAL  C ( 9 ) ,D(9) . E (9 I 

1  2  6  0  > 

D(1|:CI1I 

1270> 

D(2)=CI4) 

1  2  8  0  > 

D13I=C(7) 

1 290> 

0(4) sC (21 

1 3  0  0  > 

D  (  5 ) =C 1 5  ) 

1 31 0> 

D(6)=C(8) 

1  320> 

D ( 7 ) sC ( 3 ) 

1  3  3  0  > 

D ( 8 ) =C ( 6 ) 

1  3  4  0  > 

0(9) sC(9) 

1  3  5  0  > 

DO  1  J  =  1  ,9 

1 360> 

1 

E(J)=DtJ) 

1370> 

RETURN 

1  3  8  0  > 

END 

1 390> 

C 

1 400> 

C 

1  4 1  0  > 

C 

1  420> 

SUBROUTINE  BESJ (X,N,BJ ,D, IER) 

1  430> 

C 

1 440> 

C 

X=AR6UMENT  OF  THE  J  BESSEL  FUNCTION  DESIRED 

1  450> 

C 

N=ORDER  OF  THE  J  BESSEL  FUNCTION 

1460> 

C 

BJ=RESULTANT  J  BESSEL  FUNCTION 

1470> 

C 

D=REQUIRED  ACCURACY 

1 480> 

C 

IER=RESULTANT  ERROR  CODE 

1  4  9  0  > 

C 

0=NO  ERROR 

1 500> 

C 

1 =N  IS  NEGATIVE 

151 0> 

C 

2=X  IS  NEGATIVE 

1  520> 

C 

3  =  REQU 1  RED  ACCURACY  NOT  OBTAINED 

1  5  3  0  > 

C 

ARRANGE  OF  N  COMPARED  TO  X  NOT  CORRECT 

1  540> 

C 

1  550> 

BJ=0. 

1 560> 

IF  IN) 10,20,20 

1  570> 

10 

I ER  =  1 

1  5  8  0  > 

RETURN 

1  590> 

20 

1 F 1 X ) 33 , 30 , 3 1 

1 600> 

30 

IF(N.GT.O)  RETURN 

1  6 1  0  > 

BJ si  . 

1  6  2  0  > 

RETURN 

1 630> 

33 

IER  =  2 

1 640> 

RETURN 

1  6  S  0  > 

31 

IFIX-15. 132,32,34 

1 660> 

32 

NTESTs20.+10.#X-XlM»2/3 

1 670> 

GO  TO  36 

1 680> 

34 

NTEST=90. ♦X/2. 

1  6  9  0  > 

36 

IFIN-NTE ST  140,38,38 

1  700> 

38 

IERs4 

1 7 1  0  > 

RETURN 

1720> 

40 

IERsO 

1 730> 

N1 sN+1 

1  7  40> 

BPREV=0. 

1 750> 

C 
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1 7  6  0> 

C 

COMPUTE  STARTING  VALUE  OF  M 

1 770> 

C 

1780> 

IFIX-5. 150,60,60 

1 7»0> 

50 

MA=X*6. 

1 800> 

GO  TO  70 

1 81 0> 

60 

MAsi . 49X+60. /X 

1 8  2  0  > 

70 

MB=N*IFIXIX)/4+2 

1 8  3  0  > 

MZERO-MAXO (MA .MB) 

1  840> 

C 

1 8  5  0  > 

C 

SET  UPPER  LIMIT  OF  M 

1 860> 

C 

1 8  7  0  > 

MMAX=NTEST 

1 880> 

100 

DO  190  M:MZERO,MMAX,3 

1 8#0> 

c 

1 900> 

c 

SET  F (Ml , F (M-1 ) 

1 91 0> 

c 

1  920> 

FM1*1 .1-15 

1  930> 

FM=0. 

1  9  40  > 

ALPHAsO. 

1 9  5  0  > 

IF (M- IM/2IM2) 120, 1 10 , 120 

1 960> 

110 

J  T= - 1 

1 970> 

GO  TO  130 

1 9  8  0  > 

120 

JTsI 

1 990> 

130 

M2=M-2 

2000> 

DO  160  K=1 , M2 

201  0> 

MK=M-K 

2020> 

BMK  =  2  .  **F  LOAT IMK I  «FM1  /X-FM 

2030> 

FM=FM1 

2040> 

FM1 =BMK 

2050> 

IF(MK-N-1M50,140,150 

2060> 

140 

BJSBMK 

2070> 

150 

JTs-JT 

2080> 

S=1 ♦ JT 

2090> 

160 

ALPHA=ALPHA+BMK*S 

21 00> 

BMK=2 . #FM1 /X-FM 

2 1 1 0  > 

IF (N) 180.170.180 

21 20> 

170 

BJ=BMK 

2 1  3  0  > 

180 

ALPHA=ALPHA+BMK 

21 40> 

BJsBJ/ALPHA 

21  50> 

1 F  < ABS  (BJ-BPREVI -ABS ID*BJ 1 1200,200 ,190 

21 60> 

190 

BPREVSBJ 

21  70> 

1 ER-3 

21  80> 

200 

RETURN 

21 90> 

END 

2200> 

C 

2  2 1  0  > 

C 

2220> 

SUBROUTINE  TRIM80(NB,BI 

2230> 

LOG  1  CALttl  B  f 80) , BL 

2240> 

DATA  Bl/’  ’/ 

2250> 

NB  =  80 

2260> 

DO  1  J  =  1  ,80 

2270> 

IFIB(NB) .NE.BLI  RETURN 

2280> 

1 

NB*NB-1 

2290> 

RETURN 

2300> 

END 

2  3 1  0  > 

c 

2320> 

c 

2330> 

FUNCTION  ACOSKIAI 

2340> 

ACOSK  =  ACOS (AMAX1 1-1 . , AMIN1 C 1 . , A 1 1 ) 

23S0> 

RETURN 

2380> 

END 
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FILE  b  SYO : [ 1 1 1 ,1 2 1H0SBD . FTN 
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10 


20>  C 

H  H 

000 

SSSS 

B6BB 

DDDD 

30>  C 

H  H 

0  0 

$ 

B  B 

D  D 

40>  C 

HHHHH 

0  0 

SSS 

BBBB 

D  D 

50>  C 

H  H 

0  0 

s 

B  B 

D  D 

60>  C 

H  H 

000 

SSSS 

BBBB 

DDDD 

70>  C***» 
80>  C 

iNMHXttXMl 

90> 

SUBROUTINE  HOSBO (1  OUT, LOOP , IOE . 

1  OP, SLOBE, 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

322 

324 

330 

340 

350 

360 

370 

380 

390 

400 

410 

412 

414 

416 

418 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 


CttKN**************************************************** 


1  FREQ, POL , PERSH , P ITCH , ROLL , YAW, T I LT , PHR I 

REAL«8  T1 ,T2,T3,T4,S1 ,$2,$3,S4,T12, S22 , VI , U2 , U3 , U4 ,G1 ,Q2 
REALMS  B , N , P , RCSMAX , RCSM I N , T5 , T6 
REALX4  SLOBE (8,181) , SAVE (30) 

BYTE  BUF (1 1 

INTEGERM2  DEBUG , COMB  I N , QL I  NT 

LOG ICAL*2  GOULD , LP .OPENE , LPF . OLDGRF , MGRAPH , GDAUTO 

COMMON /SWITCH/  I OW , THETA ( 4 ) , GOULD , LP , OPENE , LPF , OLDGRF .MGRAPH , 

1  GDAUTO, DEBUG, COMB  IN, GLINT, I  BUG, LBUG, I  STATS, IZETA,STATS,NRCS, IOH 

LPsLOOP 

IF  I IOUT.EQ.O)  GO  TO  20 

COMPUTE  REQUIRED  ANGLES  FOR  GRAPH 

SL0PE  =  FL0ATINRCS-1 1 /AMAX1 1 .001 , THETA (2) -THETA M 1 1 
C0NST=1 .-SLOPEmTHETAII) 

ANGsSTATS 

NSsO 

I F  ( ANG .  GE .  THETA  Ml)  GO  TO  12 

ANG=ANG+STATS 

GO  TO  11 

ANG=ANG+STATS 

NS=NS-1 

IF ( ANG. GT. THETA (2) I  GO  TO  13 
LP= I F IX (SLOPEXANG+CONST+ . 001 ) 

NS=NS+1 


C 

C 

C 


11 


12 


C 

20 

C 


READUOE’LP)  T1  , T2  , T3  , T4 , T5  , T6 

SI  =T1 

T12=TH»T1 

S2= (T1 2-T2 ) /2 . DO 

IF IS1*S2.EQ.O. DO. AND. IOUT.EQ.O) 
IF (S1*S2 .EQ.O.DO)  GO  TO  4 


RETURN 


S3s(T12*T1-3.D0*T1*T2+2. D0*T3 I /6 .DO 

S4s(T12*T12-6. D0*T12*T2+3 . D0*T2*T2+8 . D0*T1*T3-6 . D0*T4J /24.D0 
VI  =S1 

U2=2 . D0MS2 
U3=1 2 . D0*53 
S22BS2XS2 

U4b6. DO* S22+12. DO* 51*53+132. D0»S4 
G1=8.D0»S3/(S1*S2) 

G2=-1 .500+3 .DO*S1*S3/S22+33 .D0MS4/S22 
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540> 
550> 
560> 
570> 
S80> 
590> 
600>  C 
61 0> 
620> 
630> 
640> 
650> 
680> 
670>  C 
680>  5 
690>  1 
700> 
710>  2 
720> 
730> 
740>  3 
750» 
760>  C 
770> 
780>  C 
790>  13 
800>  C 
81  0»  C 
820>  C 
825> 
830> 
832> 

8  3  5  > 
840» 
850> 
860> 
870>  C 
880> 
890> 
900> 
910> 
920> 
930>  C 
940> 
950> 
960> 
970>  C 
980» 
990> 
1000> 

1 01 0> 

1 020> 

1 030> 

1 040> 

1  050> 

1 060> 

1 070> 

1 080> 

1 090> 
1100> 


Bs2.D0«S2/1S1*S1) 

N=S1*(S1ttS2+3  .DON  S3)  /  (2  .D0*S22-3.D0MS1MS3) 

Psi .D0/!N#B) 

RCSMAXST5KN2 
RCSMINs|2.D0*T6-T5)»N2 
IF  t IOUT.EQ.O)  GO  TO  5 

SLOBEM  ,HS)sAHG 
SLOBE (2 , NS)  =B 
SL0BFI3,NS)=G1 
SLOBL (4,NS) =G2 
ANG=ANG*STATS 
GO  TO  12 

WRI TE  HOP  ,  1 )  SI  ,  S2  ,  S3  ,  S4 

FORMAT!’  SI*’, Dll. 3,’  $2=’, Dll. 3,’  S3*’ .01 1 .3,  ’  S4*’,D11.3> 
WRITE! I  OP , 2 1  VI ,U2,U3,U4,Q1 ,G2 

FORMAT!'  Vis’, Dll. 3,’  U2»’,D11.3,’  03s't011.3,’  U4a*,D11.3. 

1  ’  G1 = ' ,011.3,’  G2s ’ ,011.3) 

WRITE! IOP , 3 )  B ,  N ,  P , RCSMAX , RCSM I N 

FORMAT!’  W= ’ ,011 . 3 , ’  Ns’, Dll. 3,’  Ps • ,D1 1 . 3 , ’  RCSMAXs’ , D1 1 .3 , 
1  ’  RCSM I  Ns’ ,011.3/) 

RETURN 

CONTINUE 

PREPARE  SKEWNESS  GRAPH 

I F 1 . NOT .GDAUTO)  CALL  GRMODEI’ ALPHA’ ,0) 

IF (GDAUTO. AND. IZETA.EQ.O)  CALL  QRTEKQ 
IFIIZETA.EQ.O. AND. .NOT. GDAUTO)  READ ( I  OR . 1 8 )  IJK 
CALL  GRSAV I  SAVE ) 

IFIIZETA.EQ.O)  CALL  GRMODE t ’ ERASE ’ , 3 ) 

IF  I IZETA.EQ.O)  CALL  GR»nDE( ’ERASE’ ,3) 

IF! IZETA.GT.O)  CALL  GF  '  TT  !1 1  .  .  8 . 5 , -1  > 

CALL  QRHV ( ’ V * ) 

CALL  GRAXISIO.I .3,1 .5,5. ,0. ,0. , 1 . ,6,2, ’WIDTH* ,5,2) 

CALL  GRAX I S (3 , 1 . 5 , 1 . 5 , 7 . 5 , 90 . , 0 . , 2 . , 1 1 , 2 ,’ SKEWNESS’ ,8,2) 

CALL  GRP  1C ( 1.5, 1.5, 6. 5,9.) 

CALL  GRSCL ( 0 . ,0. ,1 . ,2. ) 

CALL  GRSXYI4. , 9 . 5 , X . Y ) 

CALL  GRTXT IBUF ,NB) 

CALL  GRPRNT (0 , X.O.Y, .15,0.) 

CALL  GRSXY (6.8,8.5,X,Y) 

CALL  GRSXY !6.B,8.3,X,DY) 

DYsY-DY 

CALL  GRTXT (’FREQUENCY  ’,30) 

CALL  GRNUMI’F’ ,-1,2, FREQ) 

CALL  GRTXT!’  GHZ’ ,4) 

CALL  GRPRNT ( -1 , X , 0 , Y , . 08 , 0 . ) 

YsY-DY 

CALL  GRTXT! ’POLARIZATION  ’,50) 

CALL  GRTXT ( POL , 3 J 

CALL  GRPRNT ( -1 , X , 0 , Y , . 08 , 0 . ) 

Y=Y-DY 

CALL  GRTXT! ’X  SHADOWED  ’,20) 
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1 1 1 0> 

1 1 20> 

1 1 3  0  > 

1 1 40> 

1 1  50> 

1 1 60> 

1 1 7  0> 

1 1  6  0  > 

1 1 90> 

1 200> 

1 2 1 0  > 

1 220> 

1 230> 

1 240> 

1  2  9  0  > 

1  2  6  0  > 
1270> 

1 2  8  0  > 

1  290> 

1 300> 

1  3 1 0  > 

T  3  2  0  > 

1  3  3  0  >  C 
1  340> 

1 350> 

1 3  8  0  > 

1 370> 

1 380> 

1 390>  10 
1 400> 

1  4 1  0  > 

1 420>  C 
1430> 

1 440> 

1 450> 

1 460> 

1 470>  14 
1480>  C 
1 490>  C 
1 500>  C 
1  503  > 

1 51 0> 

1  5 1 2  > 

1 51 4>  18 
1  5  2  0  > 

1 530> 

1 540>  C 
1  5  5  0  > 

1  560> 

1  570> 

1 5  8  0  >  C 
1 580> 

1  600> 

1 51 0> 

1 820>  C 
1  6  3  0  > 

1 6  40  > 

1 8  5  0  > 

1 680> 

1 870> 


CALL  QRNUMI ' F * ,-1 ,2,PERSH) 

CALL  QRPRNTI-1 , X , 0 , Y , . 08 , 0 . ) 

Ysy-DY 

CALL  QRTXT I ’ PITCH  ’ , 20) 

CALL  QRNUMI ’ F ’ ,-1 ,2, PITCH) 

CALL  QRPRNTI-1 ,X,0,Y, .08,0.) 

YSY-DY 

CALL  QRTXT I ’ ROLL  ’ ,20) 

CALL  QRNUMI ’F’ ,-1 .2.R0LL) 

CALL  QRPRNTI-1 ,X , 0 , Y , .08 ,0 . ) 

YsY-OY 

CALL  QRTXT I' YAW  ' ,20) 

CALL  QRNUMI ' F ’ ,-1 ,2, YAW) 

CALL  QRPRNTI-1 ,X,0,Y, .08,0.) 

YsY-DY 

CALL  QRTXTl’TILT  ' ,20) 

CALL  QRNUMI ’F’, -1,2, TILT) 

CALL  QRPRNTI-1 ,X,0,Y, .08,0.) 

YsY-OY 

CALL  QRTXT I ’ PHIR  ' ,20) 

CALL  QRNUMI 'F’ ,-1 ,2,PHR) 

CALL  QRPRNTI-1, X,0,Y, .08,0.) 

CALL  QRMOVEIO. ,0. ) 

DO  10  Jsl.50 
TsFLOAT I J ) / 2 . 5 
WsT»IT*2. J/IT+1 .1**2 
Qs2.*T»IT+3.)/l(T*1 .)*IT*2.)) 

CALL  QRDRAWIW.Q) 

CALL  QRDRAWI1 . ,2.1 
CALL  QRDRAWI .5,0. ) 

DO  14  Jsl.NS 

CALL  QRSYM I SLOBE 1 2 , J ) , SLOBE 1 3 , J ) ,  .1,0. ,42) 

CALL  QRTXT  I  ’  ’  ,1) 

CALL  QRNUMI’ I’ , -1 , I F IX I SLOBE 1 1 , J ) ) , DUMMY ) 

CALL  QRPRNTI-1 , SLOBE 12 , J ) , 0 , SLOBE 1 3 , J ) , .08,0.) 

PREPARE  KURTOSIS  QRAPH 

I F I . NOT . QDAUTO)  CALL  QRMODE I ’ ALPHA ’ , 0) 

I F (QDAUTO. AND . IZETA.EQ.O)  CALL  QRTEKQ 
IFIIZETA.EQ.O. AND. .NOT. QDAUTO)  READ ( I OR . 1 8 1  UK 
FORMAT  I A2) 

IFIIZETA.EQ.O)  CALL  QRMODE I ’ ERASE*  , 3) 

IFIIZETA.QT.OJ  CALL  QRPLOT (8.5,0. , -1 ) 

CALL  QRAXISI0.1 .5,1 .5,5. ,0. ,0. ,1 . ,6 ,2 , ’WIDTH’ ,5.2) 

CALL  QRAX 15(3,1 . 5 , 1 . 5 , 7 . 5 , 90 . , -2 . , 6 . , 8 , 2 , ’ KURTOS I S ’ ,8,2) 
CALL  QRSCLIO. ,-2.  ,1  .  ,6.  ) 

CALL  QRSXYI4. , 9 . 5  , X . Y) 

CALL  QRTXT IBUF , NB) 

CALL  QRPRNT(0,X,0,Y, .15,0. 1 

CALL  QRSXY I6.8,8.5,X,Y) 

CALL  QRSXY (6.8,8.3,X,DY) 

DYSY-DY 

CALL  QRTXT (’FREQUENCY  ’,30) 

CALL  QRNUMI ’F’ ,-1,2, FREQ) 
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1 680> 

1 690> 
J700> 
1710> 
1720> 
t  730» 

1 740> 
1750> 

1 760> 

1 770> 

1 780> 

1  790> 

1 800> 

1 81 0> 

1 820> 

1 830> 

1 840> 

1 850> 

1  860> 

1 870> 

1  880> 

1 890> 
1900> 

1 9 1  0  > 

1 920> 

1 9  3  0  > 

1 940> 

1 950> 

1 960> 
1970> 

1 980>  C 
1 990> 
2000> 

201 0> 
2020> 
2030> 
2040>  15 
2050> 
2060> 
2070> 
2080> 
2090> 

21 00> 

21 1 0> 

2 1 2  0  > 

21 30>  18 
21 40> 

2 1  SO >  C 

2ieo> 

21 70> 

2 1 8  0  > 

21 90> 
2200>  17 
2205> 

2  2 1 0  > 
2220> 
2230> 


CALL  QRTXT ( ’  GHZ’ ,4) 

CALL  GRPRNTI-1, X.O.Y, .08,0.) 

Y=Y-DY 

CALL  GRTXT I ’POLARIZATION  ’,50) 

CALL  GRTXTIPOL ,31 

CALL  GRPRNTI-1 ,X.O, Y, .08,0. ) 

Y*Y-DY 

CALL  QRTXT I ’%  SHADOWED  ’,20) 

CALL  GRNUMI’F’ ,-1,2 , PERSH) 

CALL  GRPRNTI-1 ,X,0,Y, .08,0.) 

YsY-DY 

CALL  GRTXTI’PITCH  ’ ,20) 

CALL  GRNUMI’F' ,-1 ,2, PITCH) 

CALL  GRPRNT 1-1, X,0,Y, .08,0.) 

Y=Y-DY 

CALL  QRTXT I ’ ROLL  ' ,20) 

CALL  GRNUMI’F’ ,-1 ,2, ROLL) 

CALL  GRPRNTI-1 ,X ,0 , Y , . 08 , 0 . ) 

Y=Y-DY 

CALL  QRTXT (’YAW  ’ ,20) 

CALL  GRNUMI’F’ ,-1 ,2, YAW) 

CALL  GRPRNTI-1 ,X , 0 , Y , . 08 , 0 . ) 

Y=Y-DY 

CALL  QRTXT (’TILT  ' ,20) 

CALL  GRNUMI’F’ ,-1,2, TILT) 

CALL  GRPRNTI-1, X,0,Y, .08,0.) 

Y=Y-DY 

CALL  QRTXT I ' PH IR  ’ ,20) 

CALL  GRNUMI’F’ ,-1,2,PHR) 

CALL  GRPRNTI-1 ,X,0,Y, .08,0. ) 

CALL  GRMOVEIO. ,0.) 

DO  15  J*1 ,50 
T=FLOAT(J)/2.5 
WST»IT*2.)/IT*1.)N»2 
Q:6.*T#IT+4. ) / I T ♦ 2 . ) MM2 
CALL  GRDRAWIW.G) 

CALL  GRDRAWl 1 . ,6. ) 

Wsl. 

DO  16  J  =  1 ,50 

WsW-.OI 

ZS4.NW-2. 

ZN= (Z+2 . ) / 12 .NW-Z) 

ZN2*ZN*ZN 

Gs6 . *ZN2tt (W+1 . / (ZN2*W) -3 . /ZN-1 . /ZN2) / l(ZN*1.)*(ZN+2.)) 
CALL  GRDRAWIW.G) 

CALL  GRDRAWl 0. ,-1 .5) 

DO  17  Jal , NS 

CALL  QRSYMI SLOBE (2  ,  J )  , SLOBE 1 4 , J ) , . 1 , 0 . , 42) 

CALL  QRTXT I ’  ’  ,1 ) 

CALL  QRNUMI ’ I ’ , -1 , 1 F I X I SLOBE 1 1 , J ) ) , DUMMY) 

CALL  GRPRNTI-1 , SLOBE (2 , J ) , 0 , SLOBE 14 , J ) , .08,0.) 

CALL  GRRESISAVE) 

CALL  QRHV(’H’) 

RETURN 

END 


>>>  FILE  «  SYO : [ 1 1 1 , 1 2 ] SPAWNS . FTN  27-SEP-83  14:36:42  <<< 

10>  SUBROUTINE  SPAWNS 

20>  INTEGERM2  EXSTAT 1 8 ) 

30>  DATA  TASK/6RMCR. . ./ 

40>  CALL  SPAWN  I  TASK, , ,11 , .EXSTAT, , ’RUN  LOBE/TASK=LOBE' ,18, , IDSW) 

50>  CALL  STOPFRI1 1 ) 

80>  RETURN 

70>  END 
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>>> 


FILE  s  SYO : t 1 1 1 , 1 2] LOBE .FTN 


10 

20 


30>  C 

L 

000 

BBSS 

EEEEE 

40>  C 

L 

0  0 

B  B 

E 

50>  C 

L 

0  0 

BBBB 

EEEE 

60>  C 

L 

0  0 

B  B 

E 

70>  C 

LLLLL 

000 

BBBB 

EEEEE 

80 

90 

100 

110 

120 

130 

140 

ISO 

160 

170 

175 

180 

190 

200 

210 

220 

230 

240 

2S0 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 


C« 

C 


27 -SEP-83  14:36:48  <<< 


C 

C 

C 

C 

C« 

C 


COMPUTE  RCS  LOBE  STATISTICS 


REALM  S L OBE (4,181 ) ,S 10(4,1 80) , NF , ND 1 1 81 ,3) .SAVE (30) 

REALM  THETA (41 .ALPHA  14) , DELTA ( 4 ) .SCALE  16) .BI2.2) .CI2.2) 
REALM  012,2,2,2) .AMNI2) ,DRS (2) , TMPARR 12 , 2 ) , QSIZEI5) 
INTEGERS  QLINT.RR.SS 
LOQICALK2  QDAUTO 
BYTE  BUF (80) 

DID(X):44S77. NCOS  1 1 . 5707963-ATAN IX/685 . 8 ) ) 

IOW= 2 

CALL  ASNLUNI IOW, ’ T I  * , 0) 

IOLS1 

CALL  ASNLUNI IOL, ’SY' ,0) 

IOT  =  5 

CALL  ASNLUNI IOT, ’ SY‘ .0) 

OPEN  I  UN IT= IOL, NAMES’ LOBE. TMP’ .TYPEs’OLD* .RECORDS  I ZE  =  1 81 , 

1  ACCESS-DIRECT' .READONLY, SHARED) 

READ  I IOL ’ 1 )  I  OUT  ,  KLOBE  .HUNT,  NRCS ,  NLOBE ,  THETA  .ALPHA ,  DELTA , 

1  BETA, SCALE. IZETA, I PCOOE.WLEN.LBUQ.NB, BUF. QDAUTO. 

2  FREQ.POL .PERSH, PITCH, ROLL , YAW, TILT, PHR.QS I ZE 
LBUG4=LBUG. AND . 4 

LBUG8=LBUG . AND . 8 


362 

364 

366 

368 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 


C 

IF  t  DELTA  1 3 ) .EQ.O. )  WRITE! I OW , 8 ) 

8  FORMAT  I ’  LOBE>  --  DELTA  HAS  NOT  BEEN  DEFINED’) 

IFIDELTAI3) .EQ.O.)  QO  TO  38 
C 

IF  I LBUQ4. NE . 0)  WRITEI2.227)  IOUT , KLOBE ,GL I  NT .NRCS .NLOBE . IZETA , IPCODE . 
1  THETA, DELTA, ALPHA .SCALE, BETA, WLEN 
227  FORMAT I1X,7I4/1X,8F8.2/1X,8F8.2/1X,8F8.2) 

C 

1 0= I OW 

IF! I0UT.EQ.1 )  QO  TO  4 
C 

C  INITIALIZE  LOBE  FOR  PRINTED  OUTPUT 
C 

IF! IPCODE. EQ.O)  QO  TO  7 
10=3 

CALL  ASNLUNI 10, ’QD’ ,0) 

QO  TO  7 
C 

C  INITIALIZE  LOBE  FOR  GRAPHICAL  OUTPUT 
C 

4  100*4 
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550> 

560> 

570> 

580> 

590> 

c 

600> 

c 

610> 

c 

620> 

c 

630> 

c 

640> 

c 

650> 

c 

660> 

c 

670> 

680> 

7 

690> 

222 

700> 

C 

710> 

C 

720> 

730> 

740> 

750* 

760> 

770> 

780> 

790> 

C 

8  00> 

C 

8 1 0  > 

C 

820> 

C 

830> 

C 

840> 

C 

850> 

860> 

C 

870> 
880> 
890> 
900> 
91 0> 

21 

920> 

C 

930> 

940> 

950> 

960> 

22 

970> 

228 

980> 

990> 

C 

1000> 

C 

1010> 

1 020> 

1  030> 

C 

1 040> 

C 

1050> 

C 

1 060> 

1 070> 

1 080* 

1 090> 
1100> 

C 

1 1 1 0> 

C 

*«*»» 

OOO 

V  V  V 

1 

CALL  GRINITIO,  IZETA,  100,  IER) 

CALL  GRHV(’H’) 

CALL  GRPICIGSIZEH)  .GSIZEI2)  .QSIZE  (3)  ,GSIZE(4J) 
IF  I IZETA. EQ.O)  CALL  GRDEVI’TEK’ , ’HARD' , IER) 


LBUQ  DEFINITIONS: 

LBUQsl a>TOLOBE  TRANSFORMATIONS 
LBUGs2s>T0L0BE  OUTPUT 
LBUGs4*> INPUT  A  LOBE  OUTPUT 
LBUGs8s>GLINT  OUTPUT 

NANGs  (DELTA  (2)  -DELTA  ID)  /DELTA  1 3 1  ♦  1 . 0001 
IFIIOUT. EQ.O. OR. IZETA. NE.O)  WRITE! 10,222)  NANG 

FORMAT  I ’ 0>>  COMPUTE  RCS  LOBE  STATISTICS  FOR’,13.’  THETA  ANGLES  < 

COMPUTE  CONSTANT  TERMS 

SNBET=SIN(BETA) 

CSBET1 =COS (BETA) +1 . 

SOL= . 2998E9 
Plsj. 1415926 
TWOP I =P 1*2 . 

RADsPI/180. 

IFIIOUT. EQ.1)  CALL  GRSCLISCALEM ) .SCALE (2) .SCALE  13) .SCALE (4) ) 


PERFORM  ENTIRE  SET  OF  LOBE  STATISTIC  COMPUTATIONS 
(READ  ENTIRE  TEMPORARY  FILE  EACH  TIME) 

FOR  EACH  REQUIRED  DELTA  ANGLE 

ANQ8DELTAI1) 

IF (ANG .QE .THETA (1 ) )  GO  TO  22 

ANG8ANG*DELTA(3I 

NANQsNANG-1 

IF(NANG.LE.O)  STOP  ’LOBE  ERROR’ 

GO  TO  21 

FIRSTA=ANG 

SLOPE  =  FLOAT INRCS-1 1/AMAX1 I .01 .THETA (2) ‘THETA  1 1 ) ) 
CONSTsI .‘SLOPEmTHETAID 

IF (LBUG4.NE .0)  WRITEI2.228)  SLOPE .CONST, ANG 
FORMAT!’  SLOPE, CONST, ANGs’ .3F10.3) 

DO  100  LLsI.NANQ 

DETERMINE  NEXT  INDEX  INTO  LOBE  ARRAY 
I F (ANG .GT . THETA  1 2 ) )  GO  TO  101 
L= I F I X I SLOPEMANG+CONST ♦ . 001 ) 

COMPUTE  SUM  AND  STANDARD  DEVIATION 

NLBs-3 
SUM=0 . 

STDsO, 

N*0 


NLB=NLB+4 

IFINLB.GT.NLOBE)  GO  TO  3 

READ ( IOL ’ NLB  +  1 )  ( SLOBE 1 1 , J ) , J  =  1 , L ) 


<’) 
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1 1  50> 

1 F ( SLOBE ( 1 , L 1 . EQ.Q . J  GO  TO  1 

1 1 60> 

DO  19  K  =2 , 4 

1 1 70> 

19 

READ ( I0L’ NLB+K)  ( SLOBE (K , J ) , J =1 , L) 

1 1  80> 

C 

1 1 90> 

IF ILBUG4.NE .0)  «R 1 TE (2 , 221 )  (SLOBE ( 1 , L) , 1 *1 , 4) 

1  200> 

221 

FORMAT!’  RCS.POSHs' ,4E11 .3) 

1 21 0> 

C 

1 220> 

N=N*  1 

1  230> 

DO  2  J  =  1 ,4 

1 240> 

2 

S IG ( J , N) = SLOBE ( J , L ) 

1 2  5  0  > 

C 

1 260> 

$UM=SUM+SIGI1 ,N) 

1 270> 

STD=STD*SIG(1 ,N)N«2 

1 280> 

IFIN.LT.180)  GO  TO  1 

1290> 

IFtN.EQ.180)  WRITEI2.9) 

1 300> 

9 

FORMAT ( ’  RCS>  —  LOBE  WILL  USE  ONLY  THE  FIRST  180  COMPONENTS’) 

1 31 0> 

C 

1  3  2  0  > 

C 

1  330> 

3 

IF(N.LE.I)  GO  TO  100 

1  3  40  > 

IFKIOUT.EQ.O. OR. IZETA.NE.01.AND.LL.EQ.il  WRITEI2.24)  N 

1  3  5  0  > 

24 

FORMAT ( 10X , ' >>  USING’, 14,’  SUBCOMPONENTS  <<’> 

1 3  0  0  > 

STD=SQRT ( SUM*SUM-STD) 

1 370> 

IF  1 (LBUG4! . NE .0)  WRITEI2.223)  N, SUM, STD 

1  3  8  0  > 

223 

FORMAT!'  N, SUM, STDs’ , I5.2E11 .3) 

1 3  9  0  > 

C 

1 400> 

C 

1  41  0> 

C 

COMPUTE  R  C  S  LOBE  1 

1  4  2  0  > 

C 

1  43  0  > 

C 

1 440> 

GsO. 

1  450» 

N1=N-1 

1 40O> 

C 

1 470> 

DO  5  1 =1 ,N1 

1480> 

11=1*1 

1  490> 

DO  5  J  = 1 1  ,N 

1  500> 

S 

G  =  G*SIG(1  ,  imSIGII  ,  J)#(  (SIGI4,  1  )-SIG(4,  J)  1 NCSBET1  ♦ 

1 5 1 0  > 

9  (SIGI2, H-SIGI2, J))NSNBIT)»»2 

1  S  2  0  > 

NFs SORT (Q) / (SOLMSTD) 

1  53  0> 

IF  1 (LBUG4) .NE .0)  WR 1 TE ( 2 , 224)  NF 

1 540> 

224 

FORMAT!’  NFS’ ,4E11 .31 

1  550> 

NF=NF«1 . E9 

1  5  8  0  > 

C 

1570» 

C 

1  5S0> 

C 

COMPUTE  R  C  S  LOBE  1 1 

1  590> 

C 

16  00> 

C 

1 6 1 0  > 

C 

1 620» 

C 

COMPUTE  A  MATRIX  (A ( 1 , 2 1 :A (2 , 1 ) I 

1 630> 

C 

1  8  40  > 

A1 1 sQ . 

1  6  5  0  > 

A12*0. 

1 860> 

A22*0. 

1 670> 

C 

1  8  8  0  > 

DO  6  1=1 ,N1 

1 690> 

11=1*1 

1  700> 

DO  6  J  =  1 1  ,N 

1 7 1 0  > 

C 

1 720> 

TIsSIGM  ,  1 ) MS  1 G 1 1 ,  J ) 

1 730> 

C 

1 740> 

T2= ( S IG (2 .  1  )-SIG(2,J>  JMCSBETt- (5IQ(4,  1 J-SIQI4, J) IftSNBET 
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1 7  50  > 

C 

1 760> 

T3= (S IQ(3 , 1 ) -S 1 Q (3 , J 1 IftCSBETI 

1 770> 

C 

1 780> 

A1 1 =A1 1 ♦T1*T2*T2 

1 790» 

A12=A12*T1*T2*T3 

1 800> 

A22=A22+T1«T3«T3 

1 8 1 0  > 

C 

1 820> 

6 

CONTINUE 

1 830> 

A21=A12 

1 840> 

C 

1 850> 

C 

1 860> 

C 

COMPUTE  F (DELTA)  VS  DELTA 

1 870> 

C 

1 880> 

DELT=360./DELTA<4) 

1 890> 

C 

DELT=2. 

1900> 

DEL*0. 

1  9 1 0  > 

IDELsl 

1 9  2  0  > 

1 F ( LBUQ4 . NE . 0 )  WR 1 TE 1 2 . 2 2 5 1  DELT , A1 1 ,A1 2 , A21 , A22 

1 930> 

223 

FORMAT!’  DELTs’ , Ell .3/’  A* ’ , 2E1 1 . 3/3X , 2E1 1 . 3) 

1940> 

C 

1 850> 

10 

IF (DEL .QT .339 .9)  QO  TO  11 

1 960> 

ND(IDEL,1)sDEL 

1 970> 

D1=SIN(DELmRAD) 

1 980> 

D2= -COS (DEL MR ADI 

1  990> 

T1=A1H»D1*A12«D2 

2000> 

T2=A21<*D1  +A22MD2 

201 0> 

A=D1«T1*D2*T2 

2020> 

ND( IDEL ,2) sSQRT 1 A) / (WLENMSTD) 

2030> 

C 

2040> 

DEL=DEL*DELT 

2050> 

1  DEL  =  1  DEL* 1 

2000> 

QO  TO  10 

2070* 

C 

2090> 

11 

ND 1 IDEL , 1  1 =360 . 

2090> 

ND( 1  DEL ,2) =ND1 1 ,2) 

21  00> 

C 

21 1 0> 

DO  18  J=1 , IDEL 

21 20> 

18 

ND(J,1)=AMODIND(J,1)*90.  ,360.1 

21 30> 

C 

21 40> 

IFIIOUT.EQ.il  QO  TO  IS 

21  50> 

C 

21  60> 

C 

PRINT  NF  AND  ND 

21 70> 

C 

2 1  8  0  > 

WRITE! 10,12)  ANQ.NF 

21 90> 

12 

FORMATCORCS  LOBE  OUTPUT’  ,  5X ,’  THETAs ’, F6 . 1 ,  ’  MF (GHZ)  =  ’  ,F8 . 3/ 

2200> 

8  5X, ’DELTA’ ,5X,’ND’/J 

2  2 1  0  > 

C 

2220> 

DO  13  3=1 , IDEL 

2230> 

13 

WR ITE ( 10, 1 4)  NDIJ.I) , ND 1 J , 2 1 

2240> 

14 

FORMAT (1X,F9.2,F10.3) 

2250> 

QO  TO  100 

2280> 

C 

2270> 

C 

GRAPH  NF  AND  1/ND 

2280> 

C 

2290> 

IS 

CONTINUE 

2300> 

FMAX=0 . 

2  3 1  0  > 

FMIN=1000. 

2320> 

DO  16  J=1 , IDEL 

2330> 

ND I J , 2  I =1 . / ( ND ( J , 2 IttRAO) 

2340> 

FMAX=AMAX1 (FMAX.NDI J ,2)1 
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2950 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2935 

2940 

2942 


I F ( FMAX . EQ . ND I J , 2 ) )  FAMQsND  I J  ,  1 ) 

FMI  N=AMI  N1  IFMIN.NDN.2n 
NDC J , 21 =  ND  I J ,21 

16  CONTINUE 

I F  IFANQ .QE .160.)  FANQ  =  FANQ-1 80 . 

IFIFANQ.QT.90. )  FANG=FANG-180. 

C 

CALL  GRRXY ( ANQ , 0 . ,X,Y) 

YsQS I ZE  1 2 ) -  .  95 

CALL  QRSAV I  SAVE  I 

CALL  QRINIT11 ,0,0, IER) 

IFILL.NE.1)  GO  TO  17 
C 

CALL  GRTXT ( ’ NF (GHZ) ’ , 0) 

CALL  GRPRNT I +1 , 2 . 1 , -1 , Y* . 53 .  .  1  , 0 . ) 

CALL  GRTXT I 'MAJOR' ,0) 

CALL  GRPRNT 1*1 , 2 . 1 , - 1 , Y* . 38 , . 1 , 0 . ) 

CALL  GRTXT ( ’MINOR’ ,0) 

CALL  GRPRNT (  +  1 , 2 . 1 . -1 , Y+ . 23 . . 1 , 0  . ) 

CALL  GRTXT I ' ANGLE ’ .0) 

CALL  GRPRNT (*1,2.1 , -1 , Y* . 08 , , 1 , 0 . ) 

C 

17  CALL  GRNUMI ' F ' , -1 , 3 , NF I 

CALL  GRPRNT(0,X,-1 ,  Y* . 53 , .10,0. I 
CALL  GRNUMI ’ F ’ ,-1 ,2, FMAX) 

CALL  GRPRNT (0 , X,-1,Y+.38,.1,0.) 

CALL  GRNUMI ' F ' , -1 , 2 , FM I N ) 

CALL  GRPRNT (0 , X , -1 ,  Y+ . 23 , .1 ,0.) 

CALL  GRNUMI ’ F ’ ,-1 ,2, FANG) 

CALL  GRPRNT(0,X,-1,Y+. 08, .1,0.1 
C 

CALL  GRCP I ’ POLAR' I 
YsGSIZEm 

CALL  GRPIClX.Y.Xtl  .  ,Y*1  .) 

CALL  GRSCLIO. ,0. , tSCALE 13) -SCALE! 1 11/6.5 ,360. ) 

CALL  GRP ICI0.,Y-.28,11.,Y*.5) 

CALL  GRWI ND I  -  I  DEL ,ND(1 , 2 ) , HD  1 1 , 1 ) , I DUM , 0 , DUM, DUM, I DUM) 
CALL  QRUPPN 
CALL  GRRES I  SAVE ) 

C 

100  ANG=ANG*DELTA(3) 

101  CONTINUE 
C 

IFIGLINT.EQ.O)  CALL  EXIT 
C 

C 

c 

C  COMPUTE  GLINT  LOBE  III 

C 

C 

OPEN (UN I T= IOT , NAMES’  LOBE . SCR’ .TYPEs’SCRATCH’ ) 

C 

DO  38  I A L P - 1  ,91  ,90 
C 

C  PAGE  PLOT 

IF! IOUT.EQ.O)  GO  TO  50 

IF  I . NOT.QDAUTO)  CALL  GRMODE I ’ ALPHA ' , 0 ) 

IF  I IZETA.EQ.O. AND.QDAUTO)  CALL  GRTEKG 

IF! IZETA.EQ.O. AND. .NOT.QDAUTO)  REAO(IOW,47)  UK 
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2944> 

47 

FORMAT  I A2I 

2950> 

IFI IZETA.QT.O)  CALL  GRPLOT (11. ,9. 5,-1) 

2960> 

CALL  GRIN1T11 ,0,0, IERI 

2970> 

C 

2990> 

50 

ALP=FLOATI IALP-1) 

2990> 

REWIND  I0T 

3000> 

C 

301 0> 

C 

LOOP  ON  THETA  INCREMENTS 

3020> 

ANG=THETA 1 1 ) 

3030> 

DO  35  Lsl ,NRCS 

3040> 

IFIALP.EQ.O. )  GO  TO  23 

30S0> 

C 

3060> 

READ! IOT.20)  N, XBAR , YBAR , SUM, B 

3070> 

IFIGLINT.EQ.1 .OR.N.LE.1)  GO  TO  37 

3080> 

GO  TO  39 

3090> 

C 

31  00> 

23 

NLBs-3 

3 1 1 0> 

XBARsO. 

3 1 2  0  > 

YBARSO. 

31 30> 

SUM=0. 

31 40> 

STD=0. 

31  S0> 

N:0 

3 1 60> 

C 

31 70> 

31 

NLB=NLB*4 

3 1  80> 

IF INLB.GT.NLOBE)  GO  TO  30 

3 1 9  0  > 

READ (IOL* NLB+1 1  (SLOBE (1 , J) , J=1 ,L) 

3200> 

1 F ( SLOBE ( 1 , L ) . EQ . 0 )  GO  TO  31 

32 1 0> 

DO  32  K=2 , 4 

3220> 

32 

READ  1 IOL ’ NLB*K)  1 SLOBE IK , J ) , J  =  1 , L ) 

3230> 

C 

3240> 

N=N*1 

3250> 

DO  33  J s 1 ,4 

3260> 

33 

S IG ( J , N) s SLOBE (3 , L I 

3270> 

C 

COMPUTE  E [ • i gma ] 

3280> 

SUM=SUM*SIG(1 , N) 

3290> 

C 

COMPUTE  SUMtaigmaNgraak-X) 

3300> 

XBAR=XBAR*S IQ  I 1 , N ) «S 1 Q ( 2 , N ) 

331 0» 

C 

COMPUTE  SUMIa igma«grask-YI 

3320> 

YBAR=YBAR*S IG I 1 ,N)tfSIG(3,N) 

3330> 

GO  TO  31 

3340> 

C 

3350> 

30 

IFIN.EQ.O)  WRITE! IOT , 20)  N 

3360> 

IFIN.EQ.OI  GO  TO  37 

3370> 

C 

COMPUTE  grtak-Xbar 

3390> 

XBARsXBAR/SUM 

3390> 

C 

COMPUTE  graak-Ybar 

3400> 

YBARsYBAR/SUM 

341 0> 

C 

3420> 

STD=0 . 

3430> 

IFIGLINT.EQ.1 .OR. N.EQ.1)  WRITE  1 IOT. 201 

3440> 

IFIGLINT.EQ.1 .OR. N.EQ.1 )  GO  TO  37 

3450> 

C 

3480> 

BI1 ,1  Iso. 

3470> 

B 1 1  ,21*0, 

3490> 

BI2,1)sO. 

3490> 

8(2,21=0. 

3500> 

C 

351  0> 

N1=N-1 

3520> 

DO  34  1=1 ,N1 

3S30> 

11=1*1 
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3540> 

3550> 

C 

3560> 

3570> 

3580> 

3590> 

C 

3600> 
361  0> 
3620> 
3630> 
3640> 

34 

3650> 

3660> 

20 

3670> 

C 

3680> 

39 

3690> 

3700> 

C 

3  7 1 0  > 

C 

3720> 

C 

3730> 

3740> 

3750> 

3760> 

226 

3770> 

3780> 

3790> 

C 

3800> 

37 

381 0> 
3820> 
3830> 

C 

3840> 

3850> 

3860> 

C 

3870> 

36 

3880> 

C 

3890> 

35 

3900> 

C 

391  0> 
3920> 

C 

3930> 

C 

3940> 

C 

3950> 

3960> 

53 

3970> 
3980> 
3990> 
4000> 
401 0> 

54 

4020> 

55 

4030> 

4040> 

C 

4050> 

C 

4060> 

C 

4070> 

51 

4080> 

C 

4090> 
41 00> 
41 05> 
4 1 1  0  > 
4 1  2  0  > 

DO  34  J s 1 1 ,N 

T1 =S I G 1 1 , I  IMS  IQ ( 1 , J) 

T  2  =  S I Q ( 2 , I )+SIQ(2,JI-2. MXBAR 
T  3  =  S I Q ( 3 , I )+SIG(3,J)-2. MYBAR 

BI1  ,11=8(1 ,1 J+T1MT2MT2 
B(1  ,21=8(1 , 2 ) ♦T1MT2MT3 
8(2.11=811  ,2) 

8(2,21=8(2 .2) +T1MT3MT3 
CONTINUE 

WRITE! IOT.20I  N.XBAR, YBAR, SUM, B 
FORMAT (A2, 7 A4) 

A1 =COS I ALP ) 

A  2  =  S I N ( ALP  1 

COMPUTE  GLINT  STANDARD  DEVIATION 

BB=B (1,1) MAI MAI +B ( 2 , 1 1 MA2MA1 +B (1 , 2 ) MAI MA2+B ( 2 , 2 ) mA2mA2 
STD= SORT iBB) / <2  MSUM) 

I F ( LBUQ8 . NE . 0 )  WRITE12.226)  N , SUM , STD , XBAR , YBAR , ANQ , B 
FORMAT  CON,  SUM,  STDs’  ,  I3.2E11  .3/ 

1  '  XBAR , YBAR , ANQ= ’ , 3E1 1 . 3/ 

1  '  B=’ ,4E11 .31 

ND (L , 1 ) sANQ 
ND ( L , 3 ) sSTD 

I F (ALP . EQ. 90 . )  QO  TO  36 

ND ( L , 2 ) =XBAR 
QO  TO  35 

ND ( L , 2 ) = YBAR 

ANQ  sANQ*' THETA  ( 3  J 

IF ( IOUT .EQ. 1 )  QO  TO  51 

PRINT  QLINT  DISPLACEMENT  VECTOR 

WR I T E ( 10,53)  ALP 

FORMATCORCS  QLINT  OUTPUT  ALPHA*’ ,F6 . 1/ 

1  '  ANQLE  DISPLACEMENT  D I S- 3  DIS+S'/) 

DO  54  LP=1 , NRCS 
T1=ND(LP,2)-ND(LP,3) 

T2=ND(LP ,2) >ND(LP ,3) 

WRITE! 10,55)  ND l LP , 1 ) , ND ( LP , 2 ) , T1 , T2 
FORMAT (F6. 1 ,F11.3,F10.3,F10.3) 

QO  TO  52 

GRAPH  QL INI  DISPLACEMENT  VECTOR 

IFIIZETA  EO.O)  CALL  QRMODE C ERASE ' , 3 ) 

CALL  QRPLOT IQS  I ZE (3 ) ,QS!ZE(2) ,0) 

CALL  QRPLOT (GSIZE(1),GSIZE(2),1) 

DAXsPj I ZE (3 ) -QS I ZE ( 1 ) 

CALL  QRAX I  $  1 0 , QS I ZE 1 1 ) , QS I ZE ( 2 ) - . 95 , DAX , 0 . , SCALE ( 1 ) , SCALE ( 3 ) , 
1  I F I X ( SCALE ( 5 ) I , I F I X ( SCALE ( 6 ) ), ’OBSERVATION  ANGLE', 99,1) 
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4160 

4165 

4170 

4160 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4260 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 


C 

DAY=GS IZE (4) -6SIZE (2) 

CALL  QRAX I S 1 3 , QS I ZE ( 1 ) . QS I ZE ( 2 ) , DAY . 90 . ,-360. ,360. ,5,2, 
1  ’DISPLACEMENT  IDEQI ' ,99,2) 

C 

IF( IZETA.EQ.O)  CALL  QRMODE C ALPHA’ , 2 ) 

CALL  QRTXT(BUF.NB) 

CALL  QRPRNT (0,5. 5,-1, 7., .15,0.1 

CALL  QRTXTI ’MONOSTATIC  QLINT  FOR  ALPHA:’, 0) 

CALL  QRNUMCF’  , -1,1, ALP) 

CALL  QRPRNT (0,5. 5,-1 ,6. 75,. 15,0.) 

IFI IZETA.EQ.O)  CALL  QRMODE (' ALPHA’ , 3 ) 

C 

C  DISPLAY  RADAR  PARAMETERS 
C 

X*9 . 0 
Y  =  7. 

DYs.2 

CALL  QRTXTI ’FREQUENCY  ’,30) 

CALL  QRNUMI ’ F ' ,-1 ,2, FREQ) 

CALL  QRTXTI'  GHZ' ,4) 

CALL  QRPRNT I -1 , X , 0 , Y , . 1 2 , 0 . ) 

YsY-DY 

CALL  QRTXTI ’POLARIZATION  ’,50) 

CALL  GRTXT (POL ,3) 

CALL  GRPRNTI-1 ,X,0,V, .12,0.  ) 

Y=Y-DY 

CALL  QRTXTI  X  SHADOWED  ’,20) 

CALL  QRNUMCF’  ,-1  ,2  ,  PERSH) 

CALL  QRPRNT I -1 , X , 0 , Y . . 1 2 , 0  . ) 

Y=Y-DY 

CALL  QRTXTI ’PITCH  ’ ,20) 

CALL  QRNUMCF’  ,-1,2  .PITCH) 

CALL  QRPRNT I -1 , X , 0 , Y , . 1 2 , 0 . ) 

Y=Y-DY 

CALL  QRTXT I ’ ROLL  ' ,20) 

CALL  QRNUMCF’  ,-1,2, ROLL) 

CALL  QRPRNT 1-1 ,X,0,Y, .12,0. ) 

Y=Y-DY 

CALL  QRTXT I ' YAW  ’ .20) 

CALL  QRNUMCF’  ,-1  ,2,  YAW) 

CALL  QRPRNT I -1 ,X,0,Y, .12,0. ) 

YSY-DY 

CALL  QRTXT (’TILT  ’ ,20) 

CALL  QRNUMCF’  ,-1  ,2, TILT) 

CALL  QRPRNT I -1 ,X,0,Y, . 1S,0. ) 

Y:Y-DY 

CALL  QRTXT I ' PH  I R  ’ ,20) 

CALL  QRNUMCF’  ,-1  ,2,PHR) 

CALL  QRPRNT I -1 ,X,0,Y, .12,0.  ) 

C 

C  DRAW  AVERAGE  DISPLACEMENT  CURVE 

CALL  QRP 1C (QS I ZE 1 1 ), QS I ZE 1 2 ), QS I ZE (3 ), QS I ZE 1 4) ) 

CALL  QRSCL I  SCALE  1  1  ) , -360 . , SCALE  I  3 ) , 360 . ) 

C  CALL  QRVEC I -NRCS , ND 1 1 , 1 ) , ND 1 1 , 2 ) , DUM, 0 , DUM, DUM, DUM) 

C 

CALL  QRMOVE I ND 1 1 , 1 ) , D I D I ND 1 1 , 2 )  ) ) 

DO  70  JRCS"2  NRCS 

70  CALL  QRDRAWIND I JRCS , 1  I , 0  I D I ND I JRCS , 2 ) ) ) 
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4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 


IFIQLINT.EQ.1 I  GO  TO  52 
C 

C  GRAPH  DISPLACEMENT  ♦  /-  STD 
CALL  GRMOVE (ND ( 1 ,1 ) ,0 
DO  49  1=2, NRCS 

49  CALL  G1DRAWINDIJ.1 ) ,D 

CALL  GRMOVE ( ND (1,11 ,D 
DO  59  J  =2 , NRCS 

59  CALL  GRDRAWIND1J  ,1  I  ,0 

CALL  GRUPPN 
CALL  GRRES ( SAVE ) 

C 
C 

C  COMPUTE  GLINT 

C 
C 

C  REPEAT  STD  CALCULATIONS  FOR  ELLIPSES 
C 

C  LOOP  ON  DELTA  INCREMENTS 
52  ANG  =  F  1 RSTA 

DO  45  LL=1 , NANG 
C 

C  DETERMINE  NEXT  INDEX  INTO  LOBE  ARRAY 
I F ( ANG . GT . THETA  12)1  GO  TO  38 
L=IFIX(SLOPE#ANG*CONST*.001 ) 

C 

NLB=-3 

XBARSO. 

YBAR=0. 

SUM=0 . 

STD=0 . 

N  =  0 
C 

41  NLB=NLB*4 

I F (NLB . GT . NLOBE )  GO  TO  40 
READ! IOL ’ NLB  +  1 )  ( SLOBE ( 1 , J ) , J  =  1 , L I 
I F ( SLOBE ( 1 , L 1 . EQ . 0 )  GO  TO  41 
DO  42  K  =  2 , 4 

42  READ ( IOL ’ NLB+K )  I SLOBE IK , J » , J =1 , L I 
C 

N  =  N*1 

DO  43  J=1 ,4 

43  SIGN  ,  N)  =SL08E  I  J  ,  L  ) 

C  COMPUTE  Elaigma] 

SUM=  SUM*  S I Q ( 1 , N ) 

C  COMPUTE  SUM(aigma*graek-X) 

XBAR=XBAR*S I G ( 1 , N)*SIG(2,N) 

C  COMPUTE  SUMl a i gmattg  r  eek - Y I 

YBAR  =  YBAR*S IG 1 1 , N)»SIGI3,N) 

GO  TO  41 
C 

40  IFIN.LE.1)  GO  TO  45 

C  COMPUTE  graak-Xbar 
XBAR=XBAR/SUM 
C  COMPUTE  graak-Ybar 
YBAR= YBAR/SUM 
C 

DO  48  MM=1 ,2 
DO  48  NN=1 .2 
B (MM, NN) =0 . 


IDINDI1  , 2 1 *ND (1,3)1) 

IDlNDtJ ,2)*ND(J ,3) ) I 
I  D ( ND ( 1 ,21-NDII ,3)1) 

IDlNDtJ , 2 ) -ND I J ,31 1  I 

LOBE  IV 
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5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

5610 

5620 

5630 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 


C(MM,NN)=0. 

DO  48  RR  =  1  ,2 
DO  48  SS-1  .2 

48  DIMM,NN,RR,S$)=0. 

C 

N1 =N-1 

DO  44  I =1 , HI 
11*1*1 

DO  44  J*I1  ,N 
C 

T1 =S IQ ( 1 , I IMSIQI 1 ,JJ 
T 2  =  S I G ( 2 , 1 1 *SIG (2 , J ) -2 . *XBAR 
T  3  =  S I G ( 3 , 1 1+SIG13, J)-2.*YBAR 

T4*((SIQ(4,IJ-SIG(4,J)  )*CSBET1  +  (SIG(2,I)-SIG(2,J) )*SNBET)*»2 
T5=(SIQ(2, I )-SIQ(2 , J ))»CSBET1 - (S IQ  14 , 1 ) -SIQ14 , J) IMSNBET 
T6*  C S I G ( 3 , I ) -S IQ (3  ,  J 1) ttCSBETI 
C 

TMPARRI1 ,1  )*T2*T2 
TMPARR (1,2) *T2#T3 
TMPARR (2,1) =TMPARR I 1 ,2) 

TMPARR 12,2) *T3#T3 
DO  44  MM* 1 ,2 
DO  44  NN=1 ,2 
TEMP*T1*TMPARR(MM, NN) 

B (MM, NN) =B (MM, NN) +TEMP 
C (MM , NN) *C (MM, NN) *TEMPtfT4 
C 

D (MM, NN , 1 ,1 )*D(MM,NN.1 ,1 )*TEMP»T5*T5 
D (MM, NN , 1 , 2 ) *D (MM, NN , 1 ,2) *TEMPNT5mT6 
DIMM, NN, 2,1 )*D(MM,NN,1 ,2) 
D(MM,NN,2,2)*D(MM,NN,2,2I*TEMP«T6KT6 
44  CONTINUE 

C 

A1 *COS (ALP  I 
A2=S I N ( ALP ) 

C 

C  COMPUTE  QLINT  STANDARD  DEVIATION 

BB=B ( 1 , 1 ) NA 1 MA 1 *B 1 2 , 1 J MA2NA 1 *B 1 1 , 2 )KA1 NA2 *B ( 2 , 2 1 «A2»A2 
STD* SORT  IBB) / (2 .mSUMI 

CC=C  < 1 . 1 ) *A 1 *A 1 *C 1 2 . 1 ) »A2»A 1 *C ( 1 , 2 ) #A 1 *A2  *C 1 2 , 2 ) ttA2«A2 
C 

C  COMPUTE  N(4# 1 1  a , f ) 

C 

NF*SQRT(CCI/(2.<fSTDft$0L) 

NF*NFM1 . E9 
C 

IF (LBUQ8 . NE .0)  WRITEI2.220)  N , SUM , STD , XBAR , YBAR , B , C , D 
220  FORMAT! ’ON, SUM, STD*’  ,  I3.2E11  .3/ 

1  ’  XBAR , YBAR* '  , 2E 1 1  .3/ 

1  ’  B= ’ , 4E1 1 .3/ 

1  '  C*’ , 4E1 1 .3/ 

1  4 ( ’  D*’  ,4E11  .3/)) 

C 

C 

C  COMPUTE  GLINT  LOBE  V 

C 

C 

AMNI1 ) *A1 
AMN(2)*A2 
DELT*360 . /DELTA  14) 
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5850 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 

6390 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 


DEL:0. 

I DE  L  —  1 
C 

60  I F (DEL , GT . 359 . 9 1  GO  TO  61 
ND ( I  DEL , 1 ) =DEL 

DRS I 1 ) =S I N I DELttflAD ) 

DRSI2) =-COS(DEL*RAD) 

C 

DD=0. 

DO  69  MM=1 ,2 
DO  69  NN=1 ,2 
DO  69  RR=1  ,2 
DO  69  SS  =  1  ,2 

69  DD:DD*D(M4,NN,RR,SS)»AMN(MM)*AMNINN)*DRS(RR)»DRSISS) 

C 

ND ( I  DEL , 2) =SQRT (DO ) / ( 2 . «WLEN«STD) 

C 

DEL=DEL+DELT 
I DE  L  = I DEL  +  1 
GO  TO  60 
C 

61  ND I  I  DEL , 1  ) =360 . 

ND I  I  DEL , 2) =ND 1 1 ,2) 

C 

DO  68  J:1 , IDEL 

68  ND ( J  ,  1 ) =AMOD I ND I J , 1 ) +90 . ,360.) 

C 

I F I IOUT.EQ.1 )  GO  TO  65 
C 

C  PRINT  NF  AND  ND 

C 

WRITE! 10,621  ANG , NF , ALP 

62  FORUATI’ORCS  GLINT  OUTPUT ’, 5X ,’ THETA2 ’, F8 . 1 , ’  NF (GHZ) * ’ ,F8 . 3/ 
S  5X. ’DELTA’ ,5X, ’ND1 , 5X , ' ALPHA2 ’ , F6 . 1 / 1 

C 

DO  63  J:l  ,  IDEL 

TEMP 2 1 . / (ND ( J , 2 ) NRAO ) 

63  WRITE! 10,64)  ND ! J , 1 ) , ND (J , 2 ) . TEMP 

64  FORMAT (1X.F9 .2 .2F10.3) 

GO  TO  45 

C 

C  GRAPH  NF  AND  1/ND 
C 

65  CONTINUE 
FMAX=0. 

FMIN=1000. 

DO  66  j:1 , IDEL 
NDIJ,2):1 ./INDIJ . 2 ) «RAD) 

FMAX=AMAX1 ! FMAX , ND 14 , 2 ) ) 

I F I FMAX . EQ. ND I J , 2 ) )  FANG  =  ND I J  ,  1 ) 

FMIN=AMIN1 IFMIN.ND! J , 2) ) 

C  ND! J , 2 ) 2ND ( J  ,  2  )  tt  .  1 

66  CONTINUE 

I F ( FANG . GE . 1 80 . )  FANQsFANG- 1 80 . 

IF (FANG. GT. 90.  )  F ANG  =  FANG- 1 80 . 

C 

CALL  GRRXY (ANG , 0 . ,X,Y) 

Y*GS IZE (2) - . 95 

CALL  GRINITM  ,0,0,  IER) 

IFlLL.NE.il  GO  TO  67 
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6550> 

C 

6560> 

CALL 

GRTXT ( ’ NF (GHZ)  ’  ,0) 

6570> 

CALL 

GRPRNT (  +  1 ,2.1 ,-1 ,Y*.53, .1 ,0.) 

6580> 

CALL 

GRTXT ( ’MAJOR’ ,0) 

6590> 

CALL 

QRPRNT  < ♦ 1 ,2.1 . -1 ,  Y+ .38 , .1 ,0.1 

6600> 

CALL 

GRTXTI ’MINOR’ ,0) 

6  6 1 0  > 

CALL 

QRPRNT ( +1 ,2.1 ,-1 .Y+.23, .1 ,0.) 

6620> 

CALL 

GRTXTI ’ANGLE’ ,0) 

6630> 

CALL 

QRPRNT (+1 , 2. 1,-1, Y+. 08, .1,0.) 

6640> 

C 

6850> 

67 

CALL 

GRNUMI ’ F ’  , -1  , 3 , NF ) 

6660> 

CALL 

QRPRNT (0 , X , -1 ,  Y* . 63 , .10,0.) 

6670> 

CALL 

GRNUMI ’F’ ,-1 , 2 , FMAX) 

6680> 

CALL 

QRPRNT 1 0 , X , -1 , Y+ .38 , .1 ,0. ) 

6690> 

CALL 

GRNUMI ’ F ’ ,-1 , 2.FMIN) 

6700> 

CALL 

QRPRNT (0 ,X , -1 ,  Y  +  .  23 , .1,0.) 

0  7 1 0  > 

CALL 

GRNUMI’F’ ,-1,2, FANG ) 

6720> 

CALL 

GRPRNT(0,X,-1 , Y+ . 08 , .1 ,0.) 

6730> 

C 

6740> 

CALL 

QRCP I ' POLAR' ) 

6750> 

Y=GSIZEI2) 

6760> 

CALL 

GRP  1C (X.Y.X  +  1 . , Y ♦ 1 .) 

6770> 

CALL 

GRSCLIO. ,0. , (SCALE (3) -SCALE  11)1/6.5, 360.) 

6780> 

CALL 

GRPICIX-1 . .Y-.28.X+1 . , Y* . 5) 

6790> 

CALL 

GRW INDI-IDEL.NDI1 ,2) ,ND(1 ,1) , 1 DUM , 0 , DUM , DUM , IDUM) 

6800> 

CALL 

GRUPPN 

681 0> 

CALL 

QRRES 1  SAVE ) 

6820> 

C 

6830> 

C 

6840> 

45 

ANG=ANG+DELTA  (3 ) 

6850> 

38 

CONTINUE 

6860> 

CLOSE ( UN IT=IOL) 

6870> 

CLOSE IUNIT= IOT) 

6880> 

CALL 

EXIT 

6890> 

END 
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APPENDIX 

A  PLOTTING  COORDINATE  SYSTEM  FOR  THE  APL  RCS/STATISTICS  CODE 


Background 

The  following  remarks  are  applicable  to  the  APL 
RCS/Statistics  Code  described  in  this  document.  In 
consequence  of  the  code’s  sequential  development 
over  a  rather  long  time,  several  peculiarities  exist.  Be¬ 
cause  they  are  not  explicitly  pointed  out  in  the  body 
of  the  text,  their  significance  may  be  overlooked, 
particularly  in  regard  to  recent  demands  on  the  code 
that  were  not  anticipated  during  its  development. 

The  pecularities  have  mainly  to  do  with  coordinate 
systems.  In  early  work,  a  master  coordinate  system 
(using  unprimed  variables)  was  chosen  arbitrarily 
with  the  z  axis  pointing  forward,  the  x  axis  to  port, 
and  the  y  axis  up.  That  master  system  was  connected 
to  the  component  coordinate  systems  by  transforma¬ 
tions  called  TRANS  in  the  code. 

Since  a  comparison  of  theoretical  calculations  with 
experimental  data  is  always  desirable,  and  since  the 
RATSCAT  facility  is  perhaps  the  chief  supplier  of 
such  data,  it  became  advisable  to  introduce  the  RAT- 
SCAT  standard  into  our  master  system.  RATSCAT 
uses  the  convention  z  axis  forward,  x  axis  up,  and  y 
axis  pointing  to  starboard.  Rather  than  change  the 
code  structure,  the  RATSCAT  master  system  was 
">nnected  to  our  master  system  by  a  simple  linear 
transformation.  The  shadowing  of  one  component 
by  another  was  now  defined  in  the  RATSCAT  system 
(the  instructions  are  called  VALID  in  the  code).  Since 
the  target  aspect  modifications  of  pitch,  roll,  yaw, 
and  tilt  are  also  available  in  the  RATSCAT  system  it 
was  decided  to  introduce  these  transformations  into 
the  code. 

Statistical  information  about  the  lobing  structure 
of  an  RCS  pattern  can  be  regained  provided  that  the 
component  scatterer  locations  can  be  specified. 
When  this  development  was  completed,  the  locations 
were  specified  by  giving  the  position  of  the  origin  of 
each  component  system  in  the  RATSCAT  system 
(this  instruction  is  called  ORGN  in  the  code). 

Now  we  have  the  requirement  to  describe  RCS 
over  an  arbitrary  planat  cut  or  over  a  conical  cut 
about  an  arbitrary  cone  axis. 


Modifications 

Because  shadowing  and  component  origins  are 
specified  in  the  RATSCAT,  or  radar-fixed,  coor¬ 


dinate  system,  the  pitch,  roll,  yaw,  and  tilt  transfor¬ 
mations  are  disallowed  for  other  than  small  angular 
excursions.  But  they  are  exactly  the  transformations 
one  would  like  to  use  to  satisfy  the  requirement  for 
arbitrary  planar  and  conical  cuts.  The  best  way  to 
satisfy  the  requirement  is  to  eliminate  the  original 
pitch,  roll,  yaw,  and  tilt  capability  and  replace  it  with 
a  new  system  called  PLOTTING.  The  target-fixed 
system  stands  fast  while  the  plotting  system  axes 
(originally  coincident  with  the  radar-fixed  axes)  rre 
“pitched,”  “rolled,”  or  “yawed”  to  a  new  posi¬ 
tion.  Thus,  the  original  terminology  is  retained  but 
with  the  meaning  just  defined.  (Actually,  one  can 
think  of  the  target  itself  being  pitched,  rolled,  or 
yawed,  but  in  an  angular  direction  opposite  to  that  of 
the  axes  movement.) 

In  addition  to  this  change,  the  original  unprimed 
coordinate  system  is  eliminated.  Henceforth, 
TRANS  expresses  a  transformation  from  the  primed 
to  the  RATSCAT  system  (in  a  manner  exactly 
analogous  to  the  way  TRANS  previously  expressed  a 
transformation  from  the  primed  to  the  unprimed 
system). 


Transformations 

There  are,  therefore,  three  systems  to  deal  with: 
the  plotting  or  display  system,  D ;  the  target-fixed  or 
RATSCAT  system,  R;  and  the  primed  system  of  a 
component,  P. 

D  is  connected  with  R  by  a  linear  transformation, 

C: 


Y/>  —  Gl-Gf  +  +  CjjZfl 


y n  ~  Gi  *r  "b  G2-^k  y  Gj  Zr 


Zd  ~  Gu -*«  "b  G’-Vr  b  C„  iff 


where  the  caret  denotes  unit  vector. 

The  user  specifies  a  region  given  by  angular  inter¬ 
vals  for  the  display  polar  and  azimuth  angles  6D  and 
<t>D,  respectively.  For  them,  he  must  calculate  the 
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corresponding  angles  in  the  R  system  (where 
shadowing  is  defined)  as  follows: 

cos  dR  =  c,,  sin  $n  cos  4>n  +  c23  sin  0D  sin 

+  c, 3  cos  6n 


and  U  systems  (in  the  same  order  as  R  and  D  are 
connected  above).  It  is  now  identical  with  the  process 
that,  henceforth,  connects  the  P  and  R  systems. 
Furthermore,  since  the  U  system  is  now  eliminated, 
the  transformation. 


cos  </>„ 

_  c*, ,  sin  0n  cos  <j>n  +  c\,  sin  dD  sin  <j>,}  +  c3l  cos  dp 
sin  dR 


B  = 


0  1  0 

-10  0 
0  0  1 


(note  the  subscript  R  in  the  denominator  above) 
sin  <t>R 

_  c |2  sin  dD  cos  <t>p  +  c22  sin  sin  <t>D  +  c32  cos  dn 
sin  dR 

(and  here  too). 

If 

cn  sin  dD  cos  <t>p  +  c2l  sin  dn  sin  4>n 
+  c3,  cos  dD  =  0, 


must  also  be  eliminated. 

To  recapitulate,  in  the  revised  arrangement  there 
are  only  the  three  coordinate  systems:  the  display 
system,  D ;  the  target-fixed,  or  RATSCAT,  system, 
R;  and  the  primed,  or  component-fixed,  system,  P. 
Henceforth,  the  three  instructions  TRANS,  VALID, 
and  ORGN  will  be  referred  to  the  R  system. 

Pitch,  Roil,  Yaw 

The  linear  transformation  C  can  be  expressed  as 
any  one  of,  or  a  combination  of,  pitch,  roll,  or  yaw, 
denoted  by  matrixes  P,  R,  and  Y,  respectively. 
Taking  the  original  position  of  D  to  be  coincident 
with  R  and  moving  D  while  keeping  R  fixed,  the 
angles  are  defined  as  follows: 


then  <t>R  =  ±  x/2  when 


Cjj  sin  dn  cos  <t>D  +  c22  sin  dn  sin  <t>n  +  c,2  cos  dn 


Pitch  through  an  angle  p  is  a  rotation  of  D  about 
yK,  positive  from  zH  toward  xK.  C  takes  the  special 
form  P,  where 


is  qreater  than  or  less  than  zero,  respectively. 
If 


P  = 


cos  p 
0 

sinp 


0 

1 

0 


-  sin  p 
0 

cos  p 


c,2  sin  dp  cos  0„  +  c2:  sin  dn  sin  Cp 
+  c,2  cos  6„  =  0, 


Roll  through  an  angle  r  is  a  rotation  of  D  about  zH, 
positive  from  xH  toward  yF.  C  takes  the  special  form 
R  where 


then  4>f  =  0  or  ir  when 


cn  sin  dp  cos  4>n  +  c2]  sin  dD  sin  <£„ 


R  = 


cos  r 
-  sin  r 
0 


sin  r 
cos  r 
0 


0 

0 

1 


+  Cj,  COS  dp 


is  greater  than  or  less  than  zero,  respectively. 

If  dR  =  0,  then  0„  is  indeterminate  because  the 
direction  of  interest  is  identical  with  z„  ■ 

As  was  indicated  above,  these  formulas  represent 
the  same  process  used  previously  to  connect  the  P 


Yaw  through  an  angle  y  is  a  rotation  of  D  about 
xF,  positive  from  yH  toward  zK.  C  takes  the  special 
form  Y,  where 


Y 


76 


1 

0 

0 


0  0 
cos  y  sin  y 

-sinj  cos>> 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL  MaRYLANO 


When  applied  successively,  P,  R,  and  Y  each 
represents,  in  some  order  (these  transformations  are 
not  commutative),  a  rotation  of  the  coordinate 
system  from  its  last  position  to  a  new  position.  The 
complete  transformation  C,  equal  to  a  product  of  P, 
R,  and  Y  in  some  order,  represents  a  transformation 
from  the  R  system  to  the  D  system. 


Note  that 


PR 


0 

0 

-  1 


0-1 
1  0 

0  0 


*  RP. 


Examples 

Consider  a  90°  pitch  of  axes,  followed  by  a  180° 
roll,  with  zero  yaw.  Then 


-1 

0 

0 

0 

0 

- 1 

0 

0 

1 ' 

C  =  RP  = 

0 

-  1 

0 

0 

1 

0 

= 

0 

- 1 

0 

0 

0 

1 

1 

0 

0 

1 

0 

0 

For  example,  with  the  conditions  0  =  60°,  0°  <  <t> 
<  180°,  the  transformation  RP  produces  a  portside 
viewing,  from  front  to  rear,  over  a  conical  cut  at  a 
30°  elevation  above  the  horizontal  plane.  In  contrast, 
the  transformation  PR  produces  a  portside  viewing, 
from  rear  to  front,  over  a  conical  cut  at  a  30° 
declination  below  the  horizontal  plane. 
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